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The basis of the method is selection of a "performance index" which compares the 
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performance index for analytic feasibility and practical utility. In order to provide a 
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multivariable controllers, using the Dynamic Simulation Language. The resulting 
optimum values stemming from the analysis are the controller gains which minimize 
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THESIS DISCLAIMER 



1 



The reader is cautioned that computer programs developed in this research ma3/ 
not have been exercised for all cases of interest. While every effort has been made, 
within the time available, to ensure that the programs are free of computational and 
logic errors, they cannot be considered validated. Any application of these programs 
without additional verification is at the risk of the user. 
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I. INTRODUCTION 



The design of controllers has been an engineering endeavor since the 1920's. The 
road to success was significantly changed by the introduction of the International 
Business Machine Model 360 computer in the 1960's, which allowed the control 
designer to develop numerical models for various systems and efficiently test his design 
before prototyping. The development of FORTRAN programs such as the 
Continuous System Modeling Program (CSMP), and Dynamic Simulation Language 
(DSL) greatly simplified the modeling process and allowed the designer to have highly 
accurate numerical tools at hand, rather than having to rely on dosed form or 
reduced-matrix type solutions such as the Riccatti equation. 

Along similar lines the engineering discipline of optimization has developed. 
Numerical optimization, using mainframe and micro computers, has rapidly developed 
over the past ten years with the work of G. N. Vanderplaats and others. His 
FORTRAN programs Control Program for Engineering Synthesis (COPES) and most 
recently Automated Design Synthesis (ADS) have provided routines that implement all 
of the known methods available to do numerical optimization. 

The goal of optimization has been always to provide the best for the least, and 
this credo ean be applied to the design . of all controllers. One method of accomplishing 
this task will be demonstrated by the author. It involves the modeling of single input, 
single output (SISO) or multiple input, multiple output (MIMO) systems using DSL, 
specifying an applicable performance index by which to judge the system, defining 
constraints which specify the desired system performance characteristics as well as 
verifying system stablity, and then calling upon ADS. Through iterations the 
optimization routine finds the necessary values for system gains (defined as design 
variables) which will minimize the performance index (defined as the objective 
function). The final result being the "best" controller design. 
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II. MODERN CONTROLLER DESIGN 



A. ANALYSIS TECHNIQUES 

Two general methods of controller design have developed, each having certain 
advantages, depending upon the application of the controller involved. The first is 
frequency domain analysis, in which the response of the desired system to a sinusoidal 
input is examined. The other method is time domain analysis, where the system 
response to various time dependent inputs such as impulse, step, or ramp functions are 
investigated using a state-space representation of the system. This second method 
lends itself more easily to numerical modeling techniques that can be programmed on a 
computer, and will be the method of choice in this investigation. 

B. NUMERICAL .ANALYSIS TECHNIQUES 

1. eSM? 

The CSMP-III program was developed by International Business .Vlachine 
Corporation (IBM) to provide users with a simplified method for numerically solving a 
system of ordinary differential equations and related mathematical problems. [Ref l:p. 
2]. Today it is a widely used and respected tool for the design of control systems. 

Some drawbacks to CSMP were found when the initial considerations were 
made on how to perform numerical controller optimization. The most significant areas 
of concern were that CSMP uses the FORTRAN/HX complier, it is only available 
locally in the Multiple Virtual System (.MVS) environment, it uses single precision 
variables throughout, and that graphics are available only after processing through the 
DISSPLA graphics system. These facts led to a decision to select a different method to 
do controller optimization. 

2. DSL/VM 

DSL/VM is a recent product from IBM [Ref 2:p. iii]. The program is similar 
to CSMP in that simplified source code is translated into standard FORTRAN 
statements (as output to the file FORT FORTRAN), compiled and then run in order 
to perform the desired system simulation. Following the successful run, tabled output 
is directed to the user along with any desired graphic displays. Fig. 2.1 is a flow'chart 
of this procedure. 
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Significant improvements from CSMP to DSL include use of the 
FORTIL/XN/VS compiler, Virtual Machine (VM) or MVS environment availability, on 
line graphics capability, double precision calculations throughout, and automatic 
linkage to user defined libraries. Additionally, a one time job initialization segment is 
available which is separate from the model initial segment. Thus the programmer can 
easily rerun the time dependent simulation for the the same model, after changing any 
desired parameters, within the same job specification. These facts were significant 
reasons for selecting DSL as the modeling tool to use for controller optimization. 
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III. OPTIMIZATION 



A. ANALYSIS TECHNIQUES 

1. Basic iVIethcds 

The determinaiion of the "best" value for any design requires that several 
items be considered. One is the resource that can be varied, made up of what are called 
the design variables. Also included are the mathematical relationships that describe 
how the resource relates to the design environment. These relationships define what 
are called the objective function and constraint functions. .Many diverse methods have 
been developed to determine the values of the design variables that cause the objective 
function to be at its optimum value. With the advent of the computer, nunterical 
computation methods have become the most practical ones for treatment of complex 
systems, supplanting graphical and analytic detenninations. 

2. Numerical Methods 

Many different approaches have been developed for solving the general 
optimization problem by numerical methods. Numerical Optimization Techniques for 
Engineering Design [Ref. 3] describes most of the common analytic techniques and the 
corresponding algorithms to perform this task using a computer. A major extension of 
this work was the development of the Control Program for Engineering Synthesis 
(COPES) and Constrained Minimization (CON.MIN) FORTR.^\N routines that 
accomplish the various numerical tasks necessary to successfully determine an optimal 
solution to problems that can be appropriately defined. 

a. COPES and CONMIN Programs 

COPES is a FORTRy\N program that was developed as a method for 
easily coding the parameters required to define an optimization problem for solution by 
a computer. Its use requires that a FORTR.^N subroutine called ANALIZ be written 
which evaluates the objective function and all constraint functions. Previously, the 
FORTRAN subroutines grouped under the name CON.MIN were then called to 
conduct the desired optimization analysis. Recent improvements to the CONMIN 
routines have been generated and the package has been renamed Automated Design 
Synthesis Program (ADS) [Ref 4]. To use COPES with ADS, the routine COPES A 
should be called. 
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A major drawback to the use of COPESA is the specific formatting and 
subroutine requirements. This hinders interactive optimization, and has resulted in the 
selection of an alternative method for performing controller optimization. 
b, ADS Program 

ADS can be called directly from any FORTRAN program [Ref 4:pp. .5-81. 
Prior to making the call, all input parameters must be defined, and a looping logic 
must be provided in the calling program that causes the objective and constraint 
functions to be evaluated upon the completion of each optimization iteration, as 
depicted in Fig. 3.1 . 

With this simple logic many engineering problems can benefit from the 
advanced capabilities offered by the ADS program. The specific area that has served 
as a focus for this study is the interfacing of controller design, done with DSL, with 
optimization using ADS. 
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INITIALIZE PARAMETERS 
(ADS, DSL, and model) 



s- INITIAL CALL TO ADS 

(for problem setup) 

I 

EVALUATE OBJECTIVE AND CONSTRAINTS 
(based upon design variable values) 

u 

CALL ADS 

t 

EVALUATE GRADIENTS 
(constraints and objective) 



DETERMINE NEW VALUES 
OF DESIGN VARIABLES 



Y 

RETURN TO DSL 
(conduct model simulation) 



Y 

CHECK FOR ADS COMPLETE 



OPTIMIZATION COMPLETE 



y 

GR.^PHICS POST-PROCESSING 



Figure 3.1 Program Logic Calling ADS Directly 
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IV. AUTOMATED DESIGN SYNTHESIS/DYNAMIC SIMULATION 

LANGUAGE 



A. OVERVIEW 

Linking of the Automated Design Synthesis Program and the Dynamic 
Simulation Language has been successfully demonstrated, and the acronym ADSL 
represents the close coupling of these two FORTRAN programs.- In the following 
discussion the methodology of how these two programs interface will be explained. 

B. PROGRAM SEGMENTS 

Figure 4.1 is a flowchart of how the individual segments of a DSL program are 
utilized with the ADS optimization code. The DSL Language Reference Manual 
[Ref 2] is the source document on the syntax and logic of the DSL language. Chapter 
3 of the manual explains in detail the function of each phase of DSL program 
processing. For the proposed method of controller design each phase is utilized. 

1. Program Initialization (INITLZ) - 

This portion of a DSL program is the first executed segment of code. In it 
arrays are dimensioned, non-standard variable types (fixed, complex, etc.) are defined, 
and initial values of variables are assigned. Additionally, DSL parameters which 
specify the output variables, the desired integration method, and the technique for 
plotting can be specified. 

The important fact about this section is that it is only executed one time. 
Thus the initial values that ADS requires about the status of the optimization, denoted 
by the variable INFO, and the array dimensioning can be done in this section. 

2. Run Initialization (INITIAL) 

The beginning of each system response simulation starts at this point with a 
call to ADS. On the initial call, ADS is provided with best-guess values for the design 
variables, the upper and lower bounds for each design variable, along with the number 
and types of constraints and their allowable ranges. ADS will return new values for 
the design variables to be used in the controller design. On subsequent calls to ADS 
updated values of the objective and constraint functions will be provided to ADS for 
its determination of the optimal values of the design variables. 
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The remainder of this program section sets various constants at the initial 
values that will be needed later in the dynamic section. 

3. Derivative Section (DERIVATIVE) 

This section is the heart of the simulation.. The inputs to the system are 
defined, the system transfer function is specified, and various outputs, states and the 
objective function are evaluated. Great latitude is possible in each of these 
specifications, and for this reason DSL is an outstanding method for analysis of 
controller design. 

This section also determines the integration time step that the simulation will 
take, based upon two criteria. The first is how much time is remaining before a 
specified data output time is reached, and the second is how big an integration step can 
be taken while maintaining the specified relative and absolute errors. The output time 
is determined through the use of the DELPRT and DELPLT commands. These 
control parameters indicate the desired time intervals for saving intermediate values of 
specified variables for both plotting and listing purposes. 

4. Dynamic Section (DYNAMIC) 

Once the derivative section finalizes what the size of the next integration 
interval will be, it then calculates the magnitude of the various system states, outputs, 
and objective function at the end of the integration interval. This information is then 
passed onto the dynamic section where it can be evaluated to determine the various 
constraint values and whether or not the simulation is complete (specific information 
on the evaluation of these parameters will be left to subsequent discussions of the 
single input, .single output (SISO) example). If the simulation is incomplete, 
computational flow is is directed back to the derivative secton. 

When the completion of the simulation run is determined, by either reaching a 
specified slope of the objective function, or a maximum time limit, the DSL routine 
ENDRUN is called. This procedure completes all outstanding calculations and then 
transfers control to the terminal segment. 

5. Terminal Section 

In this last program segment the ADS status variable INFO is checked to 
determine whether or not the optimization is complete. If it is, the DSL program is 
terminated by calling the routine ENDJOB. If further optimization is required, control 
is transferred back to the initial segment. 
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INITIALIZE ADS AND DSL PARAMETERS 

(INTTLZ) 



t 



CALL ADS TO PROVIDE 
VALUE OF DESIGN VARIABLE 
(INITIAL) 



EVALUATE INTEGRATION STEP 
(DERIVATIVE) 



\f 

EVALUATE CONSTRAINTS 
CHECK FOR STEADY STATE 
OBJECTIVE 
(DYNAMIC) 



V 

CHECK FOR ADS COMPLETE 
(TERMINAL) 



END RUN 



Figure 4.1 ADSL Program Segments 
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V. SINGLE INPUT SINGLE OUTPUT PROBLEM ANALYSIS 



A. CASE 1 SECOND ORDER SYSTEM 

In order to obtain a baseline for the ADSL program the well behaved, general, 
second order system was taken for examination. Ogata [Ref. 5] completely describes 
this simple system shown in Fig. 5.1 and its responses to various inputs. For simplicity 
this system will be called Case 1. In the present study, the input is taken to be the unit 
step function, R(s)= 1.0 for time greater than zero, and the magnitude O) was set equal 
to 1.0. 




Various analytic optimization methods have also been developed [Ref 5:pp. 296-306] 
for this system that can be used to verify the ADSL program. 

Objectives of this baseline investigation were as follows; 

• Determine the most efficient and correct integration method, 
o Determine the most efficient optimization strategy. 

• Determine the variations to be exoected from using different performance 
indices. 

• Determine the most overall-efficient method to obtain optimum design variable 
values. 

A sample program to accomplish these tasks is listed in Appendix C. 
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B. PROBLEM FORMULATION 



1. Objective Function 

The fundamental part of any optimization problem is definition of an 
objective function. For this problem the performance index is the logical relationship 
that defines how "good" a system is. Ogata (Ref. 5:p. 296], defines the various 
performance indices that have been used in this investigation, and these will be 
described in the following discussion. 

2. Design variable 

The design variable selected for this investigation was the damping ratio ( )• 

It is a reasonable parameter because it directly effects how the system responds to any 
input over time. Thus if is increased, the system output response will be slower and a 
larger accumulated error between the actual to desired output may be generated with 
respect to time. This larger error could then increase the value of the performance 
index/objective function. On the other hand, insufficient damping will lead to excessive 
output oscillations which also tend to increase the objective function. 

C APPLICATION OF CONSTRAINTS ' 

Ogata [Ref 5:p. 232] defines transient response specifications for the typical 
second order system. In almost all cases, these standard definitions can be directly 
applied to the response of a DSL model. A brief review of these specifications and 
how they are applied are discussed below. 

1. Maximum Overshoot 

How much the output response of an under-damped system goes beyond that 
which is desired, is the maximum overshoot. The larger this value is, the more energy 
is wasted, as a restoring force must be applied to return the system to the desired state. 
In addition, physical systems may be limited in the amount of amplitude overshoot 
that can be tolerated. A simple constraint of the maximum allowable amplitude, as a 
fraction of the intended steady state, can be defined. 

2. Delay Time 

The time the system requires to reach fifty percent of its intended response is 
the system delay time. It indicates how fast the system begins to respond, and a 
constraint of a maximum time can be established. 
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3. Rise Time 

When the system initially reaches 100 percent of the desired response, the rise 
time has been reached. Like the delay time, this parameter is also an indication of how 
quickly the system will respond. For under-damped systems, a constraint of maximum 
rise time can be programmed. 

4. Peak Time 

The time of maximum response is called the peak time. It also shows how 
fast the system responds. In cases where a peak time occurs, a maximum peak time 
constraint can be specified. 

5. Settling Time/Completion Time 

Settling time is defined [Ref 5:p. 236] as when the response falls and stag's 
within a prescribed percentage (typically two or five percent) of the total response. For 
computational considerations, the settling time was redefined for the two possible 
response conditions: 

1) Underdamped case: that time when the absolute maximum of a response 

oscillation peak is within one percent of the total response transient. 

2) Overdamped case: that time, greater than five time constants, at which the 
absolute system response is within one percent of the total response. A time 
constant is defined as that time required for the system output response to 
reach sixty three percent of the final response value. 

In order to keep the above definition separate from settling time, the new consideration 
will be called completion time. 

The different constraint criteria are illustrated in Fig. 5.2, with the constraints 
of maximum amplitude, rise time, delay time, peak time, and completion time, being 
designated by (1) to (5) repectively. 

D. INTEGRATION CONSIDERATIONS 

DSL provides nine optional integration methods which can be grouped into the 
following types: 

• Fixed-step 

• Variable-step 

• Variable-step, Variable-order. 

The default method is the Runge-Kutta Fifth Order variable-step method. Selection of 
the method which is best suited to the controller design problem is discussed below. 
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Figure 5.2 Case 1 Demonstration oF Constraints 
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1. Limits of Integration Determination 

For this analysis the objective function has been defined as the system 
performance index. The particular performance indices to be examined are weighted 
time integral errors between the actual and desired output response of the system from 
the time the input is applied to infinity. Noting that infinity is an impractical number 
for a commuter (or anything else, for that matter), a reasonable approximation must be 
found. Two criteria were considered for determining what cut-off time to use to mark 
the end of calculations, or time equivalent to infinity. 

One cut-off time criteria would be to simply select the completion time, but 
this has a significant drawback. Different values of the design variable will cause 
different cut-off times, and thus this method does not provide a consistent manner for 
evaluating the objective function for different designs. A second method is to evaluate 
the rate of growth of the objective function with time. For a stable system the 
objective function initially increases rapidly due to the mis-match of input and output. 
As time progresses this error approaches zero and the rate of change of the objective 
function slows dowm. Once the completion time is reached, if the rate of change of the 
objective function is evaluated to be essentially zero, a consistent method of 
determining the cu^-off time can be established for all values of the design variable. 

2. Integration Method Analysis 

Typical output and objective functions are shown in Figs. 5.2 and 5.3. These 
rapidly increase from zero and approach a constant value asymptotically. This type of 
behavior is described by Speckhard [Ref. l:p. 268], and is characteristic of "stiff' 
systems. In considering which integration method would best deal with this stiff 
system problem, fixed-step methods were ruled out due to their inability to accurately 
deal with rapidly changing functions. The Runge-Kutta Fifth Order (RK5), Adams 
variable-order (ADAMS), and Gear Full Jacobian (STIFF) variable-step methods were 
compared to determine which would result in the maximum computational efliciency. 

The computational efficiency in this problem can be considered in two parts; 

1) The number of individual integration steps 

2) The number of objective function evaluations requested by the optimization 

routines for gradient determination. 

The conclusion is then to use an integration method that maximizes the step size and 
an optimization routine that manimizes the number of function evaluations. The DSL 
Reference Manual [Ref 2:p. 3.16] describes the reasons for differences in execution time 
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Figure 5,3 Case I Objective Function 
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between the integration methods, and points out that for oscillatory solutions, such as 
those seen for the second order underdamped system response, all three of the methods 
are essentially equal,. Consequently, the best method to use is the one that generates 
the least number of objective function evaluation requests by the optimizer. 

Three ADSL programs were generated, each using a different integration 
method, but using the same input, system, optimization method, and objective function 
(performance index). Discussion of optimization methods and performance indices will 
follow. In the problem at hand the optimization method was to use an exterior 
penalty function strategy with the Broydon - Fletcher - Golfarb - Shanno optimizer 
and polynomial interpolation for the one dimensional search. For the performance 
index the integral square error (ISE) was employed. The results are compared in Table 
1. The applied constraints, designated by (1) to (5) rcpectively in Fig. 5.2 and Table 1, 
were: 

(1) maximum amplitude 

(2) rise time 

(3) delay time 

(4) peak time 

(5) completion time 

Interesting points from this investigation include: 

• All three integration methods found the same result within the default tolerance 
of ADS. 

® Each method, starting with the same initial design variable value, calculated a 
different initial objective function value. 

• The number of objective function calls for each method was almost the same, 
but the STIFF method, as one might have predicted, required the fewest. 

E. OPTIMIZATION METHODS 

In order to use ADS, an optimizer and one-dimensional search must be selected. 
Optionally, a strategy may be invoked. 

1. Optimization Strategy 

.An optimization strategy is simply the algorithm by which the overall 
optimization problem is transformed into one in which constraints are weighted or 
eliminated and the objective function may be redefined by use of a "pseudo-objective" 
function. Although use of a strategy is not required, some problems have shown 
significant computational improvement as a consequence. 
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TABLE 1 



CASE 1 INTEGR.ATION METHOD RESULTS 



Method 


RK5 


ADAMS 


STIFF 


Initial Design 


1.0 


1.0 


1.0 


Initial Objective 


4008.9 


12334. 


5073.9 


Final Design 


0.717 


0.718 


0.720 


Final Objective 


1.968 


1.975 


1.965 


Constraint 1 


-.361 


-.361 


-.362 


Constraint 2 


-.157 


-.190 


-.157 


Constraint 3 


-.121 


-.103 


-.077 


Constraint 4 


-.170 


-.113 


-.016 


Constraint 5 


-2.24 


-2.55 


-2.69 


Number of Runs 


49 


46 


42 



i 

I 



a. Sequential Unconstrained Minimization Techniques (SUMT) 

SUMT involves the creation of a pseudo-objective function that is the sum 
of the original objective function and a penalty function. This penalty function may be 
one of four types: 

1) Exterior Penalty 

2) Linear Extended Interior 

3) Quadratic Extended Interior 

4) Cubic Extended Interior 

All these methods cause the objective function to have its actual value when the design 
variable is within the feasible region, but when the design variable is outside the 
feasible region, the penalty function significantly increases the objective function value. 
The penalty functions listed above vary in the value that will be added to the objective 
and how the transition between the feasible and infeasible regions is conducted 
[Ref 3;p. 136]. 

b. Augmented Lagrange Multiplier {ALM) 

In this method a pseudo-objective function is generated based upon the 
exterior penalty function, but W'ith an additional term in order to create the 
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Lagrangian. It can be shown that the solution to the Lagrangian is the solution of the 
original optimization problem [Ref. 3:p. 141]. 

c. Sequential Linear Programming {SLP) 

The basis of this strategy is the simple linearization of the objective 
function and constraints. Rectangular move limits are established and then linear 
programming techniques are applied to optimize the objective function [Ref 3:p. 155]. 

d. Method of Centers 

The SLP method can produce infeasible designs with each iteration. To 
correct this problem the Method of Centers generates an inscribed hypersphere based 
upon the linearized objective function and constraints, then uses linear programming 
techniques to determine the optimum. 

e. Sequential Programming 

Sequential Programming can be of either of two types; 

1) Sequential Quadratic 

2) Sequential Convex 

Both of 'hese methods refine the method for determining the search direction for the 
objective function minimum. In the quadratic method the objective function is 
transformed into a quadratic function, while the convex method generates the 
reciprocal of the objective function. These methods create conservative 
approximations to the actual optimization problem which can then be easily evaluated 
[Ref 4:p. 2]. 

2. Optimizers 

The method by which new design variables are chosen is a result of the 
optimizer selected. Each of the possible methods predicts a new set of design variables 
which should be closer to the optimum design, based upon evaluations of the objective 
function at previous values of the design variables. 

a. Conjugate Direction Method 

This optimizer, developed by Fletcher and Reeves [Ref 6;pp. 149-154], is 
based upon finding the negative gradient of the objective function from some initial 
point and then performing a one-dimensional search in that direction for the minimum 
value of the objective function. This process is repeated from the newly found point, 
and then a new direction is determined based on' the previous two points and 
directions. 
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TABLE 2 
ADS OPTIONS 



Strategy (ISTRAT) 

0 - None 

1 - SUMT, Exterior Penalty Function 

2 - SNMT, Linear Extended Interior 

3 - SUMT, Quadratic Extended Interior 

4 - Cubic Extended Interior 

5 - Augmented Lagrange Multiplier Method 

6 - Sequential Linear Progranuning 

7 - Method of Centers 

I 8 - Sequential Quadratic Programming 
! 9 - Sequential Convex Programming 

Optimizer (I OPT) 

1 - Fletcher- Reeves 

2 - Davidon-Fletcher-Powell (DFP) 

3 - Broydon-Fletcher-Golfarb-Shanno (BFGS) 

4 - Method of Feasible Directions 

5 - Modified Method of Feasible Directions 

One dimensional Search (lONED) 

1 - Golden Section Method 

2 - Golden Section and Polynomial 

3 - Polynomial Interpolation, bounded 

4 - Polynomial Extrapolation 

5 - Golden Section Method 

6 - Golden Section and Polynomial 

7 - Polynomial Interpolation, bounded 

8 - Polynomial Extrapolation 

Note: One dimensional searches 1 through 4 are for use with 
Optimizers I through 3. and searches 5 through 8 are for 
use with Optimizers 4 and 5. 
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Description Manual [Ref. 4; pp.3-5]. In ADS, each possible combination is identified 
by use of a three digit code (XYZ) with the X digit representing the optimization 
strategy, the Y digit the optimizer, and the Z digit the one-dimensional search 
employed. 

Vanderplaats [Ref 4] points out that although there is a multitude of different 
combinations, eight combinations have, as a result of experience, become 
recommended for general applications. These eight combinations were tested using 
otherwise identical ADSL programs. The results are presented in Table 3. 



TABLE 3 

CASE 1 ADS STIL^TEGY RESULTS 



Strategy 




0^7 


057 


857 


957 


533 


233 


133 


656 


Objective 




1.966 


1.966 


1.985 


1.959 


1.964 


1.964 


1.964 


1.972 


Design Variable 


0.720 


0.719 


0.703 


0.729 


0.722 


0.722 


0.722 


0.714 


Constraint 


1 


-.362 


-.361 


-.355 


-.365 


-.362 


-.362 


-.362 


-.359 


Constraint 


2 


-.157 


-.157 


-.150 


-. 05 ^ 


-.158 


-.158 


-.158 


-.154 


Constraint 


3 


-.077 


t .089 


-.207 


-.008 


-.057 


-. 069 ^ 


-.069 


-.126 


Constraint 




-.009 


0.000 


-.036 


+ . 0^8 


-.018 


-.014 


-.014 


+ .033 


Constraint 


5 


- 2.69 


- 2 . 8 ^ 


- 2.95 


- 2.49 


- 2.64 


- 2.61 


- 2.61 


- 2.89 


Number of Runs 


10 


27 


11 


10 


47 


51 


51 
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Observations that can be made are: 

• All the combinations found essentially the same optimum value. 

» There are significant differences in the efficiency of the various methods.. If the 
criteria of minimum number of objective function evaluations is used, the 047 
and 957 strategies are the best. 

• .All the combinations found solutions that met all the constraints within the 
ADS default tolerance band. 

® More objective function evaluations did not necessarily lead to reduction in 
objective function value. In fact, one of the quickest methods (957), found the 
smallest value. 
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F. PERFORMANCE INDEX SELECTION 

Ogata [Ref. 5: pp. 296-301]. presents four performance indices that are based 

upon the time-integrated weighted error between the desired output and the actual 

system response. These indices are: 

» Integral square-error (ISE) 

« Integral-of-time-multiplied square-error (ITSE) 

• Integral absolute-error (lAE) 

• Integral-of-time-multiplied absolute-error (ITAE) 

These performance indices var\’ in complexity for analytical evaluation over the 
integration period of zero to infinity. For two of the cases, namely ISE and ITSE, 

exact analytical answers can be found as indicated in Table 4. By using ADSL the 

analytic problems are greatly simplified for the other indices. 



The four performance indices were placed into identical ADSL programs as the 
objective function for optimization, and the results are shown in Table 4. Significant 
observations that can be drawn from this table are: 

« As expected, each of the methods found a different optimum value for the 
design variable. 

• All of the methods found an optimum value that did not violate any constraint. 



TABLE 4 

CASE 1 PERFORMANCE INDEX RESULTS 



Index 

Strategy 

ADSL Obj-ective 
Exact Objective 
Design Variable { ^ ) 
Exact 

Constraint 1 
Constraint 2 
Constraint 3 
Constraint 4 
Constraint 5 

Number of Runs 



lAE ISE ITAE ITSE 

047 047 047 047 

1.605 1.014 .1.966 .707. 

n.a. 1.000 n.a. .707 

.6605 .5909 .7202 .5928 

n.a. .5000 n.a. .5946 

-.337 -.300 -.362 -.301 

-.227 -.152 -.157 -.207 

-.447 -.858 -.077 -.831 

-.299 -.717 -.009 -.646 

-3.32 -3.64 -2.69 -3.94 

23 40 10 45 
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• Within numerical tolerance, the results are the same as those found by Ogata 
• [Ref. 5; p. 300]. 

« The apparent discrepancy in optimum damping ratio for the ISE performance 
index is due to the flat minimum of the objective function in this case (a 
diiTerence of 20% in C gives a performance index variation of less than 2%). 

It should be noted that the magnitude of the resulting objective functions cannot 
be compared directly to determine which method is "better." This would be like 
comparing apples and oranges, due to the significant mathematical difierence in the 
objective functions/performance indices. 
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VI. INTEGRAL FEEDBACK CONTROL DESIGN 

A. ST/vTE SPACE ANALYSIS 

The ume and frequency domain methods discussed in the previous chapter are 
m.ost useful in the analysis of single input, single output (SISO) systems. State-space 
methods give added flexibility in controller design and allow the treatment of complex, 
multiple input, multiple output (MI MO) systems. The concept of state, which is the 
basis of modern control theory, embodies the following key elements; 

• State Variables: the smallest set of n variables, with their initial quantities 
known, that completely describes the behavior of a dynamic system for a given 
input over time. It should be noted that these states do not necessarily have to 
be physically accessible. 

• State; the condition of a dynamic system at a given time (t) that is uniquely 
determined by the initial state of the system and the inputs to the system after 
the initial time (t^). 

• State Vector: the vector made up of the n state variables needed to describe the 
system. The state vector is usually designated x(t). 

• Input Vector; the vector made up of the m variables that are inputs to the 
system. The input vector is usually designated u(t). 

• State Space: is the n dimensional space made, up of the coordinate axes 
corresponding the individual states. 

Kuo [Ref. 7] describes two types of basic systems that are typical of modern 
control systems: 

1) State feedback control 

2) State feedback with integral control 

These systems have been modeled using the combined Dynamic Simulation Language 
(DSL) and Automated Design Synthesis (ADS) optimization code ADSL. 

B. CASE 2 STATE FEEDBACK WITH INTEGRAL CONTROL 

Case 2 is the next problem to be investigated. In this system control is attained 
by feedback of the state variables through constant gains and an integrator. A signal 
flow graph representative of state feedback with integral control is shown in Fig. 6.1. 
As a mafer of interest, if the integrator is removed the resulting system is called a state 
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feedback controller. When the reference for any state feedback controller is set to zero, 
the resulting system is classified as a regulator. It can further be shown that the 
common proportional-integral-derivative (PID) and rate-feedback control systems are 
special cases of state feedback control [Ref 7:p. 522]. 

Case 2 systems are of interest because the state feedback systems are generally 
only useful as regulators that have no noise input. Integral control suppresses the 
fluctuations caused by undesirable noise and thus provides for a much more effective 
system. 

1. Model Description 

The system selected for optimization was a dc electric motor speed control 
system proposed by Kuo [Ref 7:pp. 532-536]. In this example the motor shaft angular 
velocity and armature current were selected as the state variables Xj and x^ 
respectively. Two inputs were considered: 

• Wj as a direct, torque loading acting on the motor shaft 

• W 2 as the desired constant speed set point value. 

The control and output variables can respectively be represented by the following 
relationships where g^, g 2 , and g^ are the feedback gains, and the state variables are x^ 
and X 2 : 

• u(t) = -gjXj - g2X2 - g3J c(t)dt 

• c(t) = W2 - Xj. 

The corresponding signal flow graph showing the various relations in the control 
system is shown in Fig. 6.1. The initial conditions and values assumed for the design 
are shown in Tables 5 through 8. The resulting state responses are displayed in Fig. 
6 . 2 . 

The goal of the optimization is to determine the values of the three gains - 
gp g 2 , and g^ — necessary to minimize a stated objective function without violating 
specified constraints for controller performance and stability. 

2. ADSL Program 

Appendix D contains the source code listing of the ADSL program which 
models the state feedback system with integral control for Case 2. It follows the same 
form of construction as used for Case 1. The initialization of parameters was expanded 
into three distinct groupings: 

1) Array dimensioning and parameter initialization for ADS 

2) DSL control parameter setting 
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STATE FEEDBACK WITH INTEGRAL CONTROL 



Figure 6. 1 Case 2 Signal Flow Graph 

3) Model description, constraint evaluation, and integration completion control 
parameter implementation. 

Using the results from Case 1, selection of the stiff integration method was made. An 
additional change was to reduce both the relative and absolute error criteria for the 
integration routine from l.OE-05 to l.OE-04. This was done in order to increase 
computational efficiency while still providing for an order of magnitude greater 
precision than that expected for the objective function optimization done by the ADS 
program, using the default parameter settings. 

3. Performance Index/Objective Function 

For optimizing, the classical quadratic performance index (QPI) was adopted. 
This performance index: J = qJ'^(x*Qx + u^Rii)dt allows the user to specify 

individual weights for the Q and R matrices. Thus the desired relative importance 
between the error seen in each of the actual state values to the desired state values, and 
for the energy use associated with the control variable over time can be specified. For 
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TABLE 5 

CASE 2 INITIAL CONDITIONS 

Motor Shaft Angular Velocity (xj ) = 0.0 rad/s 
Armature Current (x 2 ) = 0.0 Amps 
Feedback Control Signal (u(t)) = 0.0 Volts 
Input Disturbance Torque (wj ) = 0.0 N-m 
Input Reference Set Point Speed (w 2 ) = 1.0 rad/s 



TABLE 6 

CASE 2 INITIAL POLE SELECTION 

Initial Pole Selection 
-lO.O+lO.Oj -lO.O-lO.Oj -300.0 

Exact Resulting Initial Gains 
g| = -0.38 g2 = 0.60 g^ = -6.00 



Case 2 each state was selected to be independent of the other, and thus Q and R were 
identity matrices. Table 8 shows each of the matrices used in the Case 2 performance 
index. 

4. Exact Initial Gain Determination 

ADS requires that an initial guess for the design variables be provided in order 
to start the optimization. In conjunction with the the determination of these values a 
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TABLE 7 

CASE 2 CONSTRAINTS 

CS? = Maximum Value of the Output c(t) = 140yo of initial value 
TSP = Maximum Completion Time = 1.5 sec 
USP = Maximum Value of the Control u(t) = 10.0 V 



TABLE S 

CASE 2 PERFORMANCE INDEX MATRICES 

Q = 1.0 0.0 R = 1.0 

0.0 1.0 

X = W2-Xj u = -gjXi-g2X2-g3Jc(t)dt 

^2 



reasonable assumption has been made that the controller designer be able to specify 
stable poles for the system. Kuo [Ref 7:pp. 527-528] outlines an analytic, general 
method for the determination of initial state feedback gains without transformation to 
the phase variable canonical form of the system characteristic equation. The method 
requires defining the following terms: 

• Open loop characteristic equation: 

Aq = Det( 5 I - A ) 

• Closed loop characteristic equation with state feedback: 

= Det( .s I - A + BG ) 

• Relation Matrix: 

k (s) = [ Adj ( 5 I - A )] B 
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Where the indicated matrices are as follows for n states; 

A- the coefficient matrix for the state variables (n x n). 

B- the coefficient matrix for the control input (n x 1). 

G- the state feedback gain matrix ( 1 x n). 

I- the identity matrix. 

The relationship Gk (s) = A^(s) - Aq(s) can then be solved for G using the 
specified initial design system poles as the eigenvalues for the closed loop characteristic 
equation. In the Case 2 ADSL program, shown in Appendix D, the necessary matrix 
operations have been multiplied- out and the like terms relating the various powers of s 
have been grouped together to yield a system of three equations and three unknowns. 
Matrix solver routines from the LINPAC library' [Ref. Sj, DGEFA (double precision 
general matrix conditioning) and DGESL (double precision general matrix solver), were 
then called to respectively decompose and solve the system of equations, and thus 
generate the exact initial state feedback gains. 

0. Approximate Initial Gain Determination 

The matrix method just described above works well, and is an efficient method 
for systems of order three or less. With higher order systems the determination of 
adjoint matrices is difficult and time consuming. To counter this problem, an 
alternative method was devised. The determination of the gains was considered as an 
optimization problem, where the constraints and objective function were formulated 
such that the there would be a minimum difference between the poles corresponding to 
the gains determined by optimization using ADS, and the original poles specified by 
the designer. Two different forms of the optimization problem were developed. 

a. Determinant Method 

For the first method, the complex matrix: ( s I - .A ^ BG ) was formed by 
substituting in one of the desired poles for 5 and setting the initial values for the gains 
to 1.0. Tiie determinant of this matrix was then found using the LINP.ACK routines 
ZGEF.A (complex general matrix conditioning) and ZGEDI (complex general matrix 
derterminant). If the correct value of the gain, corresponding to the desired pole was 
substituted into the matrix equation, the resulting determinant would be equal to zero. 
Thus an equality constraint can be established equal to the absolute value of the 
resulting determinant. This process was then repeated for each of the other desired 
poles. The objective function was then defined as the sum of all the constraint values. 
For an exact solution the objective function would be equal to the value of zero. ADS 
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was then called to vary the gain values to find the minimum value of the objective 
function. The final determined gains thus should represent an approximate match to 
the desired initial ooles. 

When this algorithm was attempted, the results were not as expected. 
Tjbie 9 shews the results for various initial poles using the Case 2 model. Looking at 
possible reasons for these outcomes it can be found that the optimization problem, as 
formulated above, was ver\' ill conditioned. When the exact answers were substituted 
as initial guess for the gains, the determinant values for each gain resulted a machine 
zero of l.OE-08. A verv small deviation of 0.01 from the exact gains was then tried and 
resulted in a determinant value of greater than l.OE + 06, as did the determinant for the 
initial guess gain values of l.O. These radically different determinant values produced 
very large gradients for ADS to work with, and in turn, ADS could not find a solution 
that closely appro.ximated the desired initial poles. 
b. E:!-er.v(iiuc .Method 

In order to overcome the ill conditioning problem, a more robust method 
was aevisea. The real matrix: (A - BC) was formed, using initial guess values for the 
gains of -1.0. The eigenvalues of this matrix are the corresponding poles of the system, 
and were then found using the Eigensystem Routines (EISPACK) subroutine RG (real 
general matrix solver) [Ref 9]. The resulting largest and smallest, complex eigenvalues 
were then compared to the largest and smallest, desired complex poles using the 
FORTRAN min and max functions for both the real and imaginary parts. Four 
equality constraints were then defined as the absolute differences betw'een the 
corresponding min or max, and real or imaginaiw' parts of the desired poles and 
calculated eigenvalues. The objective function was then set equal to the sum of the 
constraint values. ADS was then called to minimize the objective, and thus find the 
gains corresponding to the desired poles. Table 10 shows the results for this method 
using the same initial conditions as in the determinant method above. It can be seen 
that this method finds reasonably accurate values for the gains, compared to the exact 
method. 

It should be pointed out that in order to prevent the the optimizer from 
finding a local minimum corresponding to a positive real pole solution the following 
mechanisms w’ere programmed into the algorithm: 

« Strategy 133 was selected. 
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TABLE 9 

CASE 2 INITIAL GAIN APPROXIMATION USING DETERMINANTS 



Desired Pole 


Derived Pole 


-10.0+ lO.Oj 

-lO.O-lO.Oj 

-300.0 


-14.33 + 9. 66j 
-14.33-9. 66j 
-294.8 


-10.0+ lO.Oj 

-lO.O-lO.Oj 

-30.0 


+ .257 
-473.2 
-30.1 


-10.0+ lO.Oj 

-lO.O-IO.Oj 

-3000.0 


-24.15 

-582.9 

-.270 


-100.0-^ lO.Oj 
-100.0- lO.Oj 
-300.0 


-99.00 
-274.7 
+ 1.08 


-10.0+ lOO.Oj 

-lO.O-lOO.Oj 

-300.0 


-66.2 

-10.59 

-201.4 



Derived Gains 


Exact Gains 


-.125 


-.380 


+ .617 


+ .600 


-8.81 


-6.00 


+ .414 


-.920 


+ 1.51 


-.750 


+ .367 


-.600 


+ .39 1 


+ 5.02 


+ 2.03 


+ 14.1 


+ .360 


-60.0 


+ .167 


+ 1.50 


+ .863 


+ 6.01 


+ 2.94 


-303.0 


+ .362 


+ .610 


+ .343 


+ .600 


-14.13 


-303.0 



ADS Strategy 957 



• The default ADS scaling of the optimization problem was disabled (to be 
discussed in a latter section). 

• A check for positive real parts of the smallest eigenvalue was done. If a 
positive eigenvalue w'as found, then the sign of the first gain value was reversed. 

An observation that should be pointed out, as a consequence of these 
determinations, is how sensitive optimization is to problem definition. Mathematically, 
both of the methods discussed above are equally good for establishing the gains 
corresponding to the desired poles, but the second method results in always dealing 
with numerical values of essentially the same order, of magnitude, and a much better 
conditioned problem for ADS to solve. The source code for both of the approximate 
gain determination methods has been included in Appendix E. 
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TABLE 10 

INITIAL GAIN DETERMINATION USING EIGENVALUES 



Desired Pole 


Derived Pole 


-10.0+ lO.Oj 

-lO.O-lO.Oj 

-300.0 


-9.99+ 10.05j 

-9.99-10.05j 

-300.0 


-10.0+ lO.Oj 
-lO.O-lO.Oj 
-30.0 • 


-7.84 + 9.48j 

-7.84-9.48j 

-29.2 


-10.0+ lO.Oj 

-lO.O-lO.Oj 

-3000.0 


-118.1 

-10.3 

-3000.0 


-100.0+ lO.Oj 

-lOO.O-lO.Oj 

-300.0 


-65.3 + 8.76j 
-65. 3-8. 76j 
-292.0 


-10.0+ 100. Oj 

-lO.O-lOO.Oj 

-300.0 


-10.0 + 99. 7j 
- 10.0-99. 7j 
-299.9 



Derived Gains 


Exact Gains 


-.384 


-.380 


+ .599 


+ .600 


-5.98 


-6.00 


-.939 


-.920 


-.775 


-.750 


-.442 


-.600 


+ 37.6 


+ 5.02 


+ 14.6 


+ 14.1 


-364.0 


-60.0 


+ 1.11 


+ 1.50 


+ 3.25 


+ 6.0 1 


-127.1 


-303.0 


+ .603 


+ .610 


+ .598 


+ .600 


-301.0 


-303.0 



ADS Strategy 133 



6. Constraint Definitions 

The maximum overshoot and completion time constraints, as previously 
described, were applied directly to this problem. Using a manner similar to the 
maximum overshoot, a new constraint was defined for the maximum allowable value of 
the feedback control signal u(t). 

An additional feature seen in the state feedback problem is that as ADS 
begins its optimization it has no way of knowing whether or not the design variable 
values proposed for the feedback gains will cause the system to become unstable. To 
account for this, two options are available: 

1) In a reverse manner to the initial gain computations the ADS proposed gains 
could be used to compute the new corresponding system poles. These poles 
could then be evaluated for stability. 
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2) An extreme upper limit could be set for the allowable system output in order 
to monitor for stability. When this condition was reached the run could be 
terminated and an unstable solution equality constraint could be invoked. 

This second method was selected as more numerically efficient as it can be assumed 
that most proposed gains would be reasonably stable, and the cost of pole evaluation 
for every set of ADS proposed gains is high due to the need for finding eigenvalues to 
the closed loop characteristic equation. 

A stability constraint was proposed, in private discussions with Vanderplaats, 
as a way to eliminate unstable solutions. A simple way to look at this method is to 
consider the total design space. The constraints related to specific performance (such 
as completion time) define the borders of the space, and the resulting interior region is 
the domain of feasible solutions. Addition of the stability constraint essentially 
generates a "hole" in the feasible region by eliminating possible solutions at specific 
points, while allowing solutions at points close by. Normally the stability constraint 
would be set to a constant, negative value of say minus one, and thus it would never 
define an active boundary of the design space. However, when an invalid design is 
proposed by ADS, through a request for evaluation of controller gains that result in 
violating the maximum allowable output (an unstable condition), the stability 
constraint w'ould be set to a positive one, and the previous value of the objective 
function would be maintained. When ADS sees a step change in the unstable solution 
constraint, it will adjust the design variables to return this constraint negative. The 
result will be that only valid values of the design variables are considered for 
evaluation. 

Nye and Tits [Ref. 10:p. 1695] proposed the definition of hard and soft 
constraints. All the previously defined constraints for system performance (maximum 
overshoot, completion time, etc.) fall into the category of soft constraints, because they 
are set by the designer to meet flexible goals through the optimization. One type of 
hard constraint they consider is that of negative response of the input. This constraint 
simply checks to see that the response curve takes off in the same direction that the 
input goes in, otherwise the solution w'ill be divergent. This constraint was 

implemented in a similar manner to the stability constraint check above. 

7. Integration Termination 

With the introduction of integral feedback and multiple inputs, the 
performance index/objective function can become a ramp, rather than having a plateau 
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as in Case 1. To account for this effect, an additional check for a constant second 
derivative of the objective function was added to the model as a possible criterion for 
the simulation integration termination. 

8. Final Pote Determinations 

Once ADS com.Dletes the optimization, the final system gains are known. In 
order to give the designer information on how much the optimization caused the 
systems poles to be changed from the originally specified ones, the program calculates 
the final system poles. This requires a two step process: 

1) Calculation of the the final closed loop system matrix equation: A - BG 

2) Determination of the resulting matrix eigenvalues by the same procedure 
described in the approximate initial pole to gain determination using 
eigenvalues. 

C. SCALING 

In work by Chow [Ref. 11] and Gordon (Ref 12], where ADS was used in the 
analysis of poie piacement techniques, scaling was "turned off' by manipulation of the 
ADS integer work matrix (IWK). The technique to perform this operation is described 
in the ADS Manual [Ref 4:pp. 12-16]. One reason for not scaling is to increase 
computational efficiency, as the scaling calculations are not required and thus directly 
reduce the time required to run the problem. It was felt that a determination should be 
made as to how scaling effects the optimization problem when using ADSL. 

1. Scaling the General Optimization Problem 

Scaling of any numerical problem is generally considered good engineering 
practice. Vanderplaats [Ref 3:pp. 136-137] describes in detail the theory of scaling the 
optimization problem by two different methods. First, the initial problem can be 
scaled to provide relatively the same order of magnitude to the design variables, 
constraints, and objective function. In Case 2 this was .done in the initial problem 
statement. The second scaling is then done to the gradients of the individual 
parameters, so that they also have the same order of magnitude. Scaling of the 
gradients, in addition to the values, is critical because the basic optimization problem 
rests on finding a minimum, or in other words, finding where the gradient goes to zero. 

2. Scaling in ADS 

Scaling wdthin ADS is done by applying separate scaling factors to the 
objective function, and each design variable and constraint. From examination of the 
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ADS source code [Ref. 13], the design variable and constraint scaling factors were 
found to be applied as follows: 

1) Gradients are found using finite differences based upon the initial ADS call 
values. 

2) The absolute value of the gradient is taken. 

3) if the gradient is less than the ADS numerical value for zero, the gradient 
value is set to 0.1. 

4) If the gradient is greater than 1000.0, the gradient value is set to 1000.0. 

5) The scale factor is set equal to the inverse of the gradient value. 

6) The product of each design variable and its respective scale factor is taken. 

7) The product of each constraint and its respective scale factor is taken. 

8) All scale factors are maintained unchanged throughout the optimization. 

9) Upon completion of optimization, each design variable and constraint is 
aivided by its respective scale factor, to provide a final value that is not scaled. 

For the objective function the scale factor is based upon the square root of the sum of 
trie squares of the design variable scale factors. When scaling is not used, ADS sets 
the scale factors to a value of 1.0. 

3. Scaling in Case 2 

In order to investigate how scaling would affect the problem posed in Case 2, 
the initial design noted in Tables 6 through 8 was assumed to provided a reasonably 
scaled problem. Various strategy combinations were then selected and optimized with 
and without scaling applied. Table 11 shows the results of these runs, and Fig. 6.3 is a 
plot of the state responses for the optimized controller, using initial poles of 
-10.0+ lO.Oj, -10.0-10. Oj, and -300.0 and the 047 strategy for optimization. 

Trends that can be seen from these results are: 

• In all cases the strategy that had scaling applied found a lower value of the 
objective function. 

• All strategies found essentially the same optimum if there was no scaling. 

• By scaling, and using the most costly strategy (133) with respect to total 
simulations required, the lowest objective function value was found. 

In Appendix F the output listing for the 057 strategy combination with and 
without scaling can be found. Close inspection of the iteration to iteration 
optimization results reveals that the first iteration is the most critical to the problem 
solution. In the run without scaling, the optimization proceeds slowly, and relatively 
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TABLE 11 

CASE 2 SCALING RESULTS 



Strategy 


Scaling 


Objective 


Final Poles 


Total 

Simulations 


057 


Yes 


.0492 


-20.6+66.2] 

-20.6-66.2j 

-31.93 


69 


057 


No 


.0658 


-11.5 + 52.2] 

-11.5-52.2] 

-21.75 


67 


047 


Yes 


.0492 


-20.6 + 66.3] 

-20.6-66.3] 

-31.93 


69 


047 


No 


.0640 


-11.1 + 52.8] 

-11.1-52.8] 

-21.77 


79 


957 


Yes 


• .0432 


-22.5 + 75.5] 

-22.5-75.5] 

-42.57 


16 


957 


No 


.0543 


uouooo 

1 1 1 


80 


133 


Yes 


.0250 


-346. + 383.] 

-346.-383.] 

-366.2 


136 


133 


No 


.0645 


-9.35 + 54.4] 

-3.35-54.4] 

-20.9 


101 



Initial Poles 

-lO.O+lO.Oj -lO.O-lO.Oj -300.0 
Initial Objective = 0.1396 



constantly towards the minimum. The gradients are small and thus provide less 
significant data for ADS to process in the determination of search direction and step 
size. When scaling is applied, a significant initial step is made. In fact, the initial step 
"jumps over" the optimum found in the case without scaling, and then simple 
refinement is done. Thus the classical case of a relative minimum being found is seen. 
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STATE VALUES 

0.0 0.2 0.4 0.6 0.8 



C\2 




Figure 6.3 Case 2 Optimal Response 
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It can be concluded that the computational elTiciency of turning scaling off is not 
beneficial when using ADSL. 

A word of precaution about the use of ADS print control parameter 
(IPRINT) [Ref. -;p. lOj is in order. The value assigned to this parameter should be 
sciecteu :o provide a printout that reveals, at least, the optimizer termination criteria. 
Use of IPRINT= 1000, for example, prints only the initial and final optimization 
results. The maximum iteration termination could then be invoked and the user would 
not know that the result obtained from ADS was questionable. While IPRINT = 3050 
produces more output than is really useful, it is valuable for checking what is really 
happening during the optimization process, the type of scaling that occurs, and to 
reveal the complete status of the solution. In general, use of IPRINT= 1010 for cases 
when no strategy is used (ISTR.A.T= 047 or 057) or IPRINT= 1100 for cases where a 
strategy is used (ISTRy\T= 133 or 957), proved out to be the most beneficial print 
control values. 

D, INITIAL POLE SELECTION 

Different initial poles were selected to determine how initial selection effects the 
optimization process. Intuition may lead one to assume that the optimizer should find 
the same result no matter where the initial poles are placed. This estimate turns out to 
be incorrect on several counts. First, the optimizer is likely to find the minimum value 
of the objective function closest to where it initially starts. This is know as a relative 
or local minimum. Without further investigating the whole spectrum of possible values 
for the design variables, the optimizer has no way of differentiating between a relative 
minimum, and the global minimum, which represents the minimum of all relative 
minimums of the objective function [Ref 3:pp. 12-13]. In Case 2, changing the initial 
poles results in corresponding gains being sent to ADS that can vary drastically. A 
significant change in the initial design then leads to solving essentially a completely 
different optimization problem, that could have its own distinct relative minimum. 
Whether or not this minimum is the global minimum is not known. Secondly, the 
poles selected may result in an initial solution in which one or more constraints are not 
met. This is know as starting in the infeasible domain. ADS will attempt to find a 
feasible solution that will meet all the constraints, but in some cases it is unsuccessful. 
In the cases where ADS cannot meet all the constraints, it will still attempt to find the 
best solution, but with one or more violations present. This means that although ADS 
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found a solution, the designer should still try a different set of initial poles in an 
attempt to start ADS in the feasible domain. A shift in the initial pole placement can 
generate a solution that is significantly better, and one that meets all constraints. 
Table 12. shows the results for both small and large changes in the initial guess for the 
system poles using the Case 2 model. Figures 6.4a and 6.4b shows the pole placement 
results graphically (not to scale). 



TABLE 12 

CASE 2 POLE SELECTION RESULTS 



Initial 

Poles 


Initial 

Gains 


Initial 

Objective 


Final 

Poles 


Final 

Gains 


Final 

Objective 


-10.0+ lO.Oj 

-lO.O-lO.Oj 

-300.0 


-.380 
+ .600 
-6.00 


.139 


-344.2 + 342.0j 

-344.2-342.0j 

-40.88 


+ 25.3 
+ 2.64 
-962.4 


.024 


-10.0+ lO.Oj 

-lO.O-lO.Oj 

-30.0 


-.920 

-.750 

-.600 


.189 


-3.17 + 25.9j 

-3.17-25.9j 

-10.2 


-.925 

-.917 

-.699 


.143 


-10.0+ lO.Oj 

-lO.O-lO.Oj 

-3000.0 


+ 5.02 
+ 14.1 
-60.0 


.133 


-75.2+ 150.8 

-75.2-150.8j 

-46.7 


+ 2.54 

-.013 

-132.8 


.029 


-100.0+ lO.Oj 
-100.0- lO.Oj 
-300.0 


+ 1.50 
+ 6.01 
-303.0 


.026 


-386.5+ 165.3j 

-386.5-165.3j 

-60.4 


+ 20.9 
+ 3.14 
-642.4 


.023 


-lO.O+lOO.Oj 

-lO.O-lOO.Oj 

-300.0 


+ .610 
+ .600 
-303.0 


.181 


-248.8 + 2I0.7j 

-248.8-210.7j 

-60.4 


+ 12.6 
+ 1.79 
-624.4 


.024 



ADS Strategy 047 



Analysis of the results shows that: 

• The assumed initial poles of: -10.0+ 10. Oj, -lO.O-lO.Oj, and -300.0 provided a 
fairly reasonable value for the optimized objective function. 

• If the complex poles are dominant, nearly the same value for the objective 
function was found for every case of selecting different initial poles. 
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Figure 6.4a Case 2 Initial Pole Selection Results 



52 








Rnal 



Figure 6.4b Case 2 Initial Pole Selection Results (cont.) 
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• There was a wide variety of initial pole combinations that resulted in essentially 
the same value for the final objective function. 

» The average values for the resulting poles are -50, -250+ 150j, and -250-150j. 

E. CONSTRAINTS 

ADS considers two general types of constraint violations; The first, side 
constraints, are based upon the upper and lower bounds applied to each design 
variable (VUB and VLB). The second, defined constraints, are those relationships of 
the design variables to specific limits, that are developed by the designer to specify the 
system. 

1. Side Constraints 

Through side constraints, ADS prevents the optimization process from using a 
value less than the specified lower bound or greater than the desired upper bound for 
each design variable. The values of these bonds are specified based upon physical 
restrictions to tne system. For Case 2, these restrictions would be derived from the 
physical circuitry used to provide the feedback gains. If at the beginning of the 
optimization, the initial value of any design variable is set outside the allowable 
bounds, the design variable is reset to the closest bound, and a warning message is 
generated. It should be noted that this situation can easily occur using the Case 2 
ADSL model. The designer will specify a stable pole combination, and in general will 
not have a good idea of what the corresponding gain values are. In reality the gain 
values will be physically limited to finite values of usually only several orders of 
magnitude. To explore how the optimization is affected by restrictive side constraints 
multiple runs of the Case 2 model, described by Tables 5 through 8, were conducted 
with varying upper and lower bound values. Table 13 shows the results for order of 
magnitude changes in both the upper and lower bounds. 

Examination of these results shows that by making the side constraints more 
restrictive, a direct effect can be seen on the objective function value. As the 
constraints are tightened, the objective function grows in value, indicating that the 
system has a greater amount of error associated with it over time, and that a larger 
feedback control signal is required to achieve the desired response. 

As a word of caution, it should be noted that in order to suppress any design 
variables from the model, the programmer must directly reduce the number of design 
variables (NDV). Simply putting the upper and lower bounds (VUB and VLB) to the 
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TABLE 13 



CASE 2 RESULTS FOR REDUCING DESIGN VARIABLE BOUNDS 



^'ppe}\ ■ 

Bound 


Lower 

Bound 


Active 

Constraint 


Final 

Objective 


Final 

Poles 


Final 

Gains 


4- 1000.0 


-1000.0 


none 


0.0233 


-381 +247j 

-381-247j 

-47.2 


23.3 

3.05 

-978.1 


+ 100.0 


-100.0 


§3 


0.0305 


-61.9+ 135j 
-61.9-135i 
-44.9 ^ 


1.78 

-0.156 

-100.0 


+ 10.0 


-10.0 


§3 


0.0561 


-12.2 + 62.2] 

-12.2-62.2] 

-24.8 


-0.536 

-0.752 

-10.0 



Initial Poles 

-lO.O+lO.Oj -lO.O-lO.Oj -300.0 
ADS Strategy 047 



same value may seem to eliminate the variable, but it will still be taken into account. 
ADS will calculate the corresponding gradients and include them in the optimization 
with unpredictable results. 

2. Defined Constraints 

The next area to be considered was how does reducing individual constraints 
influence the optimization process. It should first be pointed out that ADS keeps 
tracks of constraints in two ways: 

1) Violated constraints: a constraint is considered as violated when its value 
exceeds 0.0 + CTMIN, where CTMIN is defaulted to a value of 0.01 for 
non-linear constraints, and 0.001 for linear constraints. 

2) Active constraints: a constraint is considered as active when its value exceeds 
0.0 + CT, where CT is defaulted to -0.03 for non-linear constraints, and 
-0.005 for linear constraints. 
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CTMIN and CT may be modified by changing their corresponding WK array values 
[Ref. 4:pp. 12-14]. 

Recalling that all equality constraints, whether linear or non-linear, must be 
equal to zero, and all inequality constraints must be equal to, or less than zero, the CT 
and CTMIN values can then be considered to generate a tolerance band for ADS to 
work within. Thus ADS never exactly meets an equality constraint of zero, but rather 
finds a value in the band between CT and CTMIN. In the the CONMIN program this 
band was known as the "constraint thickness." 

Constraint thickness comes into play for the Case 2 model with the 
specifications for maximum response (CSP), completion time (TSP), and maximum 
feedback control (USP). If these specifications are made too small, they will essentially 
become lost within the constraint thickness boundaries. In order to look at this effect 
with Case 2 an additional constraint was added after noting the findings about the X 2 
(motor current) state in the initial pole investigation. In a manner similar to that used 
for the maximum feedback control signal constraint, a maximum value for state X 2 
(X2SP) was established. This constraint was then tightened in order to observe its 
effect on the objective function and the other constraints. For this investigation, the 
CSP, TSP, and USP values were reduced to where they resulted in almost active 
constraints for the initial conditions and design described by Tables 5 through 8. The 
exact values are listed in Table 14, and Table 15 documents the results of the 
simulations conducted for reduced values of X2SP. 



TABLE 14 

CASE 2 CONSTRAINT SPECIFICATIONS FOR X2SP REDUCTION 



Constraint 

Specification 



Initial 

Model 

140% 

1.5 sec 

10 V 



Model 

5% 



Constraint Reduction 



CSP 

TSP 

USP 



0.85 sec 
1.05 V 
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TABLE 15 

CASE 2 X2SP REDUCTION RESULTS 





Final 


Final 


Final 




X2SP 


Objective 


Gains 


Poles 


Constraints 






21.55 


-772.0 


USP active 


0.75 


0.0235 


4.23 


-215.0 


lower bound 






-1000.0 


-60.2 


g^ active 






27.99 


-721.7 


USP, CSP active 


0.70 


0.0235 


4.17 


-259.9 


lower bound 






-1000.0 


-53.3 


g^ active 






25.83 


-1089.0 


X2SP active 


0.60 


0.0248 


5.53 


-189.7 








-996.0 


-48.2 








24.64 


-2146.0 


X2SP active 


0.40 


0.0350 


10.32 


-16.5 


CSP active 






-681.1 


-41.4 





Initial Poles 

-lO.O+lO.Oj -lO.O-lO.Oj -300.0 
ADS Strategy 047 



From the results in Table 15, it appears that reducing the maximum allowable 
value for the X 2 state will have varying effects on the other constraints, but in general 
will cause the objective function to increase in value. This is as would be expected, 
because with a smaller allowable maximum motor current, the controller will see more 
error over time, and require a larger feedback signal to produce the desired output 
response. An additional note of interest is that for every case investigated, the final 
gains resulted in an optimum solution with all poles being real. 

When the other constraints were tightened, individually, similar results were 
seen. It should be pointed out that once a constraint is pushed too far the model 
collapses. This collapse can be explained if one looks at the design variable space. 
With over restrictive constraints, the feasible design region disappears, and an optimum 
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design is impossible. An example of this phenomenon can be seen by looking at 
reducing the completion time constraint (TSP). Tables 16 and 17, respectively, show 
the initial conditions that were assumed and typical results. • 



TABLE 16 

CASE 2 CONSTRAINT SPECIFICATIONS FOR TSP REDUCTION 



Constraint 


Initial 


Constraint Reduction 


Specification 


Model 


Model 


CSP 


140% 


5% 


X2SP 


n.a. 


1.0 A 


USP 


10 V 


1.05 V 



TABLE 17 

CASE 2 TSP REDUCTION RESULTS 



TSP 


Final 

Objective 


Final 

Gains 


Final 

Poles 


Constraints 


1.0 


0.235 


21.5 

4.23 

-1000.0 


-772.0 

-215.0 

-60.2 


USP active 
lower bound 
g^ active 


0.06 


0.139 


-.38 

0.60 

-6.0 


-10+ lOj 

-lO-lOj 

-300.0 


TSP violated 








Initial Poles 








-10.0+ lO.Oj 


-lO.O-lO.Oj 


-300.0 



ADS Strategy 047 
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From the results of Table 17 it is apparent that if a constraint is tightened to 
the point where it becomes violated, the optimization will fail to find an improved 
solution. In most cases, .ADS will stop after three iterations and return the initial 
solution as the final solution. This fact can be used as an aid when the design process 
is bemg conducted. The designer would see no improvement for the run and could 
then immediately use this result to relax the constraint specification and rerun the 
model. Typical process time for this cycle, when used with the Case 2 example, was 
one minute (including graphic display time for the IBM 3279 terminal). 

F. INITIAL POLE SELECTION WITH ACTIVE CONSTRAINTS 

In the previous discussion of initial pole selection, none of the proposed 
constraints were active due to the large values of respective specifications that were 
initially proposed. In order to compare the essentially unconstrained results of Table 
12 to a more realistic problem, the Case 2 initial pole investigation was run again, but 
with the constraint specifications of Table 18 being applied. The results of these runs 
are shown in Table 19 and graphically in Figs. 6.5a and 6.5b (not to scale). 



TABLE 18 

CASE 2 CONSTRAINT SPECIFICATIONS FOR INITIAL POLE 

SELECTION 



Constraint 

Specification 


Initial 

Model 


Constraint Reduction 
Model 


CSP 


140% 


5% 


X2SP 


n.a. 


0.6 A 


TSP 


1.5 sec 


0.85 sec 


USP 


10 V 


1.05 V 



From comparing the results of Tables 12 and 19, one finds that the results are 
very similar with respect final objective function values. In most cases the constrained 
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TABLE 19 

CASE 2 CONSTFLAINED INITIAL POLE SELECTION RESULTS 



Initial 

Poles 


Final 

Poles 


Final 

Gains 


Final 

Objective 


Constraint 

Status 


-lO.O+lO.Oj 

-10.0-10.0] 

-300.0 


-1089 

-189.7 

-48.21 


25.8 

5.56 

-996 


.0248 


X2SP 

Active 


-10.0+ lO.Oj 

-lO.O-lO.Oj 

-30.0 


-4.08 + 29.3] 

-4.08-29.3] 

-10.6 


-.903 

-.906 

-.932 


.1229 


None 

Active 


-10.0+ 10.0] 

-10.0-10.0] 

-3000.0 


-121.7+ 130.3] 

-121.7-130.3] 

-44.6 


3.26 

0.439 

-141.8 


.0299 


X2SP 

Active 


-100.0+ 10.0] 
-100.0-10.0] 
-300.0 


-197.6 

-1070 

-47.27 


26.1 

5.57 

-.1000 


.0249 


X2SP and 
gT lower 
Axtive 


-10.0+100.0] 

-10.0-100.0] 

-300.0 


-61.7 + 39.3] 

-61.7-39.3] 

-1624 


19.57 

7.73 

-869.5 


.0264 


X2SP 

Active 
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case values were slightly higher, and this would be expected due to the optimum being 
constrained. In terms of dynamic response, however, the tightening of constraints is 
seen to have a considerable effect upon the resultant poles of the optimum designs. 

These results indicate that the sensitivity of the final design to initial pole 
selection remains, in spite of the tightened constraints. It is appparent that for the 
Case 2 system local minima abound in the design spaces that have been investigated 
thus far. Restated in another way, there are many system responses, with or without 
overshoot, that lead to a relatively small value of the objective function. Factors that 
must be futher explored in an effort to find the global minimum, or at least more 
consistent results would be: 

• investigation of more exact integration methods 
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Figure 6.5a Case 2 Constrained Initial Pole Selection Results 
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Figure 6.5b Case 2 Constrained Initial Pole Selection Results (cont.) 
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• increased precision of the optimization process through reduced tolerances of 
the ADS control parameters 

• conduct a sensitivity analysis of the various optimization parameters. 
Additionally, selection of a more practical application, together with known parameters 
associated with component hardware, would no doubt eliminate much of the design 
fiexibilty that characterizes Case 2. 

G. NOISE INPUT 

The next inquiry into the Case 2 model, as describe by Tables 5 through 8 and 
with the maximum X 2 state constraint equal to 10.0 Amps, was to extend the initial 
problem to multiple inputs. This simply required changing the Wj input from zero to a 
positive torque loading value, eliminating the negative response constraint, and 
establishing new initial and final, response and feedback control signal values. The 
resulting state response curves, for the assumed initial conditions and a constant torque 
loading, are shown in Fig. 6,6. 

In order to see how the final system gains responded to different noise inputs,, the 
comstant noise torque load (w^) of 1.0 N-m was modified to simulate a more realistic 
load condition. Two different signals were superimposed onto the constant load: 

(1) 1000 Hz sinusoid of amplitude 0.1 N-m 

(2) A random variable from a normal distribution with a standard deviation of 0.1 
N-m 

These inputs were' generated by single calls using DSL functions and can be found in 
the program listed in Appendix D. 

The results of these three runs are presented in Table 21. Figure 6.7 shows the 
final state responses for the optimized system with a constant torque loading, using 
strategy 133. Observations that can be drawn are: 

• When the constant torque noise loading is modified, with either of the 
superimposed load signals, the resulting imaginary and real poles are reduced in 
magnitude and thus the system response will be slower. 

• The final objective functions are higher in value. This would be expected as the 
system is further in error, due to the applied noise input, and additional energy 
is needed to return the system to the desired state. 

Overall, the application of a noise input demonstrated the robustness of the 
ADSL program. By the addition of one. or two simple DSL function calls, a much 
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TABLE 21 

CASE 2 NOISE INPUT RESULTS 



Noise 

Type 


Final 

Poles 


Final 

Gains 


Objective 


Total 

Simulations 


Constant 

Torque 


- 17.75 + 82.6j 

-l7.75-82.6j 

-21.29 


-.209 

-.716 

-15.2 


.399 


109 


Normal 

Dist. 


-9.602+25.41] 

-9.602-25.41] 

-180.6 


-.579 

-.001 

-13.3 


.484 


53 


Sinusoid 


-9.98 + 3.79]' 

-9.98-3.79] 

-300.8 


-.379 
+ .603 
-6.02 


1.19 


19 



Initial Poles 

-lO.O+lO.Oj -lO.O-lO.Oj -300.0 
Initial Objective = 1.920 
ADS Strategy 133 



more complicated, and realistic noise input could be investigated by the control 
designer. Any attempt at this type of analysis, using conventional methods, would be 
virtually impossible. 
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VII. MICROCOMPUTER IMPLEMENTATION OF OPTIMIZATION FOR 

CONTROLLER DESIGN 



OVERVIEW 

The ADSL program, discussed in detail in the previous chapters, was run on a 
mainframe computer (see Appendix A for a specific description). The speed and size of 
this type of system makes the initial designing of an optimal controller fairly practical. 
The same methodology can, however, be applied to a microcomputer based program. 

One system envisioned would have a microcomputer for the system controller. 
This computer would also monitor the running system, then using stochastic 
processes, generate the time averaged values of the A, B, and C matrices. These values 
could then be given to an ADSL like program that was running concurrently. New 
optimal gains would be computed, and then provided back to update the controller. 
Currently, work is in progress on the development of practical microcomputer 
..‘ontreilers, but significant problems still remain. One such problem is that of 
obtaining sufficiently high enough sampling rates, w'hile still being able to perform the 
large number of calculations required for a complex controller. 

B. PROGRAM DEVELOPMENT 

It was assumed that a microcomputer based integral controller, similar to that 
used in Case 2 could be developed, as well as the necessary A, B, and C matrix 
generation routines. A program structured to work like ADSL was written for the 
IBM Personal Computer (IBM/PC) using the Microsoft FORTRAN?? 
(MS-FORTRAN) compiler and linker programs [Ref 15]. 

1. Microcomputer Specific.s 

The microcomputer selected to implement ADSLPC was the IBM/PC models 
XT and AT, equiped with the Intel 1808? and 18028? math co-processors respectively. 
This decision was based on local availibility, general acceptance as a standard machine, 
and speed of processing. Additionally, the degree of numerical accuracy is directly 
comparable with the IBM mainframe computer [Ref 16:p. 38]. It should be pointed 
out that significant computational time differences were found, with the mainframe 
being able to run the ADSLPC program in approximately 20 seconds while 20 minutes 
were required using the microcomputer. The ADSLPC program generated 
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approximately 380 Kbytes of executable code, and this size required the PC to have 
512 Kbytes of memory, and a hard disk drive for mass storage. 

2. ADSLPC Coding 

ADS was fully implemented using the procedures outline in the y\DS 
Reference Manual [Ref 4]. DSL was replaced with an independent FORTRAN main 
program named ADSLPC. FOR. The designation of ADSLPC was assigned to this 
program as a natural follow-on to the ADSL program for a personal microcomputer 
(PC). The five functional blocks of the ADSL program, described in Chapter 4, were 
closely followed in the design of the ADSLPC program. Control of time was added to 
the main program, as well as integration routines for solving the system of differential 
equations. The following programs from the International Mathematical and 
Statistical Library (IMSL) [Ref 17] were used for determination of the din’erential 
equation solutions and eigenvalues: 

• DGEAR 

• DGRCS 

• DGRPS 

• DGRST 

• LEQTIB 

• LUDATF 

• LELMF 

• UERTST 

• EIGRF 

® EBALAF 

• EBBCKF 

• EHESSF 

• EQRH3F 

The ADS code (Version 1.12) was obtained in a format suitable for processing 
by the IBM/PC using the Microsoft compiler. It should be noted that not all 
FORTRAN compilers can be used with this version of ADS, due to the source code 
using FORTRAN Level 66 dynamic dimensioning standards. The IMSL routines 
required some slight modifications to eliminate use of unique IBM extensions to the 
FORTRAN 77 standard, and thus make them compatible with the Microsoft compiler. 
The specific changes are noted in Appendix G. 
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The main program was developed by modifying the translated DSL output file 
(FORT FORTRAN) for the Case 2 model using the eigenvalue method for the initial 
pole to gain determinations. The DSL generated code was first converted to single 
precision expressions, to be compatible with ADS, and the integration routines were 
replaced by a call to the IMSL subroutine DGEAR. This methology required that the 
derivative block of the program be broken out into its own subroutine, and that the 
Adams integration method be selected. Additionally, the EISPACK routine were 
replaced by the single precision IMSL subroutine EIGRF. 

One major difference, from the ADSL program, was to incorporate the 
evaluation of coristraints and the integration termination criteria at the end of fixed 
intervals (for the Case 2 model this was established at 0.1 secs), rather than at the end 
of every integration step. Figure 7.1 shows the interrelationships of the program calls. 
The complete source code is listed in Appendix FI, as is a listing of the batch file used 
to control processing of the program. 

C. COMPARATIVE RESULTS 

Table 22 shows how the mainframe and PC program results compare. The 
general conclusions that can be drawn from these results are: 

• ADSLPC can be used to find optimal controller gains for an integral feed back 
controller, with results comparable, to tho.se found using ADSL. 

• The difference in final objective function values can be attributed to finding 
local minima, which can be mo.st likely traced backed to using a different, and 
less accurate integration method. 

• As with the ADSL model, further improvement of the numerical tolerances 
could result in additional reduction of the objective function. 

• The obtained results compare favorably with the unsealed ADSL results found 
in Table 1 1. 

• ADSLPC could be a powerful tool for optimal, on-line updating of 
microcomputer controller processes. 
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INITIALIZE ADS AND MODEL PARAMETERS 

nNITLZ) 



\r 

^CALL ADS FOR DESIGN VARIABLES 
(Pole to Gain Optimization) 

(INTT) 



CALL EIGRF FOR EIGENVALUES 
(Find poles and compare to desired) 



Y 

CAIX ADS FOR DESIGN VARIABLES ^ 
(Controller Optimization) 



^ ESTABLISH TIME STEP THEN 
CALL DGEAR FOR INTEGRATION 



EVALUATE INTEGRALS FOR EACH 
STATE AND THE OBJECTIVE 
(DERIV) 



V 

EVALUATE CONSTRAINTS AND 
- TERMINATION CRITERIA 
(DYNAMIC) 



If 

PRINTOUT FINAL INTEGRATION HISTORY 
THEN FIND FINAL POLES AND GAINS ' 
(TERMINAL) 



ENDRUN 



Figure 7.1 ADSL Program Segments 
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TABLE 22 

COMPAI^\TIVE RESULTS FOR ADSL AND ADSLPC PROGRyWlS 





ADSL 


ADSLPC 


Initial Poles 
(desired) 


-10.0+ lO.Oj 

-lO.O-lO.Oj 

-300.0 


-10.0+ 10.0] 
-10.0-10.0] 
-300.0 


Initial Poles 
(calculated) 


-9.99+ 10.05] 
-9.99- 10.05] 
-300.0 


-9.16+ 10.00] 

-9.16-10.00] 

-299.5 


Initial Gains 
(calculated) 


-0.384 
+ 0..599 
-5.98 


-0.431 
+ 0.589 
-5.71 


ADS Strategy/ 
pole to gains) 


133 


133 


Final Poles 


-79.7+ 154.7] 

-79.7-154.7] 

-47.7 


-15.9 + 62.3] 

-15.9-62.3] 

-28.4 


Final' Gains 


+ 0.036 
+ 2.79 
-144.2 


-0.495 

-0.698 

-11.7 


Final Objective 


.0285 


.0526 


ADS Strategy 


047 


047 


Evaluations 


135 


49 
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Vlll. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The following conclusions can be draw from the previous disccussions; 

1. A complicated, constrained, multiple input, multiple output (MIMO) controller 
model can be developed using time domain methods and standard Dynamic 
Simulation Language (DSL) techniques. This model can then be optimized 
based upon a performance index of choice using the Automated Design 
Synthesis (ADS) program. 

2. The ADSL methodology provides a robust way in which stdte feedback gains 
can be obtained without use of matrix approximations associated with the 
classical Riccati equation solution. 

3. For MIMO systems in which sinusoidal and random noise inputs were 
simulated, optimal feedback gains could be found. In general they were such as 
to slow the system down, but closely approximated the gains found for a 
constant noise input. 

4. The ADS optimization code, with an eigenvalue solver routine, can be 
employed to solve the problem of finding the system poles corresponding to 
given gains for any complex MIMO controller. 

5. From the studies described herein, the following approaches were determined to 
be well suited for developing ADSL programs: 

• Use a stiff integration method 

• For initial pole to gain determination, if an exact determination is 
impractical, use ADS and an eigenvalue solver method as developed in this 
study. 

• For basic controller optimization, application of the Method of Feasible 
Directions strategy (047) proved to be the most practical overall, but the 
Sequential Unconstrained Minimization Techinque using an external 
penalty function strategy (133), in general, provided better optimums at the 
expense of more iterations. 

• Order of magnitude changes of the initial poles proved to be a reasonable 
method for establishing the validity of the optimization results for the final 
gains. 
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• If the final solution has active constraints, these should be investigated 
individually, by reduction of the corresponding constraint specification. 

6. For the state feedback controller considered in this study, many local optimums 
can be found. The reasons for this appears to line in the definition of the 
performance index and the relative freedom from severe design constraints. 

7. The basic ADSL program methodology can, be down scaled to run on a 
personal computer, thus opening the possibility to using optimization in on-line 
applications. 

B. RECOMMENDATIONS 

The following recommendations are offered: 

1. A simplified state feedback model should be investigated as the results obtained 
for the integral state feedback model demonstrated that the selected design 
space had many local minima. An acutal controller, with known physical 
parameters, would likely have a w'ell behaved design space. This would allow 
for gaining a clearer understanding of the interactions investigated in this study. 
Additionally, the results of actual tests could then be applied to verify the 
ADSL predictions. 

2. A large scale controller model should be optimized. One such controller model 
has already been developed by NASA and has been documented by Merrill 
[Ref 17]. 

3. The methodology should be expanded to provide for online interaction, and 
recursive optimization capability. 

4. The same methodology used wdth ADSL can be applied to a micro-computer 
based design system using trimmed down codes. One possible system would be 
replacement of DSL by the IBM PC-Engineering Simulation Program and ADS 
with the Engineering Design Optimization program Microdot, which is a 
reduced version of the ADS code in a format for use on a micro-computer. 
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APPENDIX A 

ADS-DSL PROGRAM INTERFACE SPECIFICS 



1. COMPUTER DESCRIPTION 

The ADSL program was initially developed at the Naval Postgraduate School W. 
R. Church Computer Center using a IBM 3033 Attached Processor System (16 
Megabytes) loosely coupled with a IBM 3033 Model U (16 Megabytes) and a IBM 
4381 Model Ml (8 Megabytes). Interactive computing was provided under VM/SP 
CMS, and batch-processing under MVS with JES3 networking. 

2. SPECIFICATIONS USED FOR COMPILING ADS (VERSION 1.10) 

a. FORTIU\N V;S 

b. Auto-Double Precision (DBL) 

c. Optimization Level 3 

d. Language Level 66 

3. DSL/VS (VERSION 1.1) IN THE VM MODE 

The resulting TEXT file from the compilation of ADS was converted to a 
FORTRAN library file using the online library generation procedure [Ref 19:pp. 
35-36]. The library was named MYLIB TXTLIB in order to take advantage of the 
default DSL user library specification for the VM environment. 

The default DSL executive was modified to streamline the interfacing of of ADS 
and DSL output, and is shown in Appendix B as the ADSL EXEC. 

4. DSL/VS (VERSION 1.1) IN THE MVS MODE 

The ADS TEXT file was also placed onto the Mass Storage System (MSS) as a 
user library in order that DLS/VS could access ADS using the MVS environment. The 
standard job control statements for using DSL must then be modified to concatenate 
the ADS, LINPACK, and EISPACK libraries as shown in Appendix B. The job 
control for compiling ADS on MSS is shown below. 

5. .JOB CONTROL FOR FILING ADS ON THE MASS STORAGE SYSTEM 

//FILEADS JOB ( 0180 , 9999 ),' FILEADS ’ ,CLASS=C 

// EXEC FORTVCL,PARM. FORT= ' LVL( 66 ) ,AD( DBL ) , NOMAP ,OPT( 5 ) ,NOS,NOX* 

//FORT.SYSIN OD 
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REPLACE THIS LINE WITH THE ADS SOURCE CODE 



/* 

//LKED.SYSLMOD OD UNIT=3330V ,MSVGP=PUB4C , 

// DSN=KSS.F0I80.ADS( ADS), 

// DIS? = ( NEW, CATLG, DELETE ) ,SPACE = ( CYL, ( 4,4 ,5 ) ) 
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APPENDIX B 

STEPS TO RUN ADSL AT NAVAL POSTGRADUATE SCHOOL 

1. VM ENVIRONMENT 

a. Logon in the VM environment (note that the user must have have previously 
requested an increase of his default virtual machine size from 750 Kilobytes to at least 
1500 Kilobytes through the Computer Center Accounting Office). 

b. Define increased storage by using: 

DEF STORAGE 1500K 
I CMS 

c. Linkto the optimization disk with the command: 

LINKTO 0180P 191 OPTIM 

d. Link to the DSL and MATHPACK Disks with the command: 

DSLINK 

This command executes the DSLINK executive procedure described below. 

e. Run the ADSL program using: 

ADSL FN <PLOT | TEK> 

FN is the user provided filename containing the DSL input program. The 
filetype of the input program must be DSL. PLOT or TEK is optional, but required if 
any plotting is called. TEK is used when TEK-618 or IBM 3279 terminal plotting is 
desired, while PLOT is used for all other plotting devices. Also note that the user can 
not have a file name MAIN TEXT on his user disk, as this file will cause DSL to load 
abnormally. 

f The specified tabled output will be saved in a file named FN LISTING, 
g. Plotting will then be done on the specified graphic output device and printer 
plots will be concatenated to the FN LISTING file. 

2. DSLINK EXEC LISTING 

EXEC LINKTO DSL 
EXEC LINKTO GDDM<+0 191 GDOMR 
EXEC LINKTO PLI5 191 PLI5 
EXEC LINKTO 0062P 191 EISPACK 
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3. ADSL EXEC LISTING 



x/x Install REXX If not here. 
aCONTROL OFF NOMSG 
EXEC REX I 

8IF SRETCODE = 0 EXEC 0 81 82 85 85 86 87 88 89 810 311 812 813 8U 815 816 817 818 319 820 

8EXIT 3RETC0DE 



(end of REX Installation! •♦/ 



/K REXX EXEC TO PERFORM COMPLETE DSL SIMULATION «/ 

«/ 

/« dSL/VS system ... lEM Corporation 198A 

/>« ♦♦♦♦ LATEST UPDATE: 4/03/84 *♦/ 

- 

Trace value 'OFF* 



/«ARG FN FT FM DISP **(" GRAF REST«/ 

ARO FN GRAF 
IF FN 5 THEN 
TELL: Do 

Say 'The ADSL command makes a complete DSL run In the foresround. by' 

Say .' executing the following EXECs in order; DSLT • DSLC. and DSLG.' 

Say 'The plot post-processor GRaFAEL may also be Invoked at the same' 

Say 'time as an option to produce graphic output at a selected graphic' 

Say 'output device.' 

Say * ' 

SAY 'COMMAND FORMAT : ADSL <FN> <PL0T> ' 

Say * * 

Say *"FN" Is the filename of a previously created DSL program file.' 

Say ' * 

Say 'NOTE that the flletype must be DSL.' 

SAY ' ' 

Say '"PLOT" Is an optional parameter which specifies If graphic post-* 

Say 'processing Is desired. Leaving It out. means no plotting Is done..' 

Say 'Use the command "PLOT" for printer plots or "TEK" for any other* 

SAY ' displays . ' 

End 

IF FN = "" THEN EXIT 
/M Invoke DSL Translator «/ 

FT='DSL’ 

FM=' A' 

DISP='FILE* 

/X GRAF= ' GRAFAEL* •♦/ 

•EXEC DSLT' fn ft fm 
If rc >s 8 then Exit rc 
COMPILE: 

/„ ♦ Invoke Fortran Compiler »*/ 

Do 

SAY * »»»> RUNNING VS FORTRAN COMPILER <«<<<' 

End 

'FORTVS FORT (LANGLVLC77) OPTC3) NOXREF SOMAP NOSDUMP NOTF NOTERM' 

/«::: "FORTVS" INVOKES THE VS FORTRAN COMPILER (PROG, a 5748-F03 ) :::«/ 

/«::: Txtllbs assumed; VALTLIB VFQRTLIB ;::«/ 

IF RC >= 8 THEN 

Do Z** Type Error Diagnostics *•/ 

Say ' * 

Say * ERROR (S) DETECTED IN COMPILATION PHASE anana’ 

Say ' ' 

Say ' Fortran listing will be displayed by XEDIT; Enter "OUIT" to exit* 

Say 

Say • ' 

Queue 'PREFIX OFF* 

•XEDIT FORT LISTING' 

EXIT 33333 
End 
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IF RC >=B THEN EXIT 

/* *♦♦♦ Test for existence of user's MYLIS TXTLIB file »/ 

’STATE MYLI3 TXTLIB* /“tt; MYLI3 may be renamed : s : •»/ 

If rc a 0 verify the VS Fortran Library 

names assumed: VALTLIB S VFORTLIB : : : »/ 

/ - ALSO '.INKTO EISFACK AND LINKPaCX LIBRARIES «/ 

'HEN SLCSAL 'XTLI3 CSLSIH SSLBLKS VAlTLIB VFORTLIB CMSLIB. ' 

L’lSPA'SK 

else ’GLOBAL TXTLI3 MYLIB* /•*::: HYLIB may oe renamed :::«/ * 

’DSLS I M DSLSLK3 VALTLIB VFORTLIB CMSLIB 
'EISPACK LINPACK’ 

/« erase old "points” file If It exists..**/ 

•STATE PLCT DATA’ 

If rc ^ 0 then 'ERASE PLOT DATA' 

'FILEDEF 06 DISK OUTP DATA (RECFM F9A LRECL 133 3LKSIDE 665’ 

•FILcDEF 05 DISK AUX DATA' (RECFM F3 LRECL BO BLKSIZE BOO’ 

•FILEDEF 12 DISK PLTC DATA (RECFM FB LRECL BO 3LKSIZE BOO’ 

•FILEDEF 13 DISK TABL DATA (RECFM FB LRECL BO BLKSI2E BOO* 

’FILEDEF lA DISK PLOT DATA (RECFM VB 3LK3IDE 102A’ 

'FILEDEF 15 DISK DATA DATA (RECFM F3 LRECL 30 BLKSI2E BOO’ 

■FTLEOEF is disk 5RTSAV DATA (RECFM VB3 LRECL 5A0 SLKSI2E 5A4 ’ 

•FILEDEF SI DISK 5 UN DATA i RECFM =3 LRECL 30 BLKSIDE 300’ 

/«• adcltioral user 'lie definitions you may be added here :::*»/ 

Do 

SAY • • 

DAY • >>>>>>> punning DSL/VS SIMULATOR <<<<<<<’ 

End 

^UNi 

OP let EHSG CFF’ 

•GLOBAL ^OAOLIB VFL0DLI3’ 

•LOAD FORT (CLEAR RESET MAIN START NO MAP * 

slmrc = rc 

•CP SET EMSG TEXT’ 

If slmrc >= B then 
Do 

Say ’ ' 

Say 

Say ' ttnnnn ERROR(S) DETECTED IN SIMULATION PHASE «««««• 

Say 

Say ' ' 

End 
CONT ; 

•CP SET I MSG ON’ 

If slmrc >= B then Exit 33333 
IF GRAF = ’PLOT’ THEN ’GRAFAEL’ 

IF GRAF 3 ’TEK’ THEN ’GRAFAEL’ 

•ERASE ’FN ’ LISTING A’ 

IF DISP = ’PLOT’ THEN ’COPY OUTP DATA A PPLT DATA A ’ FN * LISTING A’ 

ELSE 'RENAME OUTP DATA A ’ FN ’ LISTING A' 

•X * FN • LISTING ' 

Exit /•* end of program "/ 

4. SAMPLE JOB CONTROL FOR MVS ENVIRONMENT 



//ADSL JOB (0180,9999), 'ADSL' ,CLASS=B 

//TEST EXEC DSLVS,PARM.C='N0MAP,N0X,N0L,N0S,0PT(3)' 

//^FORMAT PR,DDMAME=JESMSG,COPIES=0 

//i^FORMAT PR , DDMAME =SYSMSG , COPI ES=0 

//i^FORMAT PR,DDNAME=JESJCL,COPIES=0 

//TEST EXEC DSLVS,PARM.C='N0MAP,N0X,N0L,N0S,0PT(3)' 

//DSL. INPUT DD * 
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INSERT ADSL PROGRAM CODE HERE 



END 






STOP 






//L.SYSLIB 


DD 




// 


OD 




// 


DD 




// 


DD 




// 


DD 




// 


DD 




// 


DD 


DISP=SHR,DSN=MSS.F0180.ADS( ADS) 


// 


DD 


DISP=SHR ,DSN=MSS . SYSS . EISPACK . LOAD 


// 


DD 


DISP =SHR , DSN=MSS . SYS2 . LINPACK . LOAD 








// 
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APPENDIX C 

EXAMPLE SISO PROGRAM LISTING 



TITLE SISO 2ND ORDER SYSTEM-OGATA SECTION 7-S WITH CONSTRAINTS 
LISTING OF ALL VARIABLES 



* A 


HORK MATRIX FOR ADS 


* DF 


VECTOR OF DESIGN VARIABLE GREADIENTS 


* G 


VECTOR OF CONSTRAINT VALUES 


* H2 


VALUE OF THE RESPONSE OSCILATION 


^ I 


LOOP COUNTER 


* IC 


VECTOR OF CONSTRAINT GRADIENTS 


* IDG 


VECTOR OF CONSTRAINT TYPES 


* IGRAD 


SCALAR FOR INDICATING METHOD FOR OBTAINING GRADIENTS 


^ INFO 


SCALAR INDICATING OPTIMIZATION STATUS 


* lONED 


SCALAR INDICATING TYPE OF ONE-DIMESIONAL SEACH 


* lOPT 


SCALAR INDICATING TYPE OF OPTIMIZATION 


* IPRINT 


SCALAR INDICATING TYPE OF PRINT CONTROL FOR ADS 


* ISTRAT 


SCALAR INDICATING TYPE OF STRATEGY 


* IWK 


INTEGER WORK VECTOR FOR ADS 


* LSTOBJ 


VALUE OF THE OBJECTIVE FUNCTION AT PREVIOUS INTEGRATION STEP 



LSTTIM VALUE OF TIME AT THE PREVIOUS INTEGRATION STEP 



* NCOLA 


SCALAR COLUMN DIMENSION OF A MATRIX 


* NCON 


NUMBER OF CONSTRAINTS 


* NDV 


NUMBER OF DESIGN VARIABLES 


* NGT 


SCALAR ADS CONTROL VARIABLE 


* NRA 


SCALAR ROW DIMENSION OF A MATRIX 


* NRIWK 


SCALAR DIMENSION OF IWK VECTOR 


* NRWK 


SCALAR DIMENSION OF WK VECTOR 


* OBJ 


OBJECTIVE FUNCTION VALUE 


* R 


VALUE OF THE INPUT TO THE SYSTEM 


* TAU 


SYSTEM TIME CONSTANT 
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* 


TIME 


DSL TIMER VALUE 




* 


VLB 


VECTOR OF DESIGN 


VARIABLE LOWER BOUNDS 


* 


VUB 


VECTOR OF DESIGN 


VARIABLE UPPER BOUNDS 




WK 


ADS WORK VECTOR 




* 


X 


VECTOR OF DESIGN 


VARIABLES 




Y 


VALUE OF THE SYSTEM OUTPUT 




YMAX 


MAXIMUM VALUE OF 


THE SYSTEM OUTPUT 




YOSC 


VECTOR OF SYSTEM 


OUTPUT VALUES AT PREVIOUS INTEGRATION STEPS 



* PREAMBLE TO INITIALIZE THE ADS PARAMETERS 
FIXED ISTRAT,IOPT ,IONED ,IPRINT ,INFO ,IGRAD 

FIXED NDV ,NCON ,IDG ,NGT ,IC ,NRA ,NCOLA ,NRN:L ,IWK ,NRIWK ,I 

* DIMENSIONING REQUIREMENTS ARE FOUND IN THE ADS MANUAL TABLE 5 
ARRAYS MUST BE ASSIGNED USING THE DIMENSION STATEMENT IN DSL 

D DIMENSION A (2,6) 

VECTORS CAN BE ASSIGNED USING THE ARRAY STATEMENT IN DSL 
ARRAY WK(1000),IWK(S00),DF(2) 

ARRAY X( 2 ) ,VLB( 2 ) ,VUB( 2 ) ,G( S ) ,IDG( 5 ) ,IC( 5 ) 

INITIALIZATION OF ADS INPUT CONSTANTS (ADS MANUAL TABLE 5) 
PARAM NRA=1 ,NCOLA=l ,NRWK=1000 ,NRIWK=S00 
PARAM IGRAD=0 ,INFO=0 ,NDV=1 ,NCON=S 

SETTING ALL CONSTRAINTS TO NON-LINEAR, INEQUALITY 
TABLE IDG(1-S)=0,0,0,0,0 

INITIAL GUESS FOR DESIGN VARIABLE ZETA 
TABLE X(l)=l. 

* UPPER AND LONER BOUNDS FOR DESIGN VARIABLE 
TABLE VLB(1)=,01 , VUB(1)=2.0 

SET ADS OPTIMIZATION STRATEGY (ADS MANUAL PART 4.6) 

PARAM ISTRAT=1, IOPT=3, I0NED=3, IPRINT=2020 

*SET ADS PRItrr PARAMETER (ADS MANUAL TABLE 5) 
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PARAM IPRINT=1000 



^ SET DSL INTEGRATION OUTPUT PARAMETERS (DSL MANUAL PART 4.2) 
METHOD STIFF 
:CNTROL 'INTIM=45 
SAVE Y>OBJ 

GRAPH (DE=TEK618) TIME, Y, OBJ 
PRINT Y,03J 

INITIAL 

^ CALL TO ADS TO GET THE DESIGN VARIABLE FOR THE SYSTEM 
^ CALLING STATEMENT SPECIFIED IN ADS MANUAL PART 4.0 

CALL AOS( INFO,ISTRAT,IOPT,IONED,IPRINT,IGRAD,. . . 
NDV,NCCN,X,VLB,VUB,CBJ,G,IDG,NGT,IC, . . . 

DF ,A ,NRA ,NCOLA ,WK ,NRWK ,IWK,NRIWK ) 

^ THIS STATEMENT TURNS ON THE PRINTER OUTPUT ONLY ON THE LAST RUN 

* WHEN INFO HAS BEEN RESET TO 0 BY ADS TO INDICATE OPTIMIZATION 

* HAS BEEN COMPLETED 

IF (INFO.EQ.O) DELPRT=1.5 

*THIS STATEMENT TURNS ON THE RECORDING OF PLOTTER VALUES 
IF (INFO.E(5.0) DELPLT=.4 

* RE-INITIALI2E CONSTRAINT CHECK VALUES AFTER CALL TO ADS 

YMAX=0 . 

TAU=0. 

LST03J=0. 

LSTTIM=0. 

Y0SC1=0. 

YOSC2=0. 

YOSCS=0. 

DO 15 I=1,NC0N 
15 G(I)=1.D20 
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DERIVATIVE 



NOSORT 

^ INPUT TO THE MODEL (DSL MANUAL SECTION 5.1) 
^=STEP< 0. ) 



OUTPUT OF SECOND ORDER LAG SYSTEM (DSL MANUAL SECTION 5.4) 
Y=CMPXPL(0. ,0. ,X(1),1. ,R) 

ITAE PERFORMANCE INDEX TO MINIMIZE AMPLITUDE (OGATA P. 299) 
OBJ=INTGRL( 0 . ,TIME5^ABS( R-Y ) ) 



DYNAMIC 

UPDATE CONSTRAINT CHECK VALUES AT EACH INTEGRATION STEP 
Y0SC1=Y0SC2 
. Y0SC2=Y0SC3 
Y0SC3=Y 

H2=A8S(Y0SC2-R) 

IF(TAU.NE.O. )GOTO 20 
IF(Y.GE. (0.63212»«R) )TAU=TIME 
20 IF(Y.GT.YMAX) YMAX=Y 

MAX AMPLITUDE CONSTRAINT ( <140X) 

IF(Y.EQ.YMAX) G( 1 )=YMAX-( 1 ,4»fR ) 

MAX DELAY TIME CONSTRAINT ( <1.5 sec) 

IF((Y.LE.( .S»«R)).ANO.(Y.E(?.YMAX)) G( 2 )=TIME-1.5 

MAX RISE TIME CONSTRAINT ( <3.5 sec) 

IF( (Y.LE.R ). AND. (Y.E(3.YMAX) ) G( 3 )=TIME-3 .5 

MAX PEAK TIME CONSTRAINT ( <4.5 sec) 

IF(Y.EQ.YMAX) G( 4 )=TIME-4 . 5 

MAX COMPLETION TIME CONSTRAINTS ( <12 sec) 

IF(G(5).LT.O. )G0T0 30 
FOR UNDER DAMPED CASES 
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IF( (H2.GT.ABS(Y0SC1-R) ).AND.(H2.GT.ABS . . . 
(Y0SC3-R)).AN0.(H2.LE.( .01*R))) G( 5 )=TIME-12 . 

* FOR OVER DAMPED AND CRITICALLY DAMPED CASES 

IF( ( (Y-R).LT.( .01^R)).AND.(TIME.GE.(5.*TAU)). . . 

. AND . ( T AU . NE . 0 . ) . AND . ( Y . E Q . YMAX ) ) G ( 5 ) =TIME - 1 2 . 

* INTEGRATION COMPLETE CHECK USING THE OBJECTIVE FUNCTION SLOPE 

* CHECK FOR FIST INTEGRATION STEP 

30 IFKOBJ.LE.O. ).OR.(LSTOBJ.LE.O. ))GOTO 40 

* DETERMINE SLOPE AND CHECK FOR PAST COMPLETION TIME - IF MET 

* ENDRUN DIRECTS PROGRAM TO THE TERMINAL SEGMENT OTHERWISE THE 

* NEXT INTEGRATION STEP IS TAKEN UNTIL FINTIM IS MET 

IF( nOBJ-LSTOBJ)/( TIME -LSTTIM).LT.l.E-5). AND. . . . 

( ( (H2.LE.1 .01*R) ).AND. ( H2.GT.ABS( YOSCl-R)) . . . 
.AND.(H2.GT.ABS(YOSC3R))).OR.( (TIME. 
GE.(5.*TAU)).AND.(TAU.NE.O. ))))CALL ENDRUN 
40 LSTOBJ=OBJ 

LSTTIM=TIME 

TERMINAL 

* CHECK TO OPTIMIZATION COMPLETE TO EXIT DSL OTHERWISE RERUN WILL 
START THE PROGRAM OVER AT THE INITIAL SEGMENT 

IF (INFO.EQ.O) CALL ENDJOB 
CALL RERUN 



END 

STOP 



84 



APPENDIX D 

STATE VARIABLE PROGRAM LISTINGS 



TITLE DC MOTOR INTEGRAL CONTROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 
^ LISTING OF ALL VAPIABLES 



* AA 


MATRIX OF STATE COEFFICIENTS 


* BB 


VECTOR OF CONTROL COEFFICIENTS 


* C 


SCALAR OUTPUT VALUE 


* CC 


MATRIX OF STATE OUTPUT COEFFICEINTS 


* CF 


FINAL EXPECTED OUTPUT VALUE 


* Cl 


INITIAL OUTPUT VALUE 


* CON 


VECTOR OF CONSTRAINT VALUES 


* CSP 


MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 


* DELC 


DIFFERENCE OF THE ACTUAL OUTPUT TO THE INITIAL OUTPUT 


* DELU 


DIFFERENCE OF THE ACTUAL CONTROL TO THE INITIAL CONTROL 


* 0ELX2 


DIFFERENCE OF THE ACTUAL X2 STATE TO THE INITIAL X2 STATE 


* DV- 


VECTOR OF DESIGN VARIABLE GREADIENTS 


^ F 


VECTOR OF INPUT COEFFICIENTS 


* FVl 


EISPAC HORK VECTOR 


G 


VECTOR OF DESIGN VARIABLE GAIN VALUES 


* H 


VECTOR OF OUTPUT COEFFICEINTS FOR THE INPUTS 


^ I 


LOOP COUNTER 


^ IC 


VECTOR OF CONSTRAINT GRADIENTS 


* IDG 


VECTOR OF CONSTRAINT TYPES 


* lERl 


LINPAC ERROR INDICATOR 


* lERR 


EISPAC ERROR INDICATOR 


* IGRAD 


SCALAR FOR INDICATING METHOD FOR OBTAINING GRADIENTS 


^ II 


COUNTER OF NUMBER OF FUNCTION EVALUATIONS 


^ INFO 


SCALAR INDICATING OPTIMIZATION STATUS 


* lONED 


SCALAR INDICATING TYPE OF ONE-DIMESIONAL SEACH 


* lOPT 


SCALAR INDICATING TYPE OF OPTIMIZATION 


* IPRINT 


SCALAR INDICATING TYPE OF PRINT CONTROL FOR ADS 
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* IFVT LINPAC INTEGER WORK VECTOR 



* ISTRAT SCALAR INDICATING TYPE OF STRATEGY 

^ IVl =I3PAC I^iTHGcR WORK VECTOR 

^ :w:k :ntcSer work /hctcr for ads 

* J LOOP COUNTER 

LSLOPE SLOPE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 

* LSTOBJ VALUE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 
^ LSTTIM VALUE OF TINE AT THE PREVIOUS INTEGRATION STEP 

* MATZ EISPAC SCALAR VALUE TO INDICATOR EIGENVALUE DETERMINATION 

* MEAN DSL NORMAL DISTRIBUTION MEAN 

NCOLA SCALAR COLUMN DIMENSION OF WA MATRIX 

* NCON NUMSER OF CONSTRAINTS 

MOV NUMBER OF DESIGN VARIABLES 

^ NGT SCALAR ADS CONTROL VARIABLE 

^ NM EISPAC OIMCNSION OF MATRIX SIZE FOR EIGENVALUE DETERMIANTION 

NN EISPAC DIMENSION OF MATRIX SIZE FOR EIGENVALUE DETERMIANTION 

* NRA SCALAR ROW DIMENSION OF WA MATRIX 

* NRIWK SCALAR DIMENSION OF IWK VECTOR 

* NRWK SCALAR DIMENSION OF WK VECTOR 

* OBJ OBJECTIVE FUNCTION VALUE 

* OOBJ OBJECTIVE FUNCTION FOR INTEGRATION 

* PI DESIRED INITIAL COMPLEX POLE VALUE 

* P2 DESIRED INITIAL COMPLEX POLE VALUE 

* PI DESIRED INITIAL COMPLEX POLE VALUE 

* PA MATRIX OF INITIAL GAIN DETERMINATION COEFFICIENTS 

* Q STATE WEIGHTING MATRIX 

* R CONTROL WEIGHTING VECTOR 

* SD STANDARD DEVIATION OF DSL NORMAL FUNCTION 

* SEED RANDOM MUMBER SEED FOR DSL NORMAL FUNCTION 

* TAU SYSTEM TIME CONSTANT 

* TC TIME WHEN COMPLETION TIME SPECIFICATION IS MET 

* TIME DSL TIMER VALUE 

* TSP MAXIMUM COMPLETION TIME 

* U SCALAR FEEDBACK CONTROL VALUE 

* UI INITIAL VALUE OF THE CONTROL VECTOR 
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* UMAX 


MAXIMUM CONTROL VALUE 


5^ USP 


MAXIMUM VALUE OF THE CONTROL VECTOR 


^ UT 


TRANSLATED FEEDBACK CONTROL VALUE 


7 


MATRIX FOR FINAL POLE DETERMINATION USING EISPACK 


* 


'ECTOR OF DESIGN VARIABLE LOWER BOUNDS 


* VUB 


VECTOR OF DESIGN VARIABLE UPPER BOUNDS 


^ W1 


DISTURBANCE VALUE 


CM 


REFERENCE SETTING 


^ WA- 


ADS WORK VECTOR 


* WI 


VECTOR OF FINAL IMAGINARY POLES 


^ WL 


CIRCULAR FREQUENCY OF SINUSOIDAL DISTRUBANCE 


> WR 


VECTOR OF FINAL REAL POLES 


* XI 


MOTOR SHAFT SPEED STATE VALUE 


XIT 


TRANSLATION OF XI STATE 


^ :<2 


MOTOR CURRENT STATE VALUE 


^ X2I 


INITIAL VALUE OF THE X2 STATE 


^ X2MAX 


MAXIMUM VALUE OF THE X2 STATE 


* X2SP 


MAXIMUM SPECIFIED VALUE OF THE X2 STATE 


* XXI 


FUNCTION DESCRIBING THE XI STATE FOR INTEGRATION 


* XX2 


FUNCTION DESCRIBING THE X2 STATE FOR INTEGRATION 


* YOSCl 


VALUE OF SYSTEM OUTPUT OSCI CATION AT TWO PREVIOUS INTEGRATION STEPS 


* YOSC2 


VALUE OF SYSTEM OUTPUT OSCILATION AT PREVIOUS INTEGRATION STEP 


YOSC3 


VALUE OF SYSTEM OUTPUT OSCILATION 


* Z 


EISPACK DUMMY SCALAR 



INITIALIZE THE ADS PARAMETERS 
FIXED NRA ,NGT ,NCOLA ,NRWK ,NRIWK ,IGRAD ,NDV ,NCON,IDG ,IC ,IWK 
FIXED ISTRAT,IOPT,IONED,IPRINT,INFO,IPVT 

ARRAY DV(4),VLB(4),VUB(A),C0N(6),IDG16 ),IC(6),DF(4),WK(10DD),IWK(1D0D) 
D DIMENSION WA(4D,4D ) 

PARAM NRA=4D ,NCOLA=40 ,NRKK=1000 ,NRIWK=1DD0 ,NDV=3 ,NCON=6 ,IGRAD=D ,INFO=D 
ESTABLISH TYPES OF CONSTRAINTS (ADS MANUAL TABLE 5) 

TABLE IDG( 1-6)=35«D,2*2,D 
* LOWER BOUND FOR DESIGN VARIABLES 
TABLE VLB(l-3)=3»-lDDD. 
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* UPPER BOUND FOR DESIGN VARIABLES 



TABLE VUB( 1-3 )=3^100C. 

^ ADS OPTIMIZATION METHOD SELECTION 

PARAM ISTRAT=0, IOPT=^, IONED=7, I?RINT=1010 



INITIALIZE THE DSL PARAMETERS 
PRINT CBJ,C,X1,X2,U 
SAVE (RUT4) OBJ,Xl,X2,U,C 
CONTROL DELPLT=O.O05,DELPRT=O.Dl 
C0^4TR0L FINTIM=4. 

METHOD STIFF 
RELERR OBJ=l.E-^ 

ABSERR GBJ=l.E-4 



* INITIALIZE THE MODEL PARAMETERS 
FIXED I,II,IER1,IERR,NN,N.M,IV1,MATZ,SEED 
COMPLEX P1,P2,P3 

ARRAY BB(3),H( 2),R( 2),IPVT(3),G(3),WR(3),WI(3),IV1(3),FV1(3) 

D DIMENSION AA( 3 ,3 J ,CC( 2 , 2 ) ,Q( 2, 2 ) ,PA( 3 ,3 ) , F( 2 , 2 ) ,V( 3 ,3 ) ,Z( 3 ,3 ) 
ZERO ALL MATRICES 
D DATA AA,F,CC,Q/21^0V 
TABLE 3B( l-3)=3^0. ,R( l-2)=3^0, ,H( l-2)=3^0. 

^ SET CONSTRAINT SPECIFICATIONS 

CI=INITIAL OUTPUT VALUE; CF=FINAL EXPECTED OUTPUT VALUE 
^ CSP=MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 
TSP=MAXIMUM COMPLETION TIME 
USP=MAXI^nJM VALUE OF THE CONTROL VECTOR 
PARAM CI=1. ,CF=0. ,CSP=5. ,TSP= 1 . 00 ,USP=1 . 05 ,X2SP=1 . 0 
^ SET INITIAL INPUT VALUES ( Wl= DISTURBANCE , W2=REFERENCE SETTING) 
PARAM W1=0., W2=l., SEED=113, MEAN=1., SD=0.1, WL=157 

* SET INITIAL RUN FLAG AND EISPACK PARAMETERS 
PARAM 11=0, MATZ=0,NM=3,NN=3, 1=1, J=1,IER1=0,IERR=0 

INITIAL 

* CHECK FOR INITIAL RU^4 
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IF(II.GT.O) GOTO 200 



^ THE DESIGN VARIABLES ARE THE FEEDBACK GAINS 

^ INITIAL SYSTEM POLES TO BASE FEEDBACK GAINS UPON (COMPLEX FORM) 
P1=CMPLX( -10. 0,10.0) 

P2=CMPLX( -10.0, -10.0) 

P3=CMPLX( -300.0,0. ) 

* SET NON-ZERO COEFFICIENTS 
AA(1,2)=50. 

AA( 2,l)=-200. 

AA( 2,2)=-200. 

AA(3,1)=-1. 

B3(2)=200. 

cc(ia)=-i. 

H(2) = l. 

F(l,l)=-50. 

Q(l,l)=l. 

Q(2,2)=l. 

R(l) = l. 

^ ROUTINE TO FIND THE SYSTEM GAINS FROM KNOWN POLES (KUO SECTION 8-6) 
G( 1 )=REAL( -P1-P2-P3 )+AA( 1 ,1 )+AA( 2,2 )+AA( 3,3 ) 

G( 2 )=REAL( P1^P2+P1*P3+P2*P3 )-AA( 1,1 )^AA( 2,2)-... 

AA( 2,2)^AA(3,3)-AA( 1,1)*AA(3,3)+AA( 2,3)^AA(3,2)+AA(1,2)* . . . 
AA( 2,1)+AA(1,3)^AA(3,1) 

G(3 )=REAL(-P1^P2^P3)+AA( 1,1)*AA( 2,2)*AA(3,3)-AA( 1,1)*. . . 

AA( 2,3)*AA(3,2)+AA( 1,2)*AA(2,3)*AA(3,1)-AA(1,2)*AA(2,D*. . . 
AA(3,3)+AA( 1,3)*AA( 2,1 )*AA( 3 , 2 )-AA( 1 ,3 )*AA( 3 , 1 )*AA( 2 , 2 ) 
PA(1,1)=BB(1) 

PA(1,2)=BB(2) 

PA( 1,3)=BB(3) 

PA( 2,1 ) = ( -AA( 2,2 )-AA( 3,3 ) )*BB( 1 )+AA( 1,2 )*BB( 2 )+AA( 1,3 )*BB( 3 ) 
PA( 2,2)=AA(2,l)*BB(l)+(-AA(3,3)-AA(l,l))*BB( 2)+AA(2,3)*BB(3) 
PA(2,3)=AA(3,l)*3B(l)+AA(3,2)*BB(2) + (-AA(2,2)-AA(l,l) )*BB(3) 
PA(3,1) = (AA( 2,2 )*AA(3,3)-AA( 2,3)*AA(3,2) )*BB( 1)-(AA( 1,2)* . . . 
AA(3,3)-AA( 1,3 )*AA( 3 ,2 ) )*BB( 2 ) + ( AA( 1 ,2 )*AA( 2 ,3 )-AA( 1,3)* ... 
AA(2,2))*BB(3) 
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PA(3,2)=-( AA(2,1)*AA(3,3 J-AA( 2,3)^^AA(3,1) )*BB( l) + ( AA( 1,1 



AA( 3,3)-AA( 1,3 )^AA( 3 ,1 ) )^BB( 2 )-( AA( 1 ,1 )*AA( 2,3 )-AA( 1,3 )* . 
AA( 2,1) 

-A(3,3) = i AA( 2,1 j?«AA(3,2)-AA(3,1)^AA( 2,2 ) )*BB( 1 )-( AA( 1,1 
AA( 3*2 J-AA( 1,2)^AA( 3,1) )^BB( 2) + ( AA( 1,1 )*AA( 2,2 )-AA( 1 ,2 )* . 
AA( 2,1) )^BB(3) 

WRITE{6,297)Pi,?2,P3 

297 FORMAT! ‘0 INITIAL POLES ARE : ' ,3( ‘ ( ’ ,E11.9 , ' , ' ,E11 .9, ' J) ' 

CALL DGEFA(PA,3,3,IPVT,IER1) 

IF(IER1.NE.O)GOTO 300 
CALL 0GESL(PA,3,3,IPVT,G,0 ) 

WRITE! 6 , 299 )Gl 1 ) ,G( 2 ) ,G( 3 ) 

299 -FORMAT! *0 INITIAL GAINS G1 1 ) ,G( 2 ) ,G( 3 ) = *,3E15.9) 

WRITE! 6,290) 

290 FORMAT! ’0 STATE CONTROL MATRIX (B AB ) = ' ) 

WRITE! 6,298) BB( 1 ) , AA( 1 , 1 )*BB( 1 ) +AA( 1 , 2 )*BB( 2 ) 

WRITE! 6,298) BB! 2) ,AA! 2,1 )^BB! 1 )+AA! 2,2 )^BB! 2 ) 

298 FORMAT! 2E 15. 4) 

WRITE! 6,291) 

291 FORMAT! ’0 OUTPUT CONTROL MATRIX ! CB CAB ) = ' ) 

WRITE! 6,298) CC! 1,1 )^^BB! 1 ) +CC! 1 ,2 )^BB! 2 ) ,CC! 1 , 1 )3^! AA! 1,1)* 
BB! D+AA! 1,2 )*BB! 2 ) )+CC! 1,2 )*! AA! 2 , 1 )*BB! 1 )+AA! 2,2)*B3! 2) ) 
WRITE! 6,292) 

292 FORMAT! '0 OBSERVABILITY MATRIX ! C* A*C* ) = ’) 

WRITE! 6,298) CC! 1 , 1 ) ,AA! 1 ,1 )*CC! 1 , 1 )+AA! 2 , 1 )*CC! 1 , 2 ) 

WRITE! 6,298) CC! 1,2 ) ,AA! 1, 2 )*CC! 1 ,1 )+AA! 2,2)*CC! 1,2) 

GOTO 302 

300 WRITE! 6,301) lERl 

301 FORMAT! '0 POLE SELECTION ERROR lERl = ’,17) 

CALL ENDJOB 

302 CONTINUE 

DV!1)=G!1) 

DV! 2)=G!2) 

DV!3 )=G!3 ) 

IF! INFO.EQ.O ) GOTO 200 
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CALL ADS(INFO,ISTRAT,IOPT,IONED,IPRINT,IGRAD,NDV,NCON,DV,VLB 
VUB , OBJ ,CON , IDG ,NGT , IC , DF , WA , NRA ,NCOLA ,NRWK , IWK ,NRIWK ) 
IWK( 2)=0 
WRITE! 6,220) 

:C0 FORMAT! ’0 3CALI^4G MOT USED* ) 

* CALL TO ADS TO GET THE DESIGN VARIABLE FOR THE SYSTEM 

200 CALL ADS( INFO, ISTRAT,IOPT,IONED,IPRINT,IGRAD,NDV,NCON,DV, VLB,. 
VUB , OBJ ,CON , IDG ,NGT , IC , DF ,WA ,NRA ,NCOLA ,WK ,NRWK ,IWK ,NRIWK ) 

G( 1)=DV! 1) 

G( 2)=DV( 2) 

G(3)=DV(3) 

^ CHECK FOR TYPE OF RUN IN ORDER TO SAVE INITIAL AND FINAL RESULTS 
IFdl.iNE.O )30T0 235 

GRAPH ! INITIAL/RUN, RU=1,D£=TEK618 ) TIME,X1,X2,U 
GOTO 245 

235 IF(INFO.NE.O)GOTO 240 
DELPRT=0.01 
DELPLT=0.01 

GRAPH ( FINAL/RUN, RU=^,DE=TEK618 ) TIME,X1,X2,U 
GOTO 245 
240 DELPRT=0.0 

DELPLT=0,0 
24S 11=1 

* INITIALIZE CONSTRAINT CHECK VALUES 
LAST THREE VALUES OF THE OUTPUT 

Y0SC1=CI 

YOSC2=CI 

YOSC3=CI 

^ MAXIMUM VALUE OF THE OUTPUT 
CMAX=0. 

^ THE TIME CONSTANT OF THE OUTPUT 
TAU=D. 

^ TIME WHEN COMPLETION TIME SPECIFICATION IS MET 
TC=0. 

* INITIAL VALUE OF THE CONTROL VECTOR 
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UI=0 



* MAXIMUM VALUE OF THE CONTROL 

UMAX=0. 

^ INITIAL VALUE OF THE X2 STATE 
XJI=0. 

* MAXIMUM VALUE OF THE X2 STATE 

X2MAX=0. 

* VALUE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 

LSTOBJ=0. 

^ SLOPE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 
LSLOPE=0. 

^ VALUE OF TIME AT END OF LAST INTEGRATION INTERVAL 
LSTTIM=0. 

DERIVATIVE 

* PERFORM INTEGRATIONS TO FIND NEW STATES 

X1=INTGRL(0, ,XX1) 

X2=INTGRL(0. ,XX2 ) 

XX1=AA( 1,1 )9^X1+AA( 1,2 )*X2+BB( 1 )*U + F( 1,1 )9^W1 + F( 1,2 )^^W2 
XX2=AA( 2,1 )*X1+AA( 2,2 }*X2+BB( 2 )*U+F( 2,1 )^^W1 + F( 2 ,2 )^W2 

* FIND NEW CONTROL VALUE 

U=-G( 1 )*X1-G( 2 )*X2-G( 3 )^INTGRL( 0 . ,C ) 

* FIND NEW OUTPUT VALUE 

C=CC( 1,1 )*X1+CC( 1 ,2 )*X2+H( 1 )*W1+H( 2 )*W2 
FIND NEW OBJECTIVE FUNCTION VALUE USING QUADRATIC PERFORMANCE INDEX 

* AFTER TRANSLATION OF OUTPUT AND CONTROL VALUES 

X1T=X1-1, 

UT=U-1. 

OBJ=INTGRL(0. ,OOBJ) 

003J=(X1T^Q( 1,1 )+X2^^Q( 1,2) )*XlT + ( X1T^^Q( 2 , 1 )+X2^^Q( 2 , 2 ) )*X2+ . , . 
UT9^R( 1 ywj 

DYNAMIC 

* DISTURBANCE FUNCTION 

* Wl=NORMAL( SEED, MEAN, SD ) 
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^ H1 = 1.+0.1^SIN(WL^TIME ) 

^ CONSTRAINT CHECKS AT EACH INTEGRATION STEP 
^ ESTABLISH MAX OUTPUT OSCILLATION CHECKS ( YOSC ) 

^ STEP FORWARD THE LAST THREE OUTPUT VALUES 
YOSCl=YOSC2 
Y0SC2=Y0SC3 
Y0SC3=C 

^ CALCULATE THE DIFFERENCE OF THE ACTUAL OUTPUTS TO THE DESIRED OUTPUTS 
H1=H2 
H2=H3 

H3=ABS(Y0SC3-CF) 

^ CALCULATE THE DIFFERENCE OF THE ACTUAL OUTPUT TO THE INITIAL OUTPUT 
DELC=ABS(CI-C) 

^ CALCULATE THE DIFFERENCE OF THE ACTUAL CONTROL TO THE INITIAL CONTROL 
DELU=ABS(UI-U) 

^ FIND MAXIMUM CONTROL VALUE 

I F ( DELU . GE . UMAX )UMAX=DELU 

^ CALCULATE THE DIFFERENCE OF THE ACTUAL CONTROL TO THE INITIAL CONTROL 
DELX2=ABS(X2I-X2) 

^ FIND MAXIMUM X2 STATE VALUE 

IF( DELX2 . GE . X2MAX )X2MAX=DELX2 
^ ESTABLISH TIME CONSTANT CHECK (TAU) 

IF(TAU.NE.O. )GOTO 20 

I F ( DELC . GE . ABS( 0 . 63212* ( Cl -CF ) ) )TAU=TIME 
20 IF(DELC.GE.CMAX) CMAX=DELC 

* MAX AMPLITUDE CONSTRAINT (CON(D) 

* (+) IF ABOVE SPECIFICATION (-) IF BELOW SPECIFICATION 

IF( DELC.EQ.CMAX) CON( 1 )=DELC-ABS( CI-CF )*( CSP/100 + 1 . ) 

* MAX COMPLETION TIME CONSTRAINT (C0N(2)) 

* SKIP CHECK IF COMPLETION TIME ALREADY MET 

IF(TC.NE.O. )GOTO 150 

* FOR UNDER DAMPED CASE 

* CHECK FOR A OSCILLATION REVERSAL, AND IF THE PEAK MAGNITUDE IS LESS 

* THAN IX OF THE DIFFERENCE BETWEEN THE INITIAL AND FINAL OUTPUT VALUES 

IF( (H2.GE*H1).AND.(H2.GE.H3) ... 
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•AND.(H2.LE.ABS( (CI-CF)*.01))) TC=TIME 
CHECK FCR THE LAST THREE OUTPUT VALUES TO BE DECREASING AND ALL LESS 
THAN V/. OF THE DIFFERENCE BETWEEN THE INITIAL AND FINAL OUTPUT VALUES 
:F( ;H3.LT.rf2).AND. ( HE.LT.HD.AND. (Hl.LT. . . . 
A6G((CI-CF)*.01))) TC=TIME 
FOR OVER DAMPED CASE 

* CHECK FOR THE CLTTPUT AT 5 TIMES THE TIME CONSTANT TO BE LESS 

THAN 17. OF THE DIFFERENCE BETWEEN THE INITIAL AND FINAL OUTPUT VALUES 
IF( (C.GT. (CF-ABS(CI-CF]^(.01 )). AND .( TIME .GE .( S. ^fTAU ) ) . . . 
.AND.(TAU.NE.O. ) . AND . ( C . EQ . CMAX n TC=TIME 
^ IF COMPLETION TIME NOT MET SET CONSTRAINT POSITIVE 
CCN( 2 )=10. 

SET CONSTRAINT VALUE NEGATIVE 3Y VALUE FROM TIME SPECIFICATION 
IFdC.NE.O.) CON( 2)=TC-TSP 

* MAX ALLOWABLE CONTROL VECTOR CONSTRAINT ( CON( 3 ) ) 

( + ) IF ABOVE SPECIFICATION; (-) IF BELOW SPECIFICATION 
150 IF(DELU.EQ.UMAX) CON( 3 )=DELU-USP 

* UNSTABLE CONSTRAINT CHECK (C0N(4)) 

* (+1) IF OUTPUT IS TEN TIMES THE MAX ALLOWABLE (-1) OTHERWISE 

CON( <♦)=-!. 

IF(ABS(C ).GT.ABS(CI-CF)i«10.it(CSP/100..+l. ) )C0N(4)=1. 
if TERMINATE RUN EARLY IF UNSTABLE TO PREVENT "BLOWING UP" 
IF(CON(<^).GT.O. ) CALL ENDRUN 

* NEGATIVE RESPONSE CONSTRAINT CHECK ( CON( 5 ) ) 

(+1) IF OUTPUT INITIALLY GOES NEGATIVE (-1) OTHERWISE 
C0N(5)=-1. 

IF( ( ABS( CF-C ) .GT . ABS( CF-CI ) ) . AND . ( TAU . EQ . 0 . ) )CON( 5 )=1 . 

^ TERMINATE RUN EARLY IF UNSTABLE TO PREVENT "BLOWING UP" 
IF(CON(5).GT.O. ) CALL ENDRUN 
^ MAX ALLOWABLE X2 STATE CONSTRAINT (C0N(6)) 

( + ) IF ABOVE SPECIFICATION (-) IF BELOW SPECIFICATION 
IF(DELX2.EQ.X2MAX) CON( 6 )=DELX2-X2SP 
INTEGRATION COMPLETE CHECK 

30 IF((OBJ.LE.O. ).0R. ( LSTOBJ.LE.O. ) )GOTO 40 

* FOR OBJECTIVE FUNCTIONS THAT APPROACH A STEADY STATE VALUE CHECK FOR 
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* A OBJECTIVE FUNCTION SLOPE OF LESS THAN .IX OTHERWISE CHECK FOR THE 



* SECOND DERIVATIVE OF THE OBJECTIVE ’ FUNCTION TO BE LESS THAN . IX 
^ ALONG WITH COMPLETION TIME BEING REACHED IN EITHER CASE 

SLOPE = ( OBJ-LSTOBJ )/i TIME-LSTTIM ) 

ISlCPE^f SLOPE-LSLCPE )/( TIME-LSTTIM ) 

IF( ( (SLOPE.lt. .001).OR.(DSLOPE.LE. .001)). . . 

.AND.CTC.NE.O. )) CALL ENDRUN 

ADVANCE THE LAST FUNCTION, TIME AND SLOPE COUNTERS FOR NEXT INTERVAL 
40 LSTOBJ=OBJ 

LSTTIM=TIME 

LSLOPE=SLOPE 

TERMINAL 

IF (INFO.NE.O) GOTO 450 
DO 4C0 1=1,3 
00 410 J=l,3 

410 V(I,J)=AA(I,J)-BB(I)*G( J) 

400 CONTINUE 

WRITE( 6 ,499 )G( 1 ) ,G( 2 ) ,G( 3 ) 

499 FORMAT( ‘0 FINAL GAINS G( 1 ) ,G( 2 ) ,G( 3 ) = ’,3E1S.4) 

* CALL TO FIND THE EIGNVALES WHICH ARE IN TURN THE FINAL SYSTEM POLES 

CALL RG(NM,NN,V,WR,WI,MATZ,Z,IV1,FV1,IERR) 
IF(IERR.NE.0)WRITE(6,495)IERR 

495 FORMAT( ‘0 ERROR FROM EISPACK ‘,13 ,* -FINAL POLES MAYBE WRONG*) 
WRITE(6,497)( WR(I ),WI(I),I = 1,3) 

497 FORMAT(‘0 FINAL POLES ARE : * ,3( ‘ ( ‘ ,E11 .4, * , * ,E11 .4, ‘ J) *)) 
CALL ENDJOB 
450 CALL RERUN 
END 
STOP 
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APPENDIX E 

PROGRAM LISTINGS FOR POLE TO GAIN APPROXIMATIONS 



1. DETERMINANT METHOD 

The following listing would replace the initialization and initial segments of the 
DSL program listed in Appendix D, for using the determinant evaluation method to 
find the initial gains corresponding to desired initial system poles. The DERIVATIVE, 
DYNAMIC, and TER.MINAL sections of the state variable program would remain 
unchanged. 



TITLE DC MOTOR INTEGRAL CONTROL SYSTEM-KUO EXAMPLE 8-7 

INITIALIZE THE ADS PARAMETERS 
FIXED NRA,NGT,NCOLA,NRWK,NRIWK,IGRAD,NDV,NCON,IDG,IC,IWK 
FIXED ISTRAT,IOPT,IONED,IPRINT,INFO 

ARRAY DV( 7 ) ,VL8( 7 ) ,VU3( 7 ) ,CON( S 1 ,IDG( S ) ,IC( S ) ,DF( 7 ) ,WK( 1000 ) ,IWK( 1000 ) 
D DIMENSION WA( 40,40) ' 

PARAM NRA=40 ,NCOLA=40 ,NRWK=1000 ,NRIWK=1000 ,NDV=5 ,NCON=l ,IGRAD=0 ,INFO=0 

* ESTABLISH TYPES OF CONSTRAINTS (ADS MANUAL TABLE 5) 

TABLE IDG(1-S)=3^-1,2*2 

* LOWER BOUND FOR DESIGN VARIABLES 
TABLE VLB( 1-6 )=3»-1000., 35^0, 

UPPER BOUND FOR DESIGN VARIABLES 
TABLE VUB( 1-6 )=5^^1000. ,3^(0. 

ADS OPTIMIZATION METHOD SELECTION 
PARAM ISTRAT=9, lOPTrS, I0NED=7, IPRINT=3140 

* INITIALIZE THE MODEL PARAMETERS 

* SET CONTROL PARAMETERS 

FIXED I,J,K,II,IER1,IERR,NN,NM,IV1,MATZ,SEED 
ARRAY BB(3 ),H(2),R(2),G(3 ),WR(3 ),WI(3 ),IV1(3 ),FV1(3) 

D DIMENSION AA( 3 ,3 ) ,CC( 2 , 2 ) ,Q( 2 ,2 ) ,F( 2, 2 ) ,V( 3 ,3 ) ,Z( 3 ,3 ) 
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SET INITIAL PARAMETERS FOR POLE TO GAIN OETERMINATION 



TABLE Gt l-2)=-.38,.6,-6. 

0 INTEGER NP>JOB,IPVT(3),INF 

0 COMPLEX9616 Pi 3 ) ,PA( 3 ,3 ) ,PM( 3,3 ) ,?DET( 2 ) ,PWK( 3 ) 

EXCLUDE NP , JOB ,IPVT ,INF ,P ,PA ,PM,POET ,PWK 
96 ZERO ALL MATRICES 
0 DATA AA,F,CC,Q/21960,/ 

TABLE BBC 1«3)=3960, ,R(l-2) = 2*0. ,H(l-2)=2^0. 

96 SET CONSTRAINT SPECIFICATIONS 

96 CI=INITIAL OUTPUT VALUE j CF=FINAL EXPECTED OUTPUT VALUE 
96 CSP=MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 
96 TSP=MAXIMUM COMPLETION TIME 
96 US?=MAXIMUM VALUE OF THE CONTROL VECTOR 
PARAM CI=1. ,CF=D. ,CSP=4D. ,TSP=1 .5 ,USP=1D . 

96 SET INITIAL INPUT VALUES ( W1=DISTURBANCE , W2=REFERENCE SETTING) 

PARAM W1=D., W2=l., SEED=113, MEAN=1., SD=D.l, WL=1.6 
96 SET INITIAL RUN FLAG AND EISPACK PARAMETERS 
PARAM II=D,MATZ=D,NM=3,NN=3,I=1,J=1,K=1,IER1=D,IERR=0 

96 ***^ 96 * 9696*96969696969696969696969696 96969696969696 969696 96 ^** 969696 ** 969696969696 ** 969696 * 96969696 ** 9696****96969696 96 

96 INITIALIZE THE DSL PARAMETERS 
PRINT 03J,C,X1,X2,U 
SAVE (RUN) OBJ,Xl,X2,U,C 
CONTROL DELPLT=D.DD5,DELPRT=0.D1 
CONTROL FINTIM=4. 

METHOD STIFF 
RELERR 0BJ=l.E-4 
ABSERR OBJ=l.E-<h 

*********************************************************************** 

INITIAL 

* CHECK FOR INITIAL RUN 

IF(II.GT.O) GOTO 20D 

* THE DESIGN VARIABLES ARE THE FEEDBACK GAINS 

* INITIAL SYSTEM POLES TO BASE FEEDBACK GAINS UPON (COMPLEX FORM) 

P(1)=CMPLX(-1D.D,1D.D) 

P(2)=CMPLX( -1D.D,-1D.D) 
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P(3 )=CMPLX( -300.0,0. ) 

* SET NON-ZERO COEFFICIENTS 

AA( 1»2)=50. 

AAt ca)=-zco. 

AA( C,2 ) = -Z00. 

AA( 3,1)=-1. 

S3( 2)=200. 

CC(1,1)=-1. 

H(2) = l, 

F( l,l)=-50. 

1,1 ) = i. 

Qi 2>2 1 = 1. 

R(l) = l. 

^ SET THE DESIGN VARIABLE MATRIX 
1 ) 

DV( 2i=G( 2) 

DV(3 )=G( 3) 

DV(^)=Q( 1,1) 

DV(5)=q(2,2) 

DV(6)=R(1) 

* ROUTINE TO FIND THE SYSTEM GAINS FROM KNOHN POLES USING ADS 

NP=3 

JOB=10 

INFO=0 

^ USE ADS TO FIND G*S THAT MAKE SUM( DET( RELATION MATRIX)) = 0. 

660 CALL ADS( INFO, ISTRAT,IOPT,IONED,IPRINT,IGRAD, NOV, NCON,DV, VLB, . 
VUB,03J,C0N,IDG,NGT,IC,DF,WA,NRA,NC0LA,HK,NRWK,IWK,NRIWK) 
IF(INFO.EQ.O)GOTO 650 
DET=0. 

* CALCULATE THE RELATION MATRIX: (P-A+B?fG) 

* DO 600 K=1,NP 

DO 610 1=1, NP 
DO 620 J=1,NP 
620 PA(I,J) = (-AA(I,J)) 

610 PA(I,I)=P(I)-AA(I,I) 
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DO 630 1=1, NP 



G(I)=DV(I) 

00 640 J=1,NP 

640 PM( J,I )=PA( J,I )+BB( J )^^G( I ) 

630 CONTIN'UE 

T :aLL ’0 CONDITION THE RELATION MATRIX 
CALL ZG£FA(PM,NP,NP,IPVT,INF) 

* CALL TO FIND OERTERMINATE OF RELATION MATRIX 
CALL 2GEDK PH,NP,NP,IPVT,PDET,PWK,JOS) 

FIND OBJECTIVE FUNCTION-LOOKING FOR MINIMUM TO BE ZERO 
CON( K )=ABS( REAL( PDET( 1 )*1D . >e^POET( 2 ) ) ) 

OBJ=OBJ+CON( K ) 

600 CONTINUE 

GOTO 660 

^ OUTPUT INITIAL POLE AND GAINS INFORMATION 
^50 >JRITE1d,297)P(1),P( 2),P(3) 

297 FORMAT('D DESIRED POLES WERE : ’ ,3( ’ ( ' ,E11.4, ’ , • ,E11,4, ’ J) •)) 

WRITE( 6,299 )G( 1 ) ,G( 2 ) ,G( 3 ) 

299 FORMAT( '0 INITIAL GAINS G( 1) ,G( 2 ) ,G( 3 ) = ’,3E1S,4) 

CALL TO FIND THE EIGNVALES WHICH ARE IN TURN THE INITIAL SYSTEM POLES 
DO 680 1=1,3 
DO 690 J=l,3 

690 V(I,J)=AA(I,J)-BB(I)^«G(J) 

680 CONTINUE 

CALL RG(NM,NN,V,WR,WI,MATZ,Z,IV1,FV1,IERR) 

IF( lERR.NE . 0 )WRITE( 6 ,485 )IERR 

485 FORMAT(’D ERROR FROM EISPACK ' ,13 ,' -INITIAL POLES MAYBE WRONG’) 
WRITE(6,487)(WR(I ),WI(I ),I = 1,3) 

487 FORMAT( '0 INITIAL POLES ARE : ' ,3( ' ( ’ ,E11.4, ’ , ’ ,E11.4, ' J) ')) 
WRITE( 6,290) 

290 FORMATC '0 STATE CONTROL MATRIX (B AB ) = ’ ) 

WRITE( 6,298 ) BB( 1 ) ,AA( 1 ,1 )^BB( 1 )+AA( 2 ) 

WRITE(6,298) BB( 2 ) , AA( 2 ,1 )^BB( 1 )+AA( 2 , 2 )^BB( 2 ) 

298 FORMAT( 2E15.4) 

WRITE (6,291) 
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291 FORMAK ‘0 OUTPUT CONTROL MATRIX ( CB CAB) = *) 

WRITE( 6,298) CC( 1,1 )hBB( 1 )+CC( 1,2 )^BB( 2 ) ,CC( 1,1 )^( AA( 1,1 )* 
BB( 1)+AA( 1,2)*BB( 2) )+CC(l,2)^( AA( 2,1)^BB( 1)+AA( 2,2 )*BB( 2) ) 
WRITE(6,292) 

292 FORHAK'O OBSERVABILITY MATRIX ( C^ A*C^) = •) 

WRITE( 6,290 CC( 1 ,1 ) , AA( 1 , 1 )*CC( 1 ,1 )+AA( 2,1)*CC( 1,2) 

K.-iITE( 6,298) CC( 1 , 2 ) ,AA( 1 ,2 )^CC( 1 ,1 )+AA( 2 , 2 )^CC( 1 , 2 ) 

^ RESET AOS PARAMETERS FOR MODEL OPTIMIZATION 
CALL ENDJOB 
NC0N=5 
ISTRAT=0 
IOPT=^f 
I0NED=7 
IPRINT=1010 
IOG( 1)=0 
I0G(2)=0 • 

I0G(3)=0 

CALL TO AOS TO GET THE DESIGN VARIABLE FOR THE SYSTEM 
200 CALL AOS( INFO, ISTRAT,IOPT,IONED,IPRINT,IGRAD, NOV ,NCON,OV, VLB, 
VUB ,OBJ ,CON ,IOG ,NGT , IC , OF ,NA ,NRA ,NCOLA ,WK ,NRWK ,IWK ,NRIWK ) 
G(1)=0V(1) 

G(2)=0V( 2) 

G(3)=0V(3) 

Q( 1,1)=0V(4) 

Q(2,2)=0V(5) 

R(1)=0V(6) 

CHECK FOR TYPE OF RUN IN ORDER TO SAVE INITIAL AND FINAL RESULTS 
IF(II.NE.O)GOTO 235 

GRAPH ( INITIAL/RUN, RU=1,0E=TEK618 ) TIME,X1,X2,U 
GOTO 2^5 

235 IF(INFO.NE.O)GOTO 240 
0ELPRT=0.01 
0ELPLT=0.01 

GRAPH ( FINAL/RUN, RU=^f,DE=TEK618 ) TINE,X1,X2,U 
GOTO 245 
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2^0 DELPRT=0,0 



DELPLT=0.0 
245 11=1 

^ INITIALIZE CONSTRAINT CHECK VALUES 
^ LAST THREE VALUES OF THE OUTPUT 
Y0SC1=CI 
Y0SC2=CI 
Y0SC3=CI 

^ MAXIMUM VALUE OF THE OUTPUT 
CMAX=0 . 

^ THE TIME CONSTANT OF THE OUTPUT 
TAU=0. 

* TIME WHEN COMPLETION TIME SPECIFICATION IS MET 

TC=0. 

^ INITIAL VALUE OF THE CONTROL VECTOR 
UI=0. 

* MAXIMUM VALUE OF THE CONTROL 

UMAX=0. 

* VALUE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 

LST0BJ=0. 

* SLOPE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL* 

LSL0PE=0. 

^ VALUE OF TIME AT END OF LAST INTEGRATION INTERVAL 
LSTTIM=0. 



2 . EIGENVALUE METHOD 

The following program listing would be used if the the eigenvalue approximation 
method was desired for finding the gains corresponding to the desired initial system 
poles. The DERIVATIVE, DYNAMIC, and TERMINAL sections of the state 
variable program in Appendix D remain unchanged, and would follow the listing 
below. 

TITLE DC MOTOR INTEGRAL CONTROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 
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INITIALIZE THE ADS PARAMETERS 



FIXED NRA ,NGT ,NCOLA ,NRWK ,NRIWK , IGRAD , NDV ,NCON , IDG , IC , IHK 
FIXED ISTRAT,:OPT,IQNED.IPRINT,INFO 

ARRAY 0V( 7 1 ,VLS( 7 • .VUBt 7 ) ,CCN( 5 ) ,IDG( S ) , IC( 5 ) ,DF( 7 ) ,WK( 1000 ) ,IWK( 1000 ) 

0 OIHENSION WA( 40,40) 

PARAM NRA=40,NCOLA=40,NRHK=1000,NRIWK=1000,NDV=3,NCON=4,IGRAD=0,INFO=-2 
ESTABLISH TYPES OF CONSTRAINTS (ADS MANUAL TABLE 5) 

TABLE IDG(1-S)=45«-1,2 

LONER BOUND FOR DESIGN VARIABLES 
TABLE VLB( 1-6 )=3?^-1000 . ,3^0 . 

UPPER BOUND FOR DESIGN VARIABLES 
TABLE VUB( 1-6 J=3J^1000. ,3>0. 

^ ADS OPTIMIZATION METHOD SELECTION 

PARAM ISTRAT=1, I0PT=3, ICNED=3, IPRIMT=1100 

INITIALIZE THE MODEL PARAMETERS 
^ SET CONTROL PARAMETERS 

FIXED I,J,K,II,IER1,IERR,NN,NM,IV1,MATZ,SEED 
ARRAY B3(3),H(2),R(P),G(3),NR(3),WI(3),IV1(3),FV1(3) 

D DIMENSION AA( 3 ,3 ) ,CC( 2 , 2 ) ,q( 2 ,2 ) ,F( 2 ,2 ) ,V( 3 ,3 ) ,Z( 3 ,3 ) 

D C0MPLEXX16 P(3) 

EXCLUDE P 

SET INITIAL PARAMETERS FOR POLE TO GAIN DETERMINATION 
TABLE G(l-3)=35«-l. 

^ ZERO ALL MATRICES 
D DATA AA,F,CC,Q/21^0./ 

TABLE BB(l-3)=3»(0.,R(l-2) = 2»0. ,H( 1-2 )=2?^0 . 

SET CONSTRAINT SPECIFICATIONS 
^ CI=INITIAL OUTPUT VALUE 5 CF = FINAL EXPECTED OUTPUT VALUE 
* CSP=MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 
TSP=MAXIMUM COMPLETION TIME 
USP=MAXIMUM VALUE OF THE CONTROL VECTOR 
PARAM CI=1.,CF=0. ,CSP=100. ,TSP=1. 9,USP=1. 

^ SET INITIAL INPUT VALUES ( H1=DISTURSANCE , W2=REFERENCE SETTING) 

PARAM H1=0., H2=l., SEED=113, MEAN=1. , SD=0.1, WL=1.6 
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* SET INITIAL RUN FLAG AND EISPACK PARAMETERS 



PARAM 11=0, MATZ=0,NM=3,NN=3 >1=1, J=1,K=1,IER1=0,IERR=0 

^ INITIALIZE THE DSL PARAMETERS 
^PINT CBJ,C,X1,X2,U 
SAVE (RUN) 0BJ>X1,X2,U,C 
CONTROL 0ELPLT=0 . 005 , CE LPRT=0 . 01 
CONTROL FINTIM=^. 

METHOD STIFF 
RELERR 0BJ=1.E-^ 

ABSERR 03J=l.E-4 

INITIAL 

^ CHECK FOR INITIAL RUN 

IF(II.GT.O) GOTO 200 

THE DESIGN VARIABLES ARE THE FEEDBACK GAINS 
^ INITIAL SYSTEM POLES TO BASE FEEDBACK GAINS UPON (COMPLEX FORM) 
P ( 1 )=CMPLX( -100 .0,10.0) 

P( 2)=CMPLX( -100.0,-10.0 ) 

P(3)=CMPLX( -300.0,0. ) 

SET NON-ZERO COEFFICIENTS 
AA( 1,2)=50. 

AA( 2,l)=-200. 

AA( 2,2)=-200. 

AA(3,1)=-1. 

BB(2)=200." 

CC(1,1)=-1. 

H(2)=l. 

F( l,l)=-50. 

Q(l,l)=l. 

Q( 2,2)=1. 

R(l)=l. 

* SET THE DESIGN VARIABLE MATRIX 
DV(1)=G(1) 

DV( 2)=G( 2) 
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DV(2)=G13) 



DV(4)=Q( 1,1) 

DV(5)=Q(2,2) 

DV(6)=R(1) 

CALL ADS( INFO, ISTRAT,IOPT,IONED,IPRINT,IGRAD,NDV,NCON,OV, VLB 
VUB ,OBJ ,CON ,IDG ,NGT ,IC ,DF , WA ,NRA ,NCOLA ,WK ,NRWK ,IWK ,NRIHK ) 
IWK( 2)=D 

* ROUTINE TO FIND THE SYSTEM GAINS FROM KNOWN POLES USING ADS 

66D CALL ADS( INFO, ISTRAT,IOPT,IONED,IPRINT,IGRAD,NDV,NCON,DV, VLB,. 
VUB,OBJ,CON,IDG,NGT,IC,DF,WA,NRA,NCOLA,WK,NRWK,IWK,NRIWK) 
IF(INFO.EQ.D)GOTO 650 

* CALCULATE THE RELATION MATRIX: (A-B^G) 

DO 7DD 1=1,3 
G(I)=DV(I) 

DO 71D J=l,3 

71D V(J,I)=AA(J,I)-B8( J)5^G(I) 

7D0 CONTINUE 

if CALL TO FIND THE RELATION MATRIX EIGENVALUES 

CALL RG(NM,NN,V,WR,WI,MATZ,Z,IV1,FV1,IERR) 

IF( lERR .NE . D )WRITE (6,791 )IERR 

791 FORMAK '0 ERROR FROM EISPACK ' ,13 ,' -FINAL POLES MAYBE WRONG') 
IF(MAX(WR(l),WR(2),WR(3)).GT.O. )G(1)=-G(1) 

CON( 1 )=ABS( MAX( WR( 1 ) ,WR( 2 ) ,WR( 3 ) )-MAX( REAL( P( 1 ) ) , . . . 
REAL(P(2)),REAL(P(3)))) 

CON( 2 )=ABS( MIN( WR( 1 ) ,WR( 2 ) ,WR( 3 ) )-MIN( REAL( P( 1 ) ) , . . . 
REAL(P(2)),REAL(P(3)))) 

CON( 3 )=ABS( MAX( WI( 1 ) ,WI ( 2 ) ,WI ( 3 ) )-MAX( IMAG( P( 1 ) ) , . . . 
IMAG(P(2)),IMAG(P(3)))) 

CON(^)=ABS(MIN(WI( 1),WI( 2),WI( 3))-MIN(IMAG(P(l) ), . . . 
IMAG(P(2)),IMAG(P(3)))) 

OBJ=CON( 1 )+CON( 2 )+CON( 3 )+CON( ^ ) 

GOTO 660 

* RESET ADS PARAMETERS FOR PERFOMRANCE INDEX OPTIMIZATION 

65D NC0N=5 

IDG(^)=2 
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HRITE( 6,296 )P( 1 ) , P( 2 ) ,P( 3 ) 



296 FORMAK'O DESIRED POLES ARE : * ,3( * ( * ,E11 * ,E11 * J) *)) 

WRITE( 6,297 )WR(1 ),WI( 1 ),WR( 2 ) ,HI( 2 ) , WR( 3 ) ,WI( 3 ) 

297 FORMAK'O INITIAL POLES ARE : ' ,3( ' ( ' ,E11 . ' ,E11 ' J) ’)) 

WRITE( 6 , 299 )G( 1 ) ,G( 2 ) ,G( 3 ) 

299 FORMAK'O INITIAL GAINS G( 1 ) ,G( 2 ) ,G( 3 ) = ',3E15.^) 

WRITE( 6,290) 

CALL ENDJOB 

290 FORMAK'O STATE CONTROL MATRIX (B AB ) = ') 

WRITE( 6,298 ) BB( 1 ) ,AA( 1 ,1 )^BB( 1 )+AA( 1 ,2 )^BB( 2 ) 

WRITE( 6,298) BB( 2 ) ,AA( 2 , 1 )^BB( 1 )+AA( 2 , 2 )^BB( 2 ) 

298 F0RMAT(2E1S.^) 

WRITE(6,291) 

291 FORMAK'O OUTPUT CONTROL MATRIX ( CB CAB) = ') 

WRITE(6,298) CC( 1 , 1 )*BB( 1 )+CC( 1 , 2 )^BB( 2 ) ,CC( 1 , 1 )^( AA( 1 , 1 )^^ ... 
BB( 1 )+AA( 1 ,2 )5^BB( 2 ) )+CC( 1,2 )^( AA( 2,1 )^BB( 1 )+AA( 2 , 2 )^BB( 2 ) ) 
WRITE( 6,292) 

292 FORMAK'O OBSERVABILITY MATRIX ( C^^ A^^C^^ ) = ') 

HRITE( 6 ,298 ) CC( 1, 1 ) ,AA( 1 ,1 )^CC( 1 ,1 )+AA( 2 ,1 )^CC( 1,2) 

HRITE( 6,298) CC( 1 , 2 ) , AA( 1 , 2 )^CC( 1 ,1 )+AA( 2,2)^CC( 1,2 ) 

^ CALL TO ADS TO GET THE DESIGN VARIABLE FOR THE SYSTEM 

200 CALL ADS( INFO, ISTRAT,IOPT,IONED,IPRINT,IGRAD,NDV,NCON,DV, VLB, . . . 
VUB ,OBJ , CON , IDG ,NGT ,IC , DF , WA ,NRA ,NCOLA ,WK , NRWK , INK ,NRIWK ) 
G(1)=DV(1) 

G(2)=DV(2) 

G(3)=DV(3) 

Q(1,1)=DV(4) 

Q( 2,2)=DV(5) 

R(1)=DV(6) 

3^ CHECK FOR TYPE OF RUN IN ORDER TO SAVE INITIAL AND FINAL RESULTS 
IF(II.NE.O)GOTO 235 

GRAPH ( INITIAL/RUN, RU=1,DE=TEK618 ) TIME,X1,X2,U 
GOTO 295 

235 IF(IMFO.NE.O)GOTO 290 
DELPR'T=0.01 
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OELPLT=0.01 



GRAPH ( FINAL/RUN, RU=^,0E=TEK618 ) TIME,X1,X2,U 
GOTO 2^5 
2^0 OCLPRT=0.0 

G£LPLT=0. J 
Z^5 11=1 

* INITIALIZE CONSTRAINT CHECK VALUES 

* LAST THREE VALUES OF THE OUTPUT 

Y0SC1=CI 

Y0SC2=CI 

Y0SC3=CI 

^ MAXINUN VALUE OF THE OUTPUT 
CMAX=0. 

^ THE TIME CONSTANT CF THE OUTPUT 
TAU=0. 

* TIME WHEN COMPLETION TIME SPECIFICATION IS MET 

TC=0, 

* INITIAL VALUE OF THE CONTROL VECTOR 

UI=0. 

* MAXIMUM VALUE OF THE CONTROL 

UMAX=0 . 

* VALUE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 

LSTOBJ=0, 

* SLOPE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION INTERVAL 

LSLOPE=0, 

^ VALUE OF TIME AT END OF LAST INTEGRATION INTERVAL 
LSTTIM=0. 
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APPENDIX F 

CASE 2 RESULTS FOR SCALING 
RESULTS FOR STRATEGY 057 WITH SCALING APPLIED 



DSL/VS RELEASE I MODIFICATION LEVEL 1 19B5 



>>>> DSL SIMULATION INPUT DATA <<<< 

TITLE DC MOTOR INTEGRAL CONTROL SYSTEM-NO NDISE-KUD EXAMPLE 8-7 
FIXED NRA.NGT.NCDLA.NRWK.NRIWK, IGRAD .NDV.NCCN, IDG, IC, IWK 
FIXED ISTRATr ICPT, lONED. IPRINT. INFO, IPVT 

oaRAM NRA=^iO.NCDLA=AO ,NRWK=1000 ,NRIWK=lOOO,NDV=3,NCON=5, IGRAD=0 , INFO=0 
TABLE :DG(1-5) = 3»'0.2'<2 
'■ABLE 7L3C :-6) = 3>*-i0Q0. ,3-0. 

TABLE VLB( l-ols3-l000. .3-0 . 

’ARAM :STRAT=0. ;OPT=5. IDNED=7. IPRINT=3050 
SAVE C3J,Xl,X2,U 

GRAPH (DE=SPRINT ) TI ME . OB J , XI , X2 , U 
PRINT OBJ.C,XI ,X2,U; 

CONTROL FINTIM=4. 

METHOD STIFF 
RELERR 03J=l.E-4 
ABSERR OBJsl.E-4 

FIXED I , II , lERI , lERR.NN.NM, IVl .MATZ 
TABLE BB(I-3)=3-0. ,RCI-2)=3«0. ,H(l-2)=3“0. 

PARAM CI=1. ,CF=0. .CSP^^O. , TSP = I . 5 . USP = 1 0 . 



PARAM WI=0. , W2=l. 

PARAM II=0,MATZ=0,NM=3.NN=3,I=I.J=I,IERI=0,IERR=0 



2470 



end 

-- KLENGH = 4245, KPOINT = 1775, AVAILABLE SPACE LEFT IN CDMMON/CURVAL/ = 

-» STIFF INTEGRATION METHOD USED --- 
INITIAL POLES ARE ; (-0 . I OOOE *02 , 0.1000E*02 J) (-0 , I000E*02 ,-0 . 1000E*02 J3 
INITIAL GAINS G( 1 ) ,G(2 ) ,G( 3 ) = -0.3800E*00 0.6000E+00 -0.6000E*0I 

state control matrix (B AB) = 

O.OOOOE+00 0.1000E*05 
0.2000E*03 -0.4000E+05 

OUTPUT CONTROL MATRIX (C8 CAB) = 

O.OOOOE+00 -0. 1000E*05 
OBSERVABILITY MATRIX (C- A-C- ) = 



-0. I000E*01 
0.0000E*00 



O.OOOOE+00 

-0.5000E*02 



DOUBLE-WORDS --- 
3000E+03, O.OOOOE^OO J) 



AAAAA DDDDDD 
A AD D 

A AD D 

AAAAAAA D D 

A AD D 

A AD D 

A A DDDDDD 



ssssss 

s 

s 

sssss 

s 

s 



ssssss 
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FORTRAN PROGRAM 
FOR 

automated design synthesis 

VERSION 1,10 



CONTROL parameters 
ISTRAT » 0 lOPT 


3 


5 


lONED = 


7 


IPRINT s 3050 


IGRAD = 0 NDV 


3 


3 


NCON 


5 





nnk dsl output listing, group 1 



DC MOTCR INTEGRAL CCNTROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 



TIME 


08J 


C 


XI 


X2 


U 


O.OOOOOE^OO 


O.OOOOOE‘00 


1,0000 


O.OOOOOE*O0 


0. OODOOE-^OO 


O.OOOOOE 


1 ,00000E-0C 


1.94922E-02 


0. 99486 


5. 13754E-03 


2.52547E-02 


4 , 67 1 07E* 


2,OOOOOE-02 


3. 78729E-02 


0.97439 


2.56087E-02 


5.58440E-02 


9.52865E- 


3.00000E-02 


5.-;7989E-02 


0.94014 


5.98551E-02 


8 .04192E-02 


0.15103 


4.00000E-02 


7. 00222E-02 


0.89509 


0. 10491 


9.91254E-02 


0,21202 


5.00000E-02 


8.3A120E-02 


0.84193 


0.15807 


0.11263 


0.27627 


6.00000E-02 


9.A9395E- 02 


0.78319 


0.2168 1 


0.12165 


0.34196 


7,OOOO0E-O2 


0, 10A67 


0.72092 


0.27908 


0. 12683 


0,40765 


8.00000E-02 


0.11272 


0.65689 


0.34311 


0. 128 78 


0.47215 


9.00000E-02 


0.11926 


0.59259 


0.40741 


0.12805 


0.53451 


1 . OOOOOE-01 


0. 12^^7 


0.52921 


0.47079 


0-12513 


0.59399 


0.11000 


0.1285A 


0.46775 


0.53225 


0. 12047 


0.65003 


0.12000 


0.13165 


0.40896 


0.59104 


0.11448 


0.70225 


0.13000 


0.13397 


0.35343 


0.64657 


0.10751 


0.75039 


0. lAOOO 


0.13567 


0.30157 


0.69843 


9 .98488E-02 


0.79433 


0.15000 


0.13688 


0.25365 


0.74635 


9.17725E-02 


0.83402 


0. 16000 


0.13771 


0.20983 


0.79017 


8.34994E-02 


0.86952 


0.17000 


0.13827 


0.17015 


0.82985 


7.52150E-02 


0.90094 


0. 18000 


0. 13862 


0.13459 


0.86541 


6. 70709E-02 


0.92847 


0.19000 


0. 1388A 


0.10304 


0.89696 


5.9188 7E-02 


0.95229 


0.20000 


0.13896 


7.53384E-02 


0.92466 


5.16647E-02 


0.97267 


0.21000 


0.13903 


5. 12988E-02 


0.94870 


4.45697E-02 


0.98984 


0.22000 


0.13906 


3,06886E-02 


0.96.931 


3 . 79544E-02 


1.0041 


0.23000 


0.13908 


1.32590E-02 


0. 98 6 74 


3. 18513E-02 


1,0157 


0.2A00O 


0. 13910 


-1.25109E-03 


1.0013 


2.62778E-02 


1,0249 


0.25000 


0.13911 


-1.31081E-02 


1 .0131 


2. 12392E-02 


1.0320 


0.26000 


0.13913 


-2.25784E-02 


1.0226 


1 .67296E-02 


1.0372 


0.27000 


0.13915 


-2.99229E-02 


1.0299 


1.27346E-02 


1.0408 


0.28000 


0. 13918 


-3.53949E-02 


1.0354 


9 .23358E-03 


1.0430 


0.29000 


0,13922 


-3.92343E-02 


1.0392 


6.20011E-03 


1.0440 


0,30000 


0.13925 


-4. 16677E-02 


1,0417 


3.60431E-03 


1.0441 


0,31000 


0.13929 


-4.29059E-02 


1,0429 


1.41384E-03 


1.0433 


0.32000 


0. 13933 


-4.31429E-02 


1.0431 


-4.05144E-04 


1.0419 


0.33000 


0.13936 


-4.25561E-02 


1 .0426 


-1 .88703E-03 


1.0400 


0.3A000 


0. 139^0 


-4 . 13057E-02 


1.0413 


-3.06614E-03 


1.0377 


0.35000 


0. 139A3 


-3.95339E-02 


1.0395 


-3.97597E-03 


1.0352 


0.36000 


0.13945 


-3.73683E-02 


1.0374 


-4 .64895E-03 


1 , 0324 


0.37000 


0. 13947 


-3.49191E-02 


1.0349 


-5. 11580E-03 


1 .0296 


0.38000 


0. 13949 


-3.22017E-02 


1.0323 


-5.40541E-03 


1.0268 


0.39000 


0.13951 


-2.95377E-02 


1.0295 


-5.54461E-03 


1.0240 


O.AOOOO 


0.13952 


-2.67572E-02 


1.0268 


-5.55833E-03 


1.0212 


O.AIOQO 


0.13953 


-2.39965E-02 


1.0240 


-5.46898E-03 


1 .0186 


0.«2000 


0 .13954 


-2. 13020E-02 


1.0213 


-5.29686E-03 


1.0161 


0. A3000 


0.13955 


-1 .87104E-02 


1 .0187 


-5.06008E-03 


1.0138 


0, AAOOO 


0.13955 


-1 .62500E-02 


1,0163 


-4.77465E-03 


1.0116 



. STARTING RUN 1 
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0 .A5000 


0. 13956 


-1 .39415E-02 


1.0139 


'4.45455E-03 


1 .0097 


0.A6000 


0. 13956 


-1.17992E-02 


l.OllB 


-4.111B7E-03 


1.0079 


0,47000 


0.13956 


-9.B3171E-03 


1.0098 


-3. 75689E-03 


1.0063 



SCALAR 


program parameters 














REAL parameters 
















1 ) 


ALAM02 = 


O.OOOOOE+00 


20) 


EXT RAP 




3 


0.50000E*01 




2) 


BETAMC = 


0 . 00000E*00 


21) 


FOCH 




3 


0. lOOOOE-01 




3 ) 


CT 


-0.30000E-01 


22) 


FOCHM 




3 


0. lOOOOE-02 




4) 


CTL - 


-0.50000E-02 


23) 


GMULTZ 




3 


0. 10000E*02 




5) 


CTLMIN = 


0. lOOOOE-02 


24) 


PSAIZ 




3 


0.95000E+00 




6 ) 


CTMIN a 


O.lOOOOE-01 


25 ) 


RMULT 




3 


0.50000E^01 




7) 


DABALP = 


0. lOOOOE-03 


26) 


RMVLMZ 




3 


0.20000E+00 




B) 


DABDBJ = 


0. 13956E-03 


27) 


RP 




3 


0. lOOOOE + 02 




9) 


DABDSM = 


0.27913E-03 


28) 


RPMAX 




3 


0. lOOOOE+11 




10) 


DABSTR = 


0. 13956E-03 


29) 


RMULT 


S 




0.20000E*00 




11 ) 


DEL ALP = 


0.50000E-02 


30) 


RPPMIN 


3 




0. lOOOOE-09 




12) 


DELDBJ = 


O.lOOOOE-02 


31) 


RPPRIM 


3 




0. lOOOOE + 03 




13) 


DELD8M a 


0. lOOOOE-01 


32) 


SCFD 


a 




0. lOOOOE^Ol 




14 ) 


DELSTR 3 


0. lOOOOE-02 


33) 


3CLM1N 


a 




0. 10000E-02 




15) 


DLDBJl 3 


0. lOOOOE-00 


34) 


STDL 


3 




0. lOOOOE-02 




16 ) 


DLDBJ2 3 


0. 10000E*04 


35) 


THETAZ 


3 




0. lOOOOE^OO 




17) 


0X1 


0. lOOOOc-01 


36 ) 


XMULT 


S 




0.261BOE-*-01 




IB) 


DX2 


O.ZOOOOE^OO 


37 ) 


ZRD 


3 




0. lOOOOE-04 




19) 


EPSPEN 3 -0.50000E'01 


3B) 


PMLT 


S 




0. 10000E*02 




INTEGER parameters 














1) 


ICNOIR 3 


4 4 ) 


ITRMDP 


3 3 






6 ) JONEO 3 


7 


% 


ISCAL 3 


1 5) 


ITRMST 


3 2 






7 ) JTMAX 3 


20 


3) 


ITMAX 3 


40 















ARRAY STORAGE REQUIREMENTS 

DIMENSIONED REQUIRED 
ARRAY SIZE SIZE 

WK 1000 299 

IWK 1000 196 



SCALING INFORMATION CALCULATED BY ADS 

OBJECTIVE FUNCTION SCALE FACTOR. SCFD = 0.A66B5E+02 

DESIGN VARIABLE SCALE FACTORS 

1 0.65528E^01 0.168^2E^01 0.37A7BE+00 

CONSTRAINT SCALE FACTORS 

CONSTRAINT NUMBERS 

15^23 
SCALE FACTORS 

0.6B265E+01 O.lOOOOE-02 0.10000E*02 0.101A7E+00 0.73710E^Ol 



IDPT = 5; MODIFIED METHOD OF FEASIBLE DIRECTIONS 



— INITIAL DESIGN 

OBJ s 0,13956E^00 

DECISION VARIABLES (X-VECTOR) 

1) -0.38000E+00 0.60000E+00 -0,60000E*01 

LOWER BOUNDS ON THE DECISION VARIABLES (VLB-VECTDR) 
n -O.IOOOOE^OA -O.lOOOOE+04 -O.lOOOOE+04 
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UPPER 30UN0S ON THE DECISION VARIABLES (VUB-VECTOR) 
n 0.10000E*0A O.IOOOOE^OA 0.10000E*0A 



constraint values CG-VECTOR) 

i; -a.C^CoIE*0l -G.1036B£*00 -0.6601AE*02 -0.10000E*02 -0.10000E*02 



— .EG IN : "'RAT I ON NUH8ER 

“^HERE iPE 0 active CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE 3I0E CONSTRAINTS 

GRADIENJ OF '"HE OBJECTIVE FUNCTION (OF' VECTOR) 

1) Q.67675E*00 O.lCSCBE-^Ol 0.22360E*01 

SEARCH DIRECTION CS-VECTOR) 

n -0.30266E*00 -0.56A74E+00 'O.IOOOOE^Ol 

PROPOSED AL?‘^A = 0.20000E*00 

CalCULATHO alpha = 0.32652£*01 

OBJECTIVE = 0.5709AE-01 

:ec:cxon yARiABLEc cx- vector) 

I) -0.S3031E*GC -G.A9A8BE+00 -0.1A712E*02 

ZCNSTRAINT VALUES (G-VECTOR) 

1) -O.IA373E^OI -0.I2AABE*00 -0.6AA9AE*02 -0.I0000E^02 -0.I0000E*02 



— BEGIN ITERATION NUMBER 2 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 

1) -0.A9351P*00 0.53307E*00 -0.I915AE-01 

SEARCH DIRECTION (S-VECTOR) 

1) 0.70557E*00 -O.IOOOOE^OI -0.23626E+00 

PROPOSED ALPHA * 0.I7326E+0I 

CALCULATED ALPHA = 0.25660E-0I 

OBJECTIVE 5 0.56591E-01 

DECISION VARIABLES (X-VECTOR) 

1) -0.52B05E*00 -0.51012E*00 -0.U729E*02 

CONSTRAINT VALUES (G-VcCTOR) 

1) -0.1A636E*0l -0.12500E*00 -0.6A5I1E*02 -O.lOOOOE+02 -O.IOOOOE^OC 



~ BEGIN ITERATION NUMBER 3 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

I) -0.A8796E*00 0.596A5E+00 -0.22670E-01 

SEARCH DIRECTION (S-VECTOR) 

1) 0.75708E*00 -O.lOOOOE^Ol -0.11070E*00 
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PROPOSED ALPHA a 0.16A55E+01 



CALCULATED ALPHA » 0.69B95E-02 

OBJECTIVE a 0.56A56C-01 

OECiniCN /ARIABLES ^X-V£CTCR) 

:: -C.S272Ac*00 -0.51A27E^00 -0.14731E*02 

constraint values (G-VECTOR) 

-0,1A7A5E*01 -0.12506E + 00 -0.6A519E*02 -O.lOOOOE+02 -O.lOOOOE + 02 



— BEGIN ITERATION NUMBER A 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

n -0.AB500EO0 0.52880E + 00 -0.23835E-01 

SEARCH OIPECTIQN (S-VECTOR) 

n 0.303I2E*00 -0.10000E*01 -0,61035E-01 

PROPOSED ALPHA = 0.16325E-01 

SALCULaTcD alpha a 0.22228E+00 

OBJECTIVE i 0.5386AE-01 

DECISION variables (X-VECTOR) 

1) -0.A9983E+00 -0.6A625E+00 -0.1A767E+02 

CONSTRAINT VALUES (G- VECTOR I 

1) -0.175A1E*01 -0.121B1E*00 -0.6A550E^02 -0.10000E*02 -0.10000E*02 



— BEGIN ITERATION NUMBER 5 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS ‘ 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

n -0.47B29E+00 -0.33528E+00 -0.87296E-01 

SEARCH DIRECTION (S-VECTOR) 

n O.lOOOOE+01 0.70100E+00 0.1B252E+00 

PROPOSED ALPHA a 0 . 1 !<;6<5£ + 00 

CALCULATED ALPHA =_ 0.17350E+00 

OBJECTIVE a 0.52701E-01 

DECISION VARIABLES (X-VECTOR) 

1) -0.A7335E+00 -0 . 57A0<;E*00 -0.1A6B2E+02 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.20A09E*01 -0.12899E*00 -0.65109E*02 -0,10000E+02 -O.lOOOOE+02 



— BEGIN ITERATION NUMBER $ 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 
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GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 

1) -0.2A605E+00 0.38073E+00 0.44972E-01 

SEARCH DIRECTION (S-VECTOR) 

1) 0.10000E*01 -0.34160E*00 0.13126E-01 

PRORCSEO ALPHA = 0.19789£*00 

CALCULATED ALPHA = 0.54227E-»00 



OBJECTIVE 3 0.50417E-01 



DECISION VARIABLES (X-VECTOR) 

1 ) -0.39060E*00 -0'.68402E>00 -0.14663E*02 



CONSTRAINT VALUES (G-VECTOPJ 

n -0,26250E*01 -0.12582E*00 -0.65743E^02 -O.lOOOOE+02 



— BEGIN ITERATION NUMBER 7 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 

1) -0.52639E-01 -0.95469E-01 0.10113E*00 

SEARCH OIRECTION (S-VECTOR) 

1) O.lOOOOE-^01 0.69751E*00 -0.91949E*00 

PROPOSED ALPHA ^ 0.35788E+00 

CALCULATED ALPHA = 0.99532E-01 

OBJECTIVE 3 0.50195E-01 

DECISION VARIABLES (X-VECTOR) 

1) -0.37541E*00 -0.64280E*00 -0.14908E+02 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.26885E*01 -0.13443E*00 -0.65935E+02 -0.10000E*02 



— BEGIN ITERATION NUMBER 8 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRAOIENT OF THE OBJECTIVE FUNCTION COF-VECTOR ) 

1) -0.19434E-02 0.22563E+00 0.14040E*00 

SEARCH OIRECTION (S-VECTOR) 

1) 0.76057E+00 0.38162E-01 -O.lOOOOE^Ol 

PROPOSED ALPHA = 0.32090E*00 

CALCULATED ALPHA = 0.27000E*03 

OBJECTIVE s 0.35264E-01 

OECISION variables (X- VECTOR ) 

1 ) 0.30963E + 02 0.54750E-»01 -0.73533E + 03 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.27515E*01 -0.12889E*00 -0.66352E^02 -O.lOOOOE+02 



-O.lOOOOE+02 



-O.lOOOOE+02 



0. lOOOOE+02 
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BEGIN ITERATION NUMBER 



9 



THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIOE CONSTRAINTS 

GRADIENT DF THE DBJECTIVE FUNCTION (OF-VECTOR) 

1) 0.70317E-02 0.13976E-01 0.59360E-02 

SEARCH DIRECTION (S-VECTOR) 

1) -0.50312E*00 -O.lOOOOE+01 -0.42831E*00 

PROPOSED alpha » 0.13505E*03 

CALCULATED ALPHA = O.OOOOOE^OO 

OBJECTIVE = 0.3526AE-01 

DECISION VARIABLES (X-VECTOR) 

n 0.30963E*02 0.5A750E*01 -0.73533E*03 

constraint VALUES (G- VECTOR) 

1) -0.27515E+01 -0.12889E*00 -0.66352E*02 -0.10000E*02 -0.10000E*02 



— BEGIN ITERATION NUMBER 10 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIOE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

1) 0.70317E-02 0.13976E-01 0.59B60E-02 

SEARCH DIRECTION (S-VECTOR) 

1) -0.50312E*00 -O.lOOOOE+01 -0.A2B31E*00 

PROPOSED ALPHA = 0.13500E+03 

CALCULATED ALPHA = O.OOOOOE^OO 

OBJECTIVE » 0.3526AE-01 

DECISION VARIABLES (X-VECTOR) 

1) 0.30963E+02 0.5A750E+01 -0.73533E*03 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.27515E*0l -0.12889E*00 -0.66352E+02 -O.lOOOOE+02 -O.lOOOOE+02 



~ BEGIN ITERATION NUMBER II 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIOE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

1) 0.70317E-02 0.13976E-01 0.59B60E-02 

SEARCH DIRECTION (S-VECTOR) 

1) -0.50312E+00 -0.10000E*01 -0.A2831E*00 

PROPOSED ALPHA = O.lOOOOE-02 

calculated alpha s 0.5776AE*0l 

OBJECTIVE - 0.33286E-01 
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DECISION VARIABLES (X-VECTOR) 



1) 0,30520E*02 0.20^53E*01 

CONSTRAINT VALUES (G-VECTOR) 

IJ -0-.276005*Cl -0.1296bE*00 



-0.74193E*03 



-0.66357E*02 -O.IOOOOE*02 



-0. 10000E*02 



— ?EGIN ITERATION NUMBER 12 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRAOIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

n 0.72741E-02 0.6B029E-0I 0.56003E-02 

SEARCH DIRECTION (S-VECTORI 

n -0.45297E*00 -O.lOOOOE^Ol -0.3B451E*00 

proposed alpha 5 0.2BB82E*0l 

CALCULATED ALPHA r 0.5u260E*00 

OBJECTIVE = Q.33270E-01 

OECISIGN VARIABLES (X-VECTCR ) 

,) :.3C4B5E-02 a.l74o9E*01 -0.742A5E*03 

IGNSTRAINT VALUES :0-VECT0RI 

n -0.27534E*01 -0.i2B57E*00 -0.662B5E^02 -0.10000E^02 -0.10000E^02 



— BEGIN ITERATION NUMBER 13 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION CDF-VECTOR) 

1) 0.70757E-02 0.44616E-02 0.54342E-02 

SEARCH DIRECTION (S-VECTOR) 

1) -0.10000E*01 -0.63056E+00 -0.76B01E*00 

PROPOSED ALPHA s 0.31395E+0I 

CALCULATED ALPHA = 0.10690E-01 

OBJECTIVE = 0.33101E-OI 

DECISION variables (X-VECTOR) 

1 ) 0.30322E*02 0.13466E + 01 - 0 . 7<5<56<5E* 03 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.2759BE*01 -0.12961E>00 -0.656BBE*02 -O.lOOOOE+02 -O.lOOOOE+02 



— BEGIN ITERATION NUMBER 14 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRAOIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

1) 0.72210E-02 -O.lOlOlE-01 0.59956E-02 

SEARCH OIRECTION (S-VECTOR) 

1) -O.lOOOOE^Ol 0.75990E-01 -0.78970E+00 
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PROPOSED ALPHA s 0.7B5B1E+00 



CALCULATED ALPHA = 0.58^6 !£♦ 02 

OBJECTIVE s 0.2A306E-01 

DEC:c:CN VARIABLES (X-VECTOR) 

\) 3.21AC0E*02 0.398<i3E*01 -0.B67B2E*03 

CONSTRAINT VALUES (G-VECTOR) 

1 ) -0.27575E^0I -O.UOl^E + 00 -0.663^7E*02 -0.10000E*02 -O.lOOOOE+02 

— BEGIN ITERATION NUMBER 15 

"HERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

there are 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 

1) 0.22096E-02 Q.10575E-01 0.1A623E-02 

SEARCH DIRECTION (S-VECTOR) 

n -0.100COE*01 -0.63193E^00 -O.77107E^00 

RRCPOSED alpha = 0.29765E*02 

CALCULATED ALPHA ^ O.OOOOOE*00 

OBJECTIVE = 0.2A306E-C1 

DECISION VARIABLES (X-VECTOR) 

1) 0.21AOOE^02 0.39BA3E*01 -0.B67B2E*03 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.27575E^01 -O.lAOlAE+00 -0.663A7E+02 -O.lOOOOE+02 -0.10000E^02 

— BEGIN ITERATION NUMBER 16 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 

1) 0.22096E-02 0,10575E-01 0.1^623E-02 

SEARCH DIRECTION (S-VECTOR) 

n -0.20B95E+00 -0.10000E*01 -0.13B2BE*00 

PROPOSED ALPHA « 0.29230E*02 

CALCULATED ALPHA = O.OOOOOE*00 

OBJECTIVE = 0.2A306E-01 

DECISION VARIABLES (X-VECTOR) 

1) 0.21AOOE+02 0.398A3E+01 -0,867B2E^03 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.27575E*01 -0.1^01AE*00 -0.663^7E+02 -0.10000E*02 -O.lOOOOE+02 

— BEGIN ITERATION NUMBER 17 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 
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GRAOIENT OF THE OBJECTIVE FUNCTION (OF- VECTOR) 

I) 0.22096E-02 O.I0575E-OI 0.14623E-02 

SEARCH DIRECTION (S-VECTOR) 

I) -0.20895E*00 -O.IOOOOE^OI -0.13828S*00 

PROPOSED alpha = O.lOOOOE-02 

CALCULATED ALPHA = 0.2205AE*01 

OBJECTIVE = 0.239IAE-0I 

OECISION VARIABLES (X-VECTOR ) 

I) 0.2I330E>02 0.26749E+0I -0.86863E*03 



CONSTRAINT VALUES (G-VECTOR) 

I) -0.27529E*OI -0.I390IE*00 -0.6524IE-02 -0 . 1 0000E-^02 -O.IOOOOE*02 



— BEGIN ITERATION NUMBER 18 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATEO CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CCNSTRAINTS 

GRAOIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

n 0.29993E-02 0.48979E-02 O.lHlOE-02 

f 

SEARCH OIRECTION (S-VECTOR) 

I) -0.38877E*00 -0 . I OOOOE+0 I -0.I7775E+O0 

PROPOSEO ALPHA = 0.1I027E^OI 

CALCULATEO ALPHA « 0.95684E-OI 

OBJECTIVE = 0.239I3E-0I 

OECISION VARIABLES (X-VECTOR) 

1) 0.21324E*02 0.26180E+OI -0.86868E^03 

CONSTRAINT VALUES (G-VECTOR) 

I) -0.27632E^01 -0.13983E+00 -0.65I54E>02 -O.lOOOOE+02 -O.IOOOOE*02 



— BEGIN iteration NUMBER 19 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATEO CONSTRAINTS 

THERE ARE 0 ACTIVE SIOE CONSTRAINTS 

GRAOIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 

1) 0.259I5E-02 0,25879E-0l 0.13I33E-02 

SEARCH OIRECTION (S-VECTOR) 

1) -O.I0OI4E*OO -O.lOOOOE+01 -0.50749E-OI 

PROPOSEO ALPHA = 0.11505E-»OI 

CALCULATEO ALPHA = O.OOOOOE^OO 

OBJECTIVE = 0.23913E-01 

DECISION VARIABLES (X-VECTOR) 

1) 0.21324E+02 0.26180E*0l -0.86868E+03 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.27632E*0I -0.13983E*00 -0.65I54E^02 -O.IOOOOE*02 -O.lOOOOE+02 
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BEGIN ITERATION NUMBER 20 



THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (UF-VECTOR) 

1) 0.259I5E-02 0.25B79E-01 O.I3133E-02 

SEARCH DIRECTION (S-VECTOR) 

1) -0,10014E*00 -O.IOOOOE^OI -0.50749E-0I 

PROPOSED ALPHA a 0.47B42E-01 

CALCULATED ALPHA » 0.9772BE'02 

OBJECTIVE = 0.23913E-0I 

DECISION VARIABLES (X-VECTOR) 

1) 0,21324E*02 0.26122E+01 -0.B6B6BE+03 

CONSTRAINT VALUES (G-VECTOR) 

I) -0.27635E>OI -0.139B5E*00 -0.65146E*02 -0.l00O0E*02 -O.IOOOOE>02 



FINAL OPTIMIZATION RESULTS 

NUMBER OF ITERATIONS = 20 

OBJECTIVE = 0.2391’E-0I 

decision VARIABLES (X-VECTOR) 

1) 0.21324E*02 0.26122E*01 -0.B6B6BE+03 

CONSTRAINT VALUES (G-VECTOR) 

I) -0.27635E+0I -0.I39B5E+O0 -0.65I46E+02 -O.IOOOOE^02 -0.10000E*02 
CONSTRAINT TOLERANCE* CT = -0.30000E-01 CTL = -0.50000E-02 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

TERMINATION CRITERIA 

RELATIVE CONVERGENCE CRITERION WAS MET FOR 3 CONSECUTIVE ITERATIONS 
ABSOLUTE CONVERGENCE CRITERION WAS MET FOR 3 CONSECUTIVE ITERATIONS 
MAXIMUM K-T RESIDUAL » O.OOOOOE^OO IS LESS THAN O.IOOOOE-02 



OPTIMIZATION RESULTS 



OBJECTIVE FUNCTION VALUE 0.239I3E-0I 



DESIGN VARIABLES 



VARIABLE 

1 

2 
3 



BOUND 

-0.10000E*04 
-0, I0000E*'04 
-O.lOOOOE-04 



LOWER 

VALUE 

0.21324E*02 

0.26122E*01 

-0.B6B6BE*03 



UPPER 

BOUND 

0 .IOOOOE*04 
0. I0O00E*04 
O.IOOOOE+04 
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DESIGN constraints 



1) -0.<i0<i8E>00 -0.1378E*0I -0.3838E>01 -0.1000E*0I -0.1000E*0l 



'UNCTION EVALUATIONS s 173 





<.SL 7UTPUT LISTING. 


CROUP : 








TIME 


08J 


C 


XI 


X2 


U 


, starting run 175 


O.OOOOOE‘00 


O.OOOOOE^OO 


1 . 0000 


0.00000£>00 


0.00000E*00 


O.OOOOOE+00 




i .OOOOOE-02 


1 . 19744E-02 


0.75301 


0.24099 


0.69497 


0.821 II 




2.00000E-02 


I . 91402E-02 


0147626 


0.52374 


0.42929 


0.85409 




3.000C0E-02 


2. 19849E-02 


0.30355 


0.69665 


0.27368 


0.90859 




4.000C0E-02 


2.31383E-02 


0.19342 


0.80658 


0. 17440 


0.94170 




5.00000E-C2 


2.3603IE-02 


0.12329 


0.87671 


0.11115 


0.96280 




6.00000E-C2 


2.37926E-02 


7.85797E-02 


0.92142 


7.08448E-02 


0.97630 




7.00000E-02 


2.3S697E-02 


5.00829E-02 


0.94992 


4.51536E-02 


0.98489 




8 .OOOOOE-02 


2.39009E-02 


3. 19248E-02 


0.96808 


2.8 7808E-02 


0.99037 




9. OOOOCE-02 


2.39136E-02 


2.03445E-02 


0,97966 


1 .83416E-02 


0.99386 




1 .COCOOE-01 


2.391882-02 


1 .29b20E-02 


0.98704 


1 . 16859E-02 


0.99609 




0. i 1000 


2.3Q209E-02 


3.25829E-03 


0.99174 


7.44523E-03 


0.99751 





.-INAL SAINS GC: ) .G(S) .0(3) ^ O.SI32E*OS 0.2612E*01 -0.3687E^03 

final FOLEC are : l -0 . 338 7E*03 . 0.2793E*03 J) (-0 . 338 7E*03 ,-0 . 2793E+03 J) (-0 . 4508E+02 . 0.OOOOE*OO J) 



SUMMARY: ENDING RUN NO. = 175, .>^0. OF .RERUNS REQUESTED = 175 

RUN ''ERMINATED 3Y ’ EN0J08 * 



2. RESULTS FOR STRATEGY 057 WITHOUT SCALING APPLIED 



DSL/VS RELEASE 1 MODIFICATION LEVEL 1 1985 



>>>> DSL SIMULATION INPUT DATA <<<< 

TITLE DC MOTOR INTEGRAL CONT.ROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 
FIXED NRA.NGT.NCOLA.NRWK.NRIWK. 1GRAD,NDV,NC0N , IDG , IC. IWK 
FIXED ISTRAT..IOPT, lONED.IPRlNT.INFO. IPVT 

PAR AM NRA = A0 ,NCOLA=AO » NRWK= 1 000 ,NR IWK= 1000, NDV = 3 ,NCON = 5 , IGRAD^O , INFO=-2 
TA8LE lDG(l-5) = 3«0,2*»2 
TA8LE VL8(l-6) = 3*»-1000. .3»0. 

TA8LE VU8(1-6)=3*»1000. ,3“0. 

PARAM 1STRAT=0, IOPT=5, 10NED=7, IPRINT=1050 
SAVE 08J,X1.X2,U 

GRAPH (DE=3PR1NT ) TIME.08J ,X1 , X2 , U 
PRINT 08J,C.XI ,X2,U 
CONTROL FINTIM=4. 

METHOD STIFF 
RELERR 08J=l.E-4 
A8SERR 08J=l.E-4 

FIXED I ,II,IER1 ,IERR.NN,NM,IV1 ,MAT2 
TA8LE 88(l-3)=3t0. ,R(I-*2)=3»»0. ,H(l-2) = 3**0. 

PARAM CI=l. ,CF=0. ,CSP=40 . ,TSP-1 .5,USP=10. 

PARAM W1 = 0. , W2=?I . 

PARAM II=0,MAT2=0,NM=3,NN=3, 1=1 ,J=1 ,IER1=0, IERR=0 
ENO 

KLENGH = 4245, KPOINT = 1775, AVAILABLE SPACE LEFT IN COMMON/CURVAL/ = 2470 DOUBLE-WORDS 

•••«*• STIFF INTEGRATION METHOD USED *•«•*» 

INITIAL POLES ARE ; (-0 . 1000E^02 , O.iOOOE+02 J) (-0 . 1 00 OE* 02 ,-0 . 1 OOOE+02 J) (-0 . 3000E+03 , 0.0000E*00 J) 
INITIAL GAINS G( 1 ) .G(2 ) ,G(3 ) = -0.3800E*00 0.6000E + 00 -0.6000E*0I 

STATE CONTROL MATRIX (8 AB) = 

O.OOOOE-^00 O.lOOOE + 05 
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0.2000E^03 -0.4000E+05 

OUTPUT CONTROL MATRIX (C8 CA8) = 
O.OOOOE^OO -O.IOOOE^OB 
observability matrix (C-« A#C>*) = 

-O.lOOOE-01 O.OOOOE+00 

O.OOOOE*CO -0.5000E*02 



AAAAA 
A A 

A A 

AAAAAAA 
A A 

A A 

A A 



000000 
0 0 
0 0 
0 0 
0 0 
0 0 
000000 



ssssss 

s 

s 

sssss 

s 

s 

ssssss 



FORTRAN PROGRAM 
FOR 

AUTOMATEO OESIGN SYNTHESIS 

VERSION 1.10 



CONTROL PARAMETERS 

ISTRAT = 0 lOPT = 5 lONEO = 7 IPRINT = 1050 

IGRAO = 0 NOV = 3 NCON * 5 

SCALING NOT USED 

OSL OUTPUT LiSTINGt GROUP 1 



OC MOTOR INTEGRAL CONTROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 



TIME 


D8J 


C 


XI 


X2 


U 


O.OOOOOE+00 


O.OOOOOE^OO 


1.0000 


O.OOOOOE^OO 


O.OOOOOE+00 


0 .OOOOOE + 00 


l.OOOOOE-02 


1 .9A922E-02 


0.99486 


5. 13754E-03 


2.52547E-02 


4.67107E-02 


2.00000E-02 


3.78729E-02 


0.97439 


2.56087E-02 


5.58440E-02 


9.52865E-02 


3.00000E-02 


5.A7989E-02 


0.94014 


5.98551E-02 


8.04192E-02 


0.15103 


A .OOOOOE-02 


7.00222E-02 


0.89509 


0.10491 


9.91254E-02 


0.21202 


5.00000E-02 


8.3A120E-02 


0.84193 


0.15807 


0.11263 


0.27627 


6.00000E-02 


9.A9395E-02 


0.78319 


0.21681 


0.12165 


0.34196 


7.00000E-02 


0. 10A67 


0.72092 


0.27908 


0.12683 


0.40765 


8.00000E-02 


0.11272 


0.65689 


0.34311 


0.12878 


0.47215 


9.00000E-02 


0.11926 


0.59259 


0.40741 


0.12805 


0.53451 


1 .OOOOOE-01 


0. 12AA7 


0.52921 


0.47079 


0.12513 


0.59399 


0.11000 


0.12854 


0.46775 


0.53225 


0. 12047 


0.65003 


0. 12000 


0.13165 


0.40896 


0.59104 


0.11448 


0.70225 


0.13000 


0.13397 


0.35343 


0.64657 


0.10751 


0.75039 


0. lAOOO 


0.13567 


0.30157 


0.69843 


9.98488E-02 


0.79433 


0. 15000 


0.13688 


0.25365 


0.74635 


9.17725E-02 


0.83402 


0.16000 


0.13771 


0.20983 


0.79017 


8.34994E-02 


0.86952 


0.17000 


0.13827 


0.17015 


0.82985 


7.52150E-02 


0.90094 


0.18000 


0.13862 


0.13459 


0.86541 


6.70709E-02 


0.92847 


0. 19000 


0.13384 


0.10304 


0.89696 


5.91887E-02 


0.95229 



» STARTING RUN 1 >»«*• 
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0.20000 


0.13896 


7 .5338^E-02 


0.92<566 


5 . 166A7E-02 


0.97267 


0.21000 


0.13903 


5.12988E-02 


0.9<5870 


<5.<556 9 7E-02 


0.9898<5 


0.22000 


0.13906 


3.06886E-02 


0.96931 


3.795<5<5E-02 


1 .00<51 


0.23000 


0.13908 


1 .32590E-02 


0.9867<5 


3.13513E-02 


1.0157 


0.2<i000 


0.13910 


-1 .25109E-03 


1.0013 


2.62778E-02 


1.02<59 


0.25000 


0.13911 


-1.31081E-02 


1.0131 


2. 12392E-02 


1.0320 


0.26000 


0.13913 


-2.2578<5E-02 


1.0226 


1 .67296E-02 


1.0372 


0.27000 


0.13915 


-2.99229E-02 


1.0299 


1 .273<56E-02 


1 .0<508 


0.28000 


0.13918 


-3.53949E-02 


1 .035<5 


9.23358E-03 


l.0<530 


0.29000 


0.13922 


-3.92343E-02 


1.0392 


6.20011E-03 


1 .0<5<50 


0.30000 


0.13925 


-4.I6677E-02 


1 .0<517 


3.60<531E-03 


1.0<5<5! 


0.31000 


0.13929 


-<5.2905 9E-02 


1 .0<529 


l.<5 138<5E-03 


1 .0<533 


0.32000 


0.13933 


-<5.31<529E-02 


1 .0<531 


-<5.051<5<5E-0<5 


1.0<5 19 


0.33000 


0.13936 


-<5.25561E-02 


1 .0<526 


-1 .88703E-03 


1 .0<500 


0.3A000 


0. 139<50 


-<5.I3057E-02 


1.0<513 


-3.0661<5E-03 


1.0377 


0.35000 


0. 139<53 


-3.95339E-02 


1.0395 


-3.97597E-03 


1.0352 


0.36000 


0. 139<55 


-3.73683E-02 


1 .037<5 


-<5 .6<5895E-03 


1 .032<5 


0.37000 


0. 139<57 


-3.<59191E-02 


1.03<59 


-5. 1158 0E-03 


1.0296 


0.38000 


0.139<59 


-3.22817E-02 


1.0323 


-5.<505<51E-03 


1.0268 


0.39000 


0.13951 


-2.95377E-02 


1.0295 


-5.5<5<561E-03 


1.02A0 


0. AOOOO 


0.13952 


-2.67572E-02 


1.0268 


-5.55833E«03 


1.0212 


0.^1000 


0.13953 


-2.39965E-02 


1 .02<50 


«5.<56898E«03 


1.0186 


0.^2000 


0. 1395A 


-2. 13020E~02 


1.0213 


-5.29686E-03 


1.0161 


0. <53000 


0.13955 


«1 .8710<5E-02 


1.0187 


«5.06C03E-03 


1.0138 


0. <5 <5000 


0.13955 


-1 .62500E-02 


1.0163 


-<5 . 77<565E-03 


1.0116 


0.A50O0 


0.13956 


-1.39<515E-02 


1.0139 


-<5. <55<555E<03 


1.0097 


0. <56000 


0.13956 


-1.17992E-02 


1.0118 


-<5.11I87E-03 


1.0079 


0. <57000 


0.13956 


-9.83171E-03 


1.0098 


-3. 75689E-03 


1.0063 



lOPT s 5; MODIFIED METHOD OF FEAS18LE DIRECTIONS 



— INITIAL DESIGN 
08J = 0.13956E*00 



DECISION VARIA8LES (X-VECTOR ) 

n -0.38000E*00 0.60000E+00 -0.60000E>01 

LOWER 80UNDS ON THE DECISION VAR1A8LES (VL8-VECT0R) 
n -O.IOOOOE+OA -0.10000E*04 -O.IOOOOE^OA 

UPPER 80UNDS ON THE DECISION VAR1A8LES (VU8-VECT0R) 

1) 0.10000E>OA O.IOOOOE+OA O.IOOOOE-OA 

CONSTRAINT VALUES (G-VECTOR) 

n -0.35685E+00 -0.10217E+01 -0.89558E+01 -O.lOOOOE^Ol -0,lOOOOE*Ol 



— 8EGIN ITERATION NUM8ER 1 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE 08JECTIVE FUNCTION (DF-VECTOR I 

n 0.9A990E-01 0.<5556E-01 0.17951E-01 

SEARCH DIRECTION (S-VECTOR) 

n -O.IOOOOE + OI -0.<;7959E^00 -0.18897E*00 

PROPOSED ALPHA = 0.11608E+00 

CALCULATED ALPHA = 0.22788E*00 

08JECT1VE = 0.12996E+00 
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DECISION VARIABLES (X-VECTOR) 



1) -0.607B8E*00 0.^9071E*00 -0.60A31E*01 



CONSTRAINT VALUES (G-VECTOR) 

n -0.19A85E*00 -0.72B29E*00 -0.87893E*01 -0.10000E*01 -O.lOOOOE+01 



— BEGIN ITERATION NUM8ER 2 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRAOIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 

1 ) -0.59658E-01 0.813A<;E-01 0.117A9E-01 

SEARCH DIRECTION (S-VECTOR) 

1) -0.21330E+00 -O.lOOOOE^Ol -0.22827E+00 

PROPOSED ALPHA = 0.I7198E+00 

CALCULATED ALPHA = 0.900^9E*00 

OBJECTIVE = 0.10507E*00 

DECISION VARIABLES (X-VECTOR) 

1 ) -0.79995E*00 -0 . <5 097 7E-00 -0.62^86E*01 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.69A59E-02 -0.78792E*00 -0.8580^E*01 -0.10000E*01 -O.IOOOOE^Ol 



-- BEGIN ITERATION NUM8ER 3 

THERE ARE 1 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 
CONSTRAINT NUMBERS 
1 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCT ION *( OF- VECTOR ) 

1) -0.A3352E^00 0.78925E-01 -0.22729E-02 

GRADIENT OF CONSTRAINT NUMBER 1 

1) -0.2A212E*01 0.16891E+00 -0.53977E-01 

K-T PARAMETERS. BETA = 0.50816E+00 MAX. RESIDUAL =* O.A91B^E + 00 

SEARCH DIRECTION (S-VECTOR) 

1) 0.10000E*01 -O.lSaoeE+OO 0.52^29E-02 



PROPOSED ALPHA = 0.56^18E+00 

CALCULATED ALPHA * 0.15068E*00 



OBJECTIVE = 0.80689E-01 

DECISION VARIABLES (X-VECTOR) 

1) -0.6A927E+00 -0 . 372 1E*00 -0.62A7BE*01 



CONSTRAINT VALUES (G-VECTOR) 

1 ) -0.32129E+00 -0.U902E+01 -0.B9122E*01 -O.lOOOOE^Ol -O.lOOOOE + 01 



— BEGIN ITERATION NUMBER 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATED CONSTRAINTS 



121 



THERE ARE 



0 ACTIVE SIDE CONSTRAINTS 



GRADIENT OF THE OBJECTIVE FUNCTION (DF- VECTOR) 

I) 0.10485E-01 0.5086CE-01 0.67057E-02 



SEARCH DIRECTION (S-VECTCR) 

i) -O.C:518E*03 -n. lOaoOE'^01 -0.12685£*00 



-PROPOSED ^LPHA ^ :.SZ559E*00 



CALCULATED ALPHA s 0.5CC68E-02 

OBJECTIVE - 0.30397E-01 

CECISION VARIABLES IX-VECTOR) 

I) -0.650S0E + 00 -0.<5^r^3E*00 -0.62485E + 01 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.32II6E*00 -0,1I920E*01 -0.891I9E-0I -O.lOOOOE^Ol -O.lOOOOE+01 



-- BEGIN ITERATION NUMBER 5 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

^i^ERE ARE a ACTIVE SIDE CONSTRAINTS 

DRACIENT OF ’HE OBJECTIVE FUNCTION CDF- VECTOR) 

I) 0.1 1906E-01 T.52925E-OI 0.66<;78E-02 

SEARCH DIRECTION (S-VECTOR) 

1) -0.23055E+00 -O.lOOOOE+01 -0.I2623E+00 

PROPOSED ALPHA a 0.77955E-01 

CALCULATED ALPHA ^ 0,20<509E+00 

OBJECTIVE a 0.70223E-01 

DECISION VARIABLES (X- VECTOR) 

1) -0.69756E+00 -0.6^652E*00 -0.627^3E*01 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.30233E*00 -0.I2122E-01 -0.88751E+01 -O.IOOOOE^Ol -O.IOOOOE+OI 

— BEGIN ITERATION NUMBER 6 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (DF- VECTOR) 

1) -0.31I21E-0I O.A<775E-OI 0.35<599E-02 

SEARCH DIRECTION (S-VECTOR) 

1) 0.<550<57E-0I -O.IOOOOE^OI — 0.II225E^00 

PROPOSED ALPHA = 0.I0466E*00 

CALCULATED ALPHA = 0.10A66E*00 

OBJECTIVE s 0.67007E-0I 

DECISION VARIABLES (X-VECTOR) 

1 ) -0.6928AE-^00 -0.751I8E+00 -0.62860E + 0I 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.32042E+00 -0.11341E+01 -0.88657E+01 -O.IOOOOE+01 -O.lOOOOE+01 
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BEGIN ITERATION NUMBER 



7 



there are 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

there 4R£ 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT CF “HE OBJECTIVE FUNCTION COF-VECTOR) 

— 0.57436E-0C 0.18823E-02 

SEARCr: DIRECTION (S-VECTOR) 

n 0.A5896E*00 -0.10000E*01 -0.12388E*00 

PROPOSED alpha = 0.15437E*00 

CALCULATED ALPHA a 0.20074E-01 

OBJECTIVE = 0.66636E-01 

DECISION VARIABLES (X-VECTOR) 

n -0.68363E*00 -0.77125E*00 -0.62885E+01 

CONSTRAINT VALUES (G-VECTOR) 

[) -0.3A051E*00 -0.11226E*01 -0.88804E*01 -O.lOOOOE^Ol -O.lOOOOE^Ol 



— BEGIN ITERATION NUMBER 8 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION CDF-VECTOR) 

1) -0.36555E-01 -0.67581E-02 0.19196E-02 

SEARCH DIRECTION (S-VECTOR) 

1) O.lOOOOE^Ol 0.18488E+00 -0.52513E-01 

PROPOSED ALPHA a 0.62366E-01 

CALCULATED ALPHA = 0.241 17E-01 

OBJECTIVE = 0.66174E-01 

DECISION VARIABLES (X-VECTOR) 

1 ) -0.65951E + 00 -0.76680E + 00 -0 . 62898E-^0 1 

CONSTRAINT VALUES (G-VECTOR ) 

1) -0.37897E*00 -0.11891E>01 -0.89462E*01 -0.10000E*01 -O.lOOOOE+01 



— BEGIN ITERATION NUMBER 9 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 

1) 0.82166E-02 0.94493E-02 0.38751E-02 

SEARCH DIRECTION (S-VECTOR) 

1) -0.42778E^00 -O.IOOOOE^Ol -0.47765E>00 

PROPOSED ALPHA = 0.22095E-01 

CALCULATED ALPHA = 0.67373E-02 

OBJECTIVE » 0.66144E-01 
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DECISION VARIABLES (X-VECTOR) 



1) -0.66239E+00 -0.77353E*00 -0.62930E+01 



constraint values (G-VECTOR) 

1) -0.37A21E+00 -0.11589E*01 -0.89355E+01 -0.10000E*01 -O.lOOOOE+01 



— 8EG1N ITERATION NUMBER 10 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (DF^^ VECTOR) 

n 0.78501E-03 0.27K1E^02 0.3<5 708E-02 

SEARCH DIRECTION (S-VECTOR) 

1) -0.22617E*00 -0.78196E+00 -O.IOOOOE+Ol 

PROPOSED ALPHA = 0.15A27E-01 

CALCULATED ALPHA 0.21057E-02 

OBJECTIVE = 0.66U2E-01 

DECISION VARIABLES (X-VECTOR) 

1) -0.66287E*00 -0.77518E+00 -0.62951E*01 

CONSTRAINT VALUES (G-VECTOR) 

1 ) -0.3730AE+00 -0.H599E + 01 -o’89333E*Ol -O.IOOOOE+Ol -O.lOOOOE + 01 



— BEGIN ITERATION NUMBER 11 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (DF-VECTOR) 

1) -0.86A7AE-03 -0.80A95E-01 0.33776E-02 

SEARCH DIRECTION (S-VECTOR) 

1) 0.107A3E-01 O.lOOOOE+01 -0.A1960E-01 

PROPOSED ALPHA » O.AA215E-02 

CALCULATED ALPHA = 0.11576E-01 

OBJECTIVE = -0.65590E-01 

DECISION VARIABLES (X-VECTOR) 

1) -O.66275E^00 -0.76360E+00 -0.62956E+01 

CONSTRAINT VALUES (G-VECTOR) 

1) -0.39095E+00 -0.13935E+01 -0.89388E^01 -O.lOOOOE^Ol -O.IOOOOE+Ol 



— BEGIN ITERATION NUMBER 12 

THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION CDF-VECTOR) 

1 ) 0.9207^iE-01 0.21756E-01 0.12965E-01 

SEARCH DIRECTION (S-VECTOR) 

1) -0.99A96E+00 O.lOOOOE+01 -0.19393E*00 
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PROPOSeO ALPHA a 0.68607E-02 

CALCULATED ALPHA a O.OOOOOE+00 

OBJECTIVE a 0.65590E-01 

DECISION VARIABLES (X-VECTOR) 

1) -0.66275E*00 -0.76360E*00 -0.62956E*01 

constraint values (G-VECTOR) 

1 ) -0.39095E*00 - 0 . 1 3 9 35E* 0 1 -0.B93BBE*01 -0.10000E*0I -O.lOOOOE + 01 



— BEGIN ITERATION NUMBER 13 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATEO CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

1) 0.9207<5E-01 0.21756E-01 0.12965E-01 

SEARCH DIRECTION (S-VECTCR) 

n -O.lCOOOE^Ol -0.23629E+00 -0.IA0B1E*00 

PROPOSEO alpha = 0.5787BE-02 

CALCULATEO ALPHA a 0.10B67E-03 

OBJECTIVE = 0.655B9E-0I ' 

DECISION VARIABLES (X-VECTQR) 

1) -0.662B5E*00 -0.76363E*00 -0.62956E+01 

CONSTRAINT VALUES XG- VECTOR! 

I) -0.39071E^00 -0.13935E+01 -0.B9385E*01 -O.lOOOOE+01 -O.lOOOOE+01 



— BEGIN ITERATION NUMBER 

THERE ARE 0 ACTIVE CONSTRAINTS ANO 0 VIOLATEO CONSTRAINTS 

THERE ARE 0 ACTIVE SlOE CONSTRAINTS 

GRADIENT OF THE OBJECTIVE FUNCTION (OF-VECTOR) 

n 0.92007E-01 0,21717E-01 0.12970E-0I 

SEARCH DIRECTION (S-VECTOR) 

1) -O.IOOOOE^OI -0,23616E*00 -0.1A089E*00 

PROPOSEO ALPHA = O.lOOOOE-02 

CALCULATED ALPHA = 0.72687E-04 

OBJECTIVE a 0.65589E-01 

DECISION VARIABLES (X-VECTOR) 

1) -0.66293E*00 -0.76365E+00 -0,62956E*01 

CONSTRAINT VALUES (G-VECTOR ) 

1 ) -0,3905<;E*00 -0.13935E + 01 -0.89383E*01 -0.10000E*0I -O.lOOOOE+01 



FINAL OPTIMIZATION RESULTS 
NUMBER OF ITERATIONS * K 
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OBJECTIVE a 0.65589E-01 



DECISION VARIABLES (X-VECTCR] 

1) -0.66293S-00 -0.76365£*00 -C.62956E*01 

constraint .'AL'JES 'C- ‘/ECTOR) 

♦‘J.C?C54c^0Q -0. 13935E-01 -0.89383E*0l -O.iOOOOE^Ol -0.10000E*01 

CCNSTRAINT tolerance. CT a -0.:OOOOE-Ol CTL a -0.50000E-02 
THERE ARE 0 ACTIVE CONSTRAINTS AND 0 VIOLATED CONSTRAINTS 

THERE ARE 0 ACTIVE SIDE CONSTRAINTS 

TERMINATION CRITERIA 

RELATIVE CONVERGENCE CRITERION WAS MET FOR 3 CONSECUTIVE ITERATIONS 
ABSOLUTE CONVERGENCE CRITERION WAS MET FOR 3 CONSECUTIVE ITERATIONS 



OPTIMIZATION RESULTS 



:bject:ve 


FUNCTION VALUE 


0.bS539E-01 




DESIGN variables 










LOWER 


UPPER 


variable 


BOUND 


VALUE 


BOUND 


1 


-O.IOOOOE+OA 


-0.66293E*00 


O.IOOOOE^OA 




-0. lOOOOE^OA 


-0.76365E*00 


0. lOOOOE^OA 


3 


-0. lOOOOE^OA 


-0.62956E>01 


0. lOOOOE-^OA 



DESIGN CONSTRAINTS 

n -0.3905E*00 -0.139AE*01 -0.893BE-01 -0.1000E*01 -O.lOOOE^Ol 



FUNCTION EVALUATIONS = 90 

HMM DSL OUTPUT LISTING. GROUP 1 



TIME 


OBJ 


C 


XI 


X2 


U 


O.OOOOOE-OO 


O.OOOOOE^OO 


1.0000 


0. 0OOOOE*OO 


O.OOOOOE-^OO 


O.OOOOOE-^OO 


1.00000E-02 


1 .90A51E-02 


0.99076 


9.23799E-03 


5.24768E-02 


0.10900 


2.00000E-02 


3.50A69E-02 


0.9370A 


6.29613E-02 


0.16774 


0.29362 


3.00000E-02 


A.70222E-02 


0.B2258 


0.17742 


0.28635 


0.51578 


A.OOOOOE-02 


5.51022E-02 


0.65749 


0.34251 


0.36476 


0.73190 


5.00000E-02 


6.007775-02 


0.46817 


0.53183 


0.3B167 


0.905B0 


6.00000E-02 


6.28622E-02 


0.28607 


0.71393 


0.33750 


1.0164 


7.00000E-02 


6.A2180E-02 


0.13799 


0.B6201 


0.24934 


1.0604 


B.OOOOOE-02 


6.A726AE-02 


3.96980E-02 


0.96030 


0.14261 


1.0494 


9.00000E-02 


6.A8304E-02 


-5.97585E-03 


1.0060 


4.28098E-02 


1.0042 


1 .OOOOOE-Ol 


6.ABA37E-02 


-7.61740E-03 


1.0076 


-3.09017E-02 


0.94839 


0.11000 


6.«^9337E-02 


1.89452E-C2 


0.9B105 


-6.92305E-02 


0.90176 


0.12000 


6.51313E-02 


5.50209E-C2 


0.94418 


-7.29235E-02 


0.87684 


0.13000 


6.53820E-02 


B.75761E-02 


0.91242 


-5.06554E-02 


0.87737 


FINAL GAINS 


G(1).G(2).G(3) a -0.66 


29E*00 -0, 


.7636E*00 


0.6296E*01 



FINAL POLES ARE i (-0 . 12«5E* 02 . 0-5156E>02 J) (-0 . 12A5E*02 .-0 . 5 15 6E* 02 J) (-0 
SUMMARY: ENDING RUN NO . = 92. NO. OF RERUNS REQUESTED = 92 

RUN terminated by ’ENDJ08’ 



. STARTING RUN 92 



2237E*02. O.OOOOE^OO J) 
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APPENDIX G 

ADSLPC IMSL PROGRAM CHANGES 



Changes were done to the following IMSL routines to make them compatible 
with MS-FORTR.A.N as follows: 



Routine Line Number 


Change 


DGEAR 4260 


SEPS = 1.525878906E-05 


5980 


6HDGEAR to DGEAR ' 


5990 


6MDGEAR to 'DGEAR ' 



LEQTIB 2190 


6HLEQT1B to 'LEQTIB' 



UERTST 530 


define CHARACTER-^6 NAME 


550 


define CHARACTER lEQ 


560 


define CHARACTERS NAMSET, NAMUPK 


570 


NAMSET = 'UERSET' 


580 


NAMEQ = ' 


590 


lEQ = '=' 


620 


NAMUPK = NAME 


640 


lOUNIT = 6 


880 


deleted 


• 890 


IF (NAMUPK.NE.NAMSET) GOTO 25 


1050 


6A1 to 1A6 (TWICE) 


1070 


6A1 to 1A6 (TWICE) 


1090 


6A1 to 1A6 (TWICE) 


1110 


6A1 to 1A6 (TWICE) 


1170 


deleted 


1180 


NAMEQ = NAMUPK 
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EIGRF 1030 
2740 
2770 



RDELP = 1.525878906E-05 
6HEIGRF to 'EIGRF ' 
6HEIGRF to 'EIGRF ' 



EQRH3F400 

3700 



RDELP = 1.52878906E-05 
6HEQRH3F' to 'EQRH3F 



Numerical conversions were calculated as follows: 

• Original IBM data statement: DATA SEPS /Z3C 100000/ 

• Find exponent: 3C - 40 = -3 

• The resulting hexadecimal number is: 0.1 x 10'^ 

• Converting 

• Converting 



to decimal the number is: 16.0 x 10"^ 
to exponential form the value is: 1.525878906E-05 
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APPENDIX H 

ADSLPC PROGRAM LISTING 



L ADSLPC MS-FOTRAN77 PROGRAM 

$nof loatcalls 
$large 
^storage: 2 

C DC MOTOR INTEGRAL CONTROL SYSTEM-NO NOISE-KUO EXAMPLE 8-7 

c 

C LISTING OF ALL VARIABLES 
C 



C AA 


MATRIX OF STATE COEFFICIENTS 


C BB 


VECTOR CF CONTROL COEFFICIENTS 


C C 


SCALAR OUTPUT VALUE 


C CC 


MATRIX OF STATE OUTPUT COEFFICEINTS 


C CF 


FINAL EXPECTED OUTPUT VALUE 


C Cl 


INITIAL OUTPUT VALUE 


C CON 


VECTOR OF CONSTRAINT VALUES 


C CSP 


MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 


C DELC 


DIFFERENCE OF THE ACTUAL OUTPUT TO THE INITIAL OUTPUT 


C DELT 


INITIAL INTEGRATION STEP SIZE 


C DELU 


DIFFERENCE OF THE ACTUAL CONTROL TO THE INITIAL CONTROL 


C DELX2 


DIFFERENCE OF THE ACTUAL X2 STATE TO THE INITIAL X2 STATE 


C DV 


VECTOR OF DESIGN VARIABLE GREADIENTS 


C F 


VECTOR OF INPUT COEFFICIENTS 


C FVl 


EIGENVALUE WORK VECTOR 


C G 


VECTOR OF DESIGN VARIABLE GAIN VALUES 


C H 


VECTOR OF OUTPUT COEFFICEINTS FOR THE INPUTS 


C I 


LOOP COUNTER 


C II 


EVALUATION COUNTER 


C IC 


VECTOR OF CONSTRAINT GRADIENTS 


C IDG 


VECTOR OF CONSTRAINT TYPES 
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C lERR EIGENVALUE ERROR INDICATOR 



C IGRAD 


SCALAR FOR INDICATING METHOD FOR OBTAINING GRADIENTS 


C INFO 


SCALAR INDICATING OPTIMIZATION STATUS 


C lONED 


SCALAR INDICATING TYPE OF CNE-DIMESIONAL SEACH 


C lOPT 


SCALAR INDICATING ’"YPE OF OPTIMIZATION 



C IPRINT SCALAR INDICATING TYPE OF PRINT CONTROL FOR ADS 
C IFVT LINPAC INTEGER WORK VECTOR 
C ISTRAT SCALAR INDICATING TYPE OF STRATEGY 



C IVl 


EIGENVALUE INTEGER WORK VECTOR 


C INK 


INTEGER WORK VECTOR FOR ADS 


C J , 


LOOP COUNTER 


C JOB 


STATUS INDICATOR OF OPTIMIZATION PROBLEM 



C LSLCPE SLOPE OF THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION 
C LST08J value of THE OBJECTIVE FUNCTION FOR THE LAST INTEGRATION 
C LSTTIM VALUE OF TINE AT THE PREVIOUS INTEGRATION STEP 



C MATZ 


EIGENVALUE SCALAR VALUE TO INDICATOR EIGENVALUE DETERMINATION 


C MEAN 


DSL NORMAL DISTRIBUTION MEAN 


C NCOLA 


SCALAR COLUMN DIMENSION OF WA MATRIX 


C NCON 


NUMBER OF CONSTRAINTS 


C NDV 


NUMBER OF DESIGN VARIABLES 


C NGT 


SCALAR ADS CONTROL VARIABLE 


C NM 


DIMENSION OF MATRIX SIZE FOR EIGENVALUE DETERMIANTION 


C NN 


DIMENSION OF MATRIX SIZE FOR EIGENVALUE DETERMIANTION 


C NRA 


SCALAR ROW DIMENSION OF WA MATRIX 


C NRIWK 


SCALAR DIMENSION OF IWK VECTOR 


C NRWK 


SCALAR DIMENSION OF WK VECTOR 


C OBJ 


OBJECTIVE FUNCTION VALUE 


C P 


VECTOR OF DESIRED INITIAL COMPLEX POLE VALUES 


C PA 


MATRIX OF INITIAL GAIN DETERMINATION COEFFICIENTS 


C Q 


STATE WEIGHTING MATRIX 


C R 


CONTROL WEIGHTING VECTOR 


C SD 


STANDARD DEVIATION OF DSL NORMAL FUNCTION 


C SEED 


RANDOM MUMBER SEED FOR DSL NORMAL FUNCTION 


C STEP 


EVALUATION INTERVAL FOR CONSTRAINTS 


C TAU 


SYSTEM TIME CONSTANT 
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C TC 


TIME WHEN COMPLETION TIME SPECIFICATION IS MET 


C TIMEND 


TIMER VALUE AT EACH EVALUATION INTERVAL 


C TOL 


initial' TOLERANCE FOR DGEAR 


C IS? 


MAXIMUM COMPLETION TIME 


c u 


SCALAR FEEDBACK CONTROL VALUE 


C UI 


INITIAL VALUE OF THE CONTROL VECTOR 


C UMAX 


MAXIMUM CONTROL VALUE 


C USP 


MAXIMUM VALUE OF THE CONTROL VECTOR 


C UT 


TRANSLATED FEEDBACK CONTROL VALUE 


C V 


MATRIX FOR FINAL POLE DETERMINATION USING EIGENVALUE 


C VLB 


VECTOR OF DESIGN VARIABLE LOWER BOUNDS 


C VUB 


VECTOR OF DESIGN VARIABLE UPPER BOUNDS 


C W1 


DISTURBANCE VALUE 


C N2 


REFERENCE SETTING 


C WA 


ADS WORK VECTOR 


C WL 


CIRCULAR FREQUENCY OF SINUSOIDAL DISTRUBANCE 


C WO 


COMPLEX VECTOR OF FINAL POLES 


C XI 


MOTOR SHAFT SPEED STATE VALUE 


C XIT 


TRANSLATION OF XI STATE 


C X2 


MOTOR CURRENT STATE VALUE 


C X2I 


INITIAL VALUE OF THE X2 STATE 


C X2MAX 


MAXIMUM VALUE OF THE X2 STATE 


C X2SP 


MAXIMUM SPECIFIED VALUE OF THE X2 STATE 


C Y 


VECTOR OF INTEGRATION FUNCTION VALUES 


C YOSCl 


VALUE OF SYSTEM OUTPUT OSCILATION AT TWO PREVIOUS INTEGRATIONS 


C YOSC2 


VALUE OF SYSTEM OUTPUT OSCILATION AT PREVIOUS INTEGRATION STEP 


C Y0SC3 


VALUE OF SYSTEM OUTPUT OSCILATION 



C YPRIME VECTOR OF INTEGRATION FUNCTION DERIVATIVES 
C Z EIGENVALUE DUMMY MATRIX 

C?t*3e6*********^**********9^*********?t*****************3H6*******^****^** 

c 

COMMON AA(3,3),BB(3),CC(3,3),H( 2),R( 2 ),Q( 2 ,2 ), F( 2 ,2 ) ,G( 3 ), 

$ W1,N2,U 

INTEGER IOG(5),IC( 20),INK(‘+00) 

REAL DV(3 3,VLB(3),VUB(3 ),CON(5),DF(3 3,WA( 20 ,20 3 ,WK( ^00 3 



131 



REAL LSTTIM,LSTOBJ,LSLOPE 



COMPLEX P(3),W0(3),2(3,3) 

INTEGER IGWK(^+) 

REAL FV1(3),GWK(69),V(3,3)>Y(^) 

EXTERNAL DUMMY, DERIVA 
OPEN(6,FILE='OUTPUT.DAT' ) 

OPEN( 7,FILE=‘C0N: ' ) 

C ******)t*H' INITIALIZATION SEGMENT ^6^*****^t*** 

C INITIALIZE THE AOS PARAMETERS 
NRA=20 
NC0LA=20 
NRHK=^00 
NRIWK=^00 
NDV=3 
N'':0N=5 
IGRAD=0 
INFO=-2 

C ESTABLISH TYPES OF CONSTRAINTS (ADS MANUAL TABLE 5) 

IDG(1)=-1 
IDG(2)=-1 
IDG(3)=-1 
IDG(^) = -1 
IDG(5)=2 

C LONER BOUND FOR DESIGN VARIABLES 
VUB(1)=1000, 

VUB( 2)=1000. 

VUB(3)=1000. 

C UPPER BOUND FOR DESIGN VARIABLES 
VLB(1)=-1000. 

VLB( 2)=-1000. 

VLB(3)=-1000. 

C ADS OPTIMIZATION METHOD SELECTION 
ISTRAT=1 
I0PT=3 
I0NED=3 
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IPRIMT=1000 



C^^****************^«H^*********************^^^**^*********')eHf******^**^(*^* 
C INITIALIZE THE INTEGRATION PARAMETERS 
FINTIM=1.0 
STEP=0.1 
TOL=l.E-3 

C SET CONSTRAINT SPECIFICATIONS 

C CI=INITIAL OUTPUT VALUE > CF=FINAL EXPECTED OUTPUT VALUE 
C CSP=MAXIMUM PERCENT OVERSHOOT FOR THE OUTPUT 
C TSP=MAXIMUM COMPLETION TIME 
C USP=MAXIMUM VALUE OF THE CONTROL VECTOR 
CI = 1. 

CF=0. 

CSP=100. 

TSP=1.5 

USP=10.0 

X2SP=10.0 

C SET INITIAL INPUT VALUES ( W1=DISTURBANCE , W2=REFERENCE SETTING) 

W1=0. 

W2=l. 

SEED=113. 

MEAN=1 

SD=0.1 

WL=157 

C SET INITIAL RUN FLAG AND EIGENVALUE PARAMETERS 
MATZ=0 
NM=3 
NN=3 
11=0 

C INITIAL SEGMENT ******** 

C SET ALL MATRIX AND VECTOR INPUTS 
P( 1)=CMPLX( -10.0,10.0) 

P( 2)=CMPLX( -10.0,-10.0) 

P(3)=CMPLX(-300.0,0. ) 
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AA(1,1)=0.E0 



AA(1,2)=S0.E0 
AA(1,3)=0.E0 
AA( C,1)=-200.S0 

aa( e,<: j=-eoo.EO 

AA(2,5)=0.E0 
AA(3,1)=-1.E0 
AA(3,2)=0.E0 
AA( 3,3 )=0.E0 
BB( 1)=0.E0 
BBC 2 )=200.E0 
BBC3 )=0.E0 
CCC 1,1)=-1.E0 
CCC 1,2 )=0.E0 
CCC 1,3^=O.EO 
CCC 2,1)=0.E0 
CCC 2,2)=0.E0 
CCC2,3)=0.E0 
CCC3,1)=0.E0 
CCC 3,2 )=0.E0 
CCC3,3)=0.E0 
HC 1)=0.E0 
HC 2)=1.E0 
FC 1,1)=-50.E0 
FC1,2)=0.E0 
FC 2,1)=0.E0 
FC2,2)=0.E0 
OC 1,1)=1.E0 
QC 1,2)=0.E0 
QC 2,1 )=0.E0 
QC2,2)=1.E0 
RC 1)=1.E0 
RC 2 )=0.E0 
DVC 



DVC2) = -1 



DV( 

DO 610 I=l,NCON 
610 COr4(Il=O.EO 
C POLE TO GAIN ROUTINE 

:ALL ADS( info, ISTRAT,I0PT,I0NED,IPRINT,IGRAD,NDV,NC0N,DV, VLB, 
$ VUB ,0BJ ,CON ,IDG ,NGT ,IC ,DF ,WA ,NRA ,NCOLA ,WK ,NRWK ,IWK ,NRIWK ) 

C TURN SCALING OFF 
IWK(2)=0 

660 CALL ADS( INFO, ISTRAT, lOPT ,IONED ,IPRINT,IGRAD, NDV,NCON,DV, VLB , 

$ VUB , OBJ ,CON ,IDG ,NGT ,IC , DF , WA ,NRA ,NCOLA ,WK ,NRNK , IWK ,NRIWK ) 
IF(INFO.EQ.0)GOTO650 
DO 700 1=1,3 
G(I)=DV(I) 

DO 710 J=l,3 

710 V( J ,I ) = AA( J ,I )-BB( J )*G( I ) 

700 CONTINUE 

C FIND EIGENVALUES FOR PROPOSED GAINS 

CALL EIGRF( V,NM,NN,MAT2,W0,Z,NM,FV1,IERR) 

RP1=REAL(P(D) 

RP2=REAL(P(2) ) 

RP3=REAL(P(3) ) 

RIP1=AIMAG(P( 1 ) ) 

RIP2=AIMAG(P(2)) 

RIP3=AIMAG(P(3) 1 
WR1=REAL(W0( 1) ) 

WR2=REAL(WO( 2) ) 

WR3=REAL(WO(3)) 

WI1=AIMAG(W0( 1)) 

WI2=AIHAG( WO( 2) ) 

WI3=AIMAG(WO(3)l 

IF( AMAXK WRl ,WR2 ,WR3 ) . GT . 0 . EO )G( 1 ) = -G( 1 ) 

CON( 1 )=ABS( AMAXK HRl ,WR2 ,WR3 )-AMAXl( RPl ,RP2 ,RP3 ) ) 

CON( 2 )=ABS( AMINK WR1,WR2,WR3 )-AMINl( RP1,RP2,RP3 ) ) 

CON( 3 )=ABSl AMAXK WI1,WI2 ,WI3 )-AMAXl( RIPl ,RIP2,RIP3 ) ) 

CON( A)=ABS(AMIN1(WI1,WI2,WI3)-AMIN1(RIP1,RIP2,RIP3) ) 
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OBJ=CON( 1 )+CON( 2 ) +CON( 2 )+CON( 4 ) 



GOTO 660 

C PRINTOUT INITIAL POLES AND GAINS 
650 WRITE(6,296)P(1),P(2),P(3) 

296 FORMAT(’0 DESIRED POLES: ‘ ,3( ' ( ' >E10.4, ' , ' ,E10.4, ' J) ’)) 

WRITE (6,297 )WR1 , WI 1 , WR2 , WI 2 , WRS ,HIS 

297 FORMAT('0 INITIAL POLES: ' ,S( ' ( ' ,E10 .4, * , ' >E10 .4, ' J) ')) 

WRITE(6,299)G(1),G(2),G(S) 

299 FORMAT! '0 INITIAL GAINS G( 1 ) ,G( 2 ) ,G( 2 ) = ',2E11.4) 

WRITE! 6,290) 

290 FORMAT! *0 STATE CONTROL MATRIX !B AB ) = ' ) 

WRITE! 6,298 )BB! 1 ) ,AA( 1 ,l)»fBB! 1 )+AA! 1,2 )>fBB( 2 ) 

WRITE! 6,298 )B3! 2 ) ,AA! 2,1 )*B8! 1 )+AA! 2 ,2 )>fBB! 2 ) 

298 FORMAT! 2E IS. 4) 

WRITE! 6,291) 

291 FORMAT! *0 OUTPUT CONTROL MATRIX ! CB CAB) = ') 

WRITE(6,298)CC!1,1)*BB!1)+CC! l,2)i^BB!2),CC!l,l)*!AA!l,l)»tBB!l) + 
$ AA( 1,2 )^<BB! 2) )+CC! 1,2 )*! AA! 2, 1 )^BB! 1 )+AA! 2,2 )^<BB( 2 ) ) 

WRITE! 6,292) 

292 FORMAT! '0 OBSERVABILITY MATRIX !C^ = ') 

WRITE! 6,298 )CC! 1,1 ) ,AA! 1,1 )^CC! 1 ,1 )+AA( 2,1 )^<CC( 1,2 ) 
WRITE(6,298)CC! 1,2),AA! 1,2)^CC! 1,1 )+AA!2,2))fCC(l,2) 

C CONTROLLER OPTIMIZATION ROUTINE 
C RESET OPTIMIZATION PARAMETERS 
IWK!2)=1 
ISTRAT=0 
I0PT=4 
IONED=7 
IPRINT=2020 
IDG!1)=0 
IDG!2)=0 
IDG!2)=0 
IDG(4)=2 

200 CALL A0S( INFO ,ISTRAT , lOPT ,IONED ,IPRINT , IGRAD ,NDV ,NCON ,DV ,VLB , 

$ VUB ,OBJ ,CON,IDG ,NGT ,IC ,DF ,WA ,NRA ,NCOLA ,WK ,NRWK ,IWK ,NRIWK ) 
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IF( INFO. EQ. 0 )WRITE( 6,300 ) 



300 FORMAT ( * TIME XI X2 U *, 

$ • OBJ C * ) 

C RESET CONSTRAINT AND RUN VALUES 
G(1)=DV(1) 

G(2)=DV(2) 

G(3)=DV(3) 

11 = 11+1 

JOB=0 

Y0SC1=CI 

YOSC2=CI 

YOSC3=CI 

CMAX=0.E0 

TAU=0,E0 

TC=0.E0 

U=0.E0 

UI=0.E0 

UMAX=0.E0 

LSTOBJ=O.EO 

LSL0PE=0.E0 

LSTTIM=0.E0 

C INITIALIZE THE DGEAR CONSTANTS FOR THE INTEGRATION TYPE DESIRED 
TIME=0.0E0 
TIMEND=0.E0 
DELT=0.001 
NDE=4 
METH=1 
MITER=0 
INDEX=1 

DO 201 1=1, NDE 
201 Y(I)=O.EO 

C INCREMENT THE INTEGRATION INTERVAL 
20 TIMEND=TIMEND+STEP 

WRITEt 7,455)11, TIME 

455 FORMAT ( * EVAL ^?*,I3,* TIME=‘,F5.3) 
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C CALL DIFFERENTIAL EQUATION SOLVER ROUTINE 



CALL DGEAR( NDE ,DERIVA , DUMMY , TIME ,DELT ,Y ,TIMEND ,TOL ,METH , 

$ MITER, inde:mgwk,gwk,ier) 

C UPDATE S“'ATE AND CONTROL VALUES 
Xl=V(i; 

X2=Y(2) 

0BJ=Y(4) 

C=CC( 1,1 )*X1+CC( 1,2 )*X2+H( 1 )*W1+H( 2 J*W2 
U=-G( 1 )*X1-G( 2 )*X2-G( 3 )*Y( 3 ) 

C DYNAMIC SEGMENT 9^***^*** 

C PRINT STATES AND OBJECTIVE - LAST RUN ONLY 

IF( INFO. EQ. 0 }WRITE( 6,310 JTIMEND ,X1 ,X2 ,U ,0BJ ,C 
310 F0RMAT(6G12.<+) 

C UPDATE CONSTRAINT CHECK VALUES 
Y0SC1=Y0SC2 
Y0SC2=Y0SC3 
Y0SC3=C 
H1=H2 
H2=H3 

H3=ABS(Y0SC3-CF) 

DELC=ABS(CI-C ) 

DELU=ABS(UI-U) 

C EVALUATE CONSTRAINTS 

I F ( DE LU . GE . UMAX ) UMAX =D E LU 
IF(TAU.NE.O.EO)GOTO 15 

IF( DELC .GE . ABS( 0. 63212E0*( CI-CF ) ) )TAU=TIMEND 
15 IF(DELC.GE.CMAX1CMAX=DELC 

IF(DELC.EQ.CMAX)CON( 1)=DELC-ABS(CI-CF )*( CSP/100+1 . EO ) 
IF(TC.NE.O.EO)GOTO 150 

IF( (H2.GE.HD.AND. (H2.GE.H3).AND.(H2.LE.ABS( ( CI-CF )* . OlEO ) )) 

$ TC=TIMEND 

IF( (H3.LT.H2).AND. (H2.LT,H1).AND.(H1.LT.ABS( (CI-CF)*.01E0) ) ) 

$ TC=TIMEND 

IF( (C.GT.(CF-ABS(CI-CF)*.01E0) ).AND.(TIMEND.GE .(S.EO^TAU) ).AND. 
$ (TAU.NE.O.EO).AND. (C.EQ.CMAX) )TC=TIMEND 
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CON( 2 ) = 10,E0 



IFdC.NE.O.EO )CON( 2 )=TC-TSP 
150 I F ( OE LU . EQ . UMAX )CON( 3 ) =0E LU-USP 
CCN(^)=-1.E0 

Ir( ASSJC ).GT.ABS(CI-CF )*10.E0*( CSP/100 . EO + 1 . EO ) )CON( ) = 1 . EO 

I F ( CON( ) . GT . 0 . EO ) JOB= 1 

CON(5)=-1.EO 

IF( (ABS(CF-C).GT.ABS(CF-CI ) ) . ANO . ( TAU . EQ . 0 . EO ) )CON( 5 ) = 1 . EO 
IF(CON(S).GT.O.EO )J0B=1 
C EVALUATE FOR ENO OF RUN 

IF(TIMEND.LE.STEP*3. )GOTO ^0 
SLOPE=ABS( ( OBJ-LSTOBJ )/( TIMEND-LSTTIM ) ) 

OSLOPE=A8St ( SLOPE-LSLOPE )/( TIMENO-LSTTIM ) ) 

IF{ ( (SLOPE.lt. .OOIEO ).OR.(OSLOPE.LE. .OOIEO) ).AN0.(TC.NE.0.E0) ) 

$ J0B=1 
<hO LSTOBJ=OBJ 

LSTTIM=TIMEN0 

LSLOPE=SLOPE 

C CHECK FOR END OF RUN CONOITION 
IF( JOB.EQ.DGOTO 30 
C CHECK FOR MAX TIME 

IF(TIMENO.GE.FINTIM)GOTO 30 
GOTO 20 

C ******** TERMINAL SEGMENT ******** 

C CHECK FOR ENO OF JOB 

30 IF(INFO.NE.O)GOTO 200 
00 400 1=1,3 
00 410 J=l,3 

410 V(I,J)=AA(I,J)-;BB( I)*G( J) 

400 CONTINUE 
C CALCULATE FINAL POLES 

CALL EIGRF(V,NM,NN,MATZ,H0,Z,MM,FV1,IERR) 

WRITE ( 6,497 )W0( 1 ) ,W0( 2 ) ,W0( 3 ) 

WRITE ( 6 ,499 )G( 1 ) ,G( 2 ) ,G( 3 ) 

497 FORMAT( *0 FINAL POLES: ’ ,3( M ' ,E10 .4, * , * ,E10 .4, ' J) ’)) 
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499 FORMAK *0 FINAL GAINS G( 1 ) ,G( 2 ) ,G( 3 ) 



* ,3E11.4) 



STOP 

END 

C 9f**^**** DERIVATIVE SEGMENT ******** 

SUBROUTINE OERIVA( NDE ,TIME , Y ,YPRIME ) 

COMMON AA(3,3),B8(3),CC(3,3),H(2),R( 2),Q(2,2),F(2,2 ),G(3), 
$H1,H2,U 

REAL Y(NDE),YPRIME(MDE) 

X1=Y(1) 

X2=Y(2) 

C DIFFERENTIAL EQUATIONS FOR CONTROLLER 

X1D0T=AA( 1,1 )*X1+AA( 1,2 )*X2+BB( 1 )*U+F( 1 , 1 )*W1 + F( 1 ,2 )*W2 
X2D0T=AA( 2,1 )^X1+AA( 2,2 )*X2+BB( 2 )*U+F( 2 ,1 )^W1 + F( 2 ,2 )*W2 
Y3D0T=CC( 1 ,1 )^X1+CC( 1, 2 )*X2+H( 1 )*W1+H( 2 )*W2 
U=-G( 1 J*X1-G( 2 )*X2-G( 3 )*Y( 3 ) 

UT=U-1.E0 

X1T=X1-1.E0 

OBJDOT = ( X1T*Q( 1,1 )+X2*Q( 1,2) )*XlT+( X1T*Q( 2,1) +X2*Q( 2,2) )*X2 + 

$ (UT*R(1)*UT) 

YPRIME(1)=X1D0T 

YPRIME( 2)=X2D0T 

YPRIME(3)=Y3D0T 

YPRIME(4)=0BJD0T 

RETURN 

END 

C DUMMY ROUTINE REQUIRED BY DGEAR 

SUBROUTINE DUMMY (N,X,Y,PD) 

INTEGER N 

REAL Y(N),PD(N,N),X 

RETURN 

END 
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ADSLPC BATCH FILE 



9 



echo off 

rem: ASK Is a orosram to query the user for a yes or no response 
rem; KEDIT Is the line editor prosram of choice 
rem; PUTIME displays the current clock time 

rem: PU7IME3 displays the current clock time and sounds the system bell 

rem: BLANK. DATls an empty file 

rem: ScTPRN sends escapecodes to the printer 

echo off 

path adsl util; adsl prosram: fortran; : utlltys 

cd adsl prosram 

;start 

ask Edit a Compile ADSLPC. FDR? 

If errorlevel 1 soto start2 
; b e s 1 0 

del ADSLPC. DBJ 
; star t5 

KEDIT ADSLPC. FDR 
putlme 

FORI ADSLPC. ADSLPC. NULrNUL 

putlme 

PAS2 

If exist ADSLPC. DBJ soto link 
putlmeb 

ask Re-edit a Compile ADSLPC. FDR? 
if errorlevel 1 soto start2 
goto Starts 
;start2 

ask Link and Run the ADSLPC prosram? 

If errorlevel 1 soto llnk2 

:llnk 

putlme 

LINK ADSLPC*ADS1A<^ADS1B*ADS2A*ADS2B*ADS3A*ADS3B*ADS3C«'IMSL1 + IMSL2^IMSL3.MDDEL87»NUL. FDRTRAN 8087. LIB* FORTRAN FDRTRAN.LIB 

putlme 

soto runl2 

:llnk2 

ask Run the ADSLPC model? 

If errorlevel 1 soto print 
:r unl2 

copy blank.dat outPUt.dat 

putlme 

MDCELB7 

putlme 

pause 

BRDNSE output.dat 
: print 

ask Printout results (If yes — verify printer on. paper allsned)? 

If errorlevel 1 soto donel 

setprn [271N(B1 

setprn (151 

PRINT OUTPUT.DAT 

ASK PRINT ADSLPC. FDR? 

IF ERRDRLEVEL=1 GOTD NEKT2 
PRINT ADSLPC. FDR 
:NEXT2 
: done 1 

ask Edit a Compile ADSLPC. FDR? 

If errorlevel 1 soto done 
soto besln 
: done 
echo on 
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